Spectroscopic survey of Kepler stars

Spectroscopic survey of Kepler stars.1 I. HERMES/Mercator observations of A- and F-type stars

Abstract

The Kepler space mission provided near-continuous and high-precision photometry of about 207,000 stars, which can be used for asteroseismology. However, for successful seismic modelling it is equally important to have accurate stellar physical parameters. Therefore, supplementary ground-based data are needed. We report the results of the analysis of high-resolution spectroscopic data of A- and F-type stars from the Kepler field, which were obtained with the HERMES spectrograph on the Mercator telescope. We determined spectral types, atmospheric parameters and chemical abundances for a sample of 117 stars. Hydrogen Balmer, Fe i, and Fe ii lines were used to derive effective temperatures, surface gravities, and microturbulent velocities. We determined chemical abundances and projected rotational velocities using a spectrum synthesis technique. The atmospheric parameters obtained were compared with those from the Kepler Input Catalogue (KIC), confirming that the KIC effective temperatures are underestimated for A stars. Effective temperatures calculated by spectral energy distribution fitting are in good agreement with those determined from the spectral line analysis. The analysed sample comprises stars with approximately solar chemical abundances, as well as chemically peculiar stars of the Am, Ap, and  Boo types. The distribution of the projected rotational velocity, , is typical for A and F stars and ranges from to about  km s, with a mean of  km s.

keywords:
stars: general – stars: abundances – stars: chemically peculiar – stars: rotation – space missions Kepler
24

1 Introduction

The NASA space mission Kepler was successfully launched on 7th March, 2009. The main scientific goal of the mission was to discover Earth-sized planets around Sun-like stars (Borucki et al., 2010). During May 2013 the satellite lost the second of its four reaction wheels (the first one failed in July 2012), preventing further observation of the original field. Nevertheless, during the four years of near-continuous data collection, Kepler produced photometric time series of an exceptional precision of a few ppm (parts per million), which provide a unique and excellent database for studying variability due to pulsations in thousands of stars across the Hertzsprung-Russell (H-R) diagram (Kjeldsen et al., 2010).

Asteroseismology is the only way to probe the internal structure of stars and significantly improve our understanding of stellar evolution. It provides direct tests of the modelling of processes taking place in stellar interiors, such as diffusion, transport of angular momentum (Mosser et al., 2012; Marques et al., 2013), convective overshooting, mechanisms driving pulsations (e.g. Moravveji et al., 2012; Antoci et al., 2014), and the value of the mixing-length parameter in stars of different effective temperatures (Trampedach & Stein, 2011; Bonaca et al., 2012). The key ingredients for an asteroseismic study are precise pulsation frequencies, mode identification, and strong constraints on atmospheric parameters, such as effective temperature , surface gravity , chemical abundances, and rotational velocity . Accurate values of the pulsation frequencies are provided by the ultra-precise Kepler photometry. The accuracy of atmospheric parameters available from the Kepler Input Catalogue (KIC) (Brown et al., 2011) is generally too low for asteroseismic modelling, particularly for stars of spectral types earlier than F0 (Huber et al., 2014). Also, information on the stellar chemical composition and rotational velocity is lacking. Therefore, high-resolution spectra which offer a wide range of information are needed to fully characterise Kepler targets (Uytterhoeven, 2014).

For more than a thousand Kepler stars high- and medium-resolution spectra have been collected (Uytterhoeven et al., 2010). These observations are summarized in the Kepler Asteroseismic Science Operations Center (KASOC) database5. In this paper, we focus on the analysis of high-resolution spectra of main-sequence A- and F-type stars. The main groups of pulsators among these objects are  Scuti,  Doradus, and  Dor/ Sct hybrid stars. The  Sct stars are main-sequence or evolved stars situated in the classical instability strip. They pulsate in low-order radial and non-radial pressure modes (p-modes) with typical frequencies of 5 to 50 d (see reviews by Breger, 2000; Handler, 2009; Balona, 2014). The  Dor stars are pulsating variables with spectral types A7–F5 that are on or close to the main-sequence and pulsate in high-order gravity modes (g-modes) with typical frequencies of 0.3 to 3 d (Kaye et al., 1999). The instability strip for  Dor stars partly overlaps with the  Sct instability strip. The  Dor variables are of particular interest in asteroseismology, as their g-mode pulsations probe the deep interiors of these stars, reaching their convective cores. Additionally, these pulsators occupy a region in the H-R diagram where  Dor-like g-mode pulsations,  Sct-like p-mode pulsations and possibly solar-like oscillations can co-exist in a single star. This overlap means that some individual objects show g-mode as well as p-mode pulsations. These hybrid variables permit simultaneous probing of two distinct regions of the interior: the deep layers between the convective core and the envelope, where the g-modes propagate, and the outer layers of the envelope where the p-modes propagate (Kurtz et al., 2014; Saio et al., 2015). Therefore, hybrid stars are key objects for testing stellar models with seismic techniques.

Grigahcène et al. (2010) found that the hybrid stars in the Kepler field are spread over the entire  Sct and  Dor instability domains. The analysis of Kepler data and KIC parameters for 750 A and F stars presented by Uytterhoeven et al. (2011) confirms this result. They found that 23% of their sample of A and F objects shows both  Sct pulsations and low-frequency variations. Balona & Dziembowski (2011) analysed observations of 1568  Sct pulsators taken from the Kepler archive and identified a group of  Sct stars with low-frequency peaks similar to those characteristic for  Dor stars. Contrary to the results of Uytterhoeven et al. (2011) and Grigahcène et al. (2010), these variables are located in the same region of the H-R diagram as the ground-based  Dor pulsators and none is found in the  Sct instability domain. The low frequencies are observed in the hotter stars, but their origin is not clear (Balona, 2014). A high-resolution spectroscopic analysis by Tkachenko et al. (2013b) found each of their 38  Dor stars to lie within the  Dor instability strip. Accurate stellar atmospheric parameters are clearly required to place  Dor/ Sct hybrids unambiguously in the H-R diagram.

This paper is organised as follows: in Sect. 2, we describe observational data used in the analysis. Details on the spectral classification methods are given in Sect. 3. In Sect. 4, we present all the methods applied for determining the atmospheric parameters of the investigated objects, the sources of possible errors, and discuss consistency of the results obtained from different approaches. The distribution of rotational velocities is presented and discussed in Sect. 5. The derived abundances of elements, their correlation with the other parameters and a sample of chemically peculiar (CP) stars are discussed in Sect. 6. Positions of the analysed stars in the  –  diagram, conclusions and future prospects are given in Sect. 7.

2 Spectroscopic surveys of Kepler stars

High-resolution spectra were collected for several hundred stars located in the Kepler field-of-view. Most of these spectra were obtained by échelle spectrographs attached to the Mercator Telescope6, Nordic Optical Telescope7 (NOT), Tautenburg 2-m Alfred-Jensch Telescope8 (TLS), Bernard Lyot Telescope9 and the Canada-France-Hawaii Telescope10 (CFHT).

In this paper we analyse high-resolution spectra taken with the cross-dispersed, fibre-fed échelle spectrograph HERMES (High Efficiency and Resolution Mercator Échelle Spectrograph; Raskin et al., 2011) attached to the 1.2-m Mercator Telescope located on La Palma (Canary Islands, Spain). The spectra have a resolving power and cover the spectral range from 3770 to 9000 Å. The typical signal-to-noise (S/N) ratio for an individual spectrum at 5500 Å is 80–100. The spectra have been reduced with a dedicated pipeline11, which includes bias subtraction, extraction of scattered light produced by the optical system, cosmic ray filtering, division by a normalized flat-field, wavelength calibration by a ThArNe lamp and order merging. The normalisation to the continuum was performed manually by using the standard iraf12 procedure continuum. Most of the stars in our sample have    km s. The spectra of rapid rotators are very difficult to normalize correctly in spectral regions rich in blended lines, e.g. for  Å. For this reason, we checked the quality of our normalisation by fitting synthetic spectra to all Balmer lines in the spectrum of each star. If the spectrum has been normalized correctly, all hydrogen lines can be fitted satisfactorily with the same . The process of normalization is otherwise repeated. For the few stars for which we have more than one spectrum, we analysed the averaged spectrum. The averaging process is applied after normalization and forms an additional test of the quality of the normalization process.

The data analysed in this paper were collected in spectroscopic surveys carried out in the years 2010 to 2012 (Table 1). Table 2 lists the analysed stars and gives their KIC numbers, dates of the observations, number of spectra, visual () magnitudes, and spectral types from the literature and determined in this paper. We also provide notes on chemical peculiarity, binarity and pulsation characteristics. Chemically peculiar stars in our sample were identified in the process of spectroscopic classification and confirmed by the detailed spectroscopic analysis afterwards. Twenty stars were found to be double-lined spectroscopic binaries (SB2) and the atmospheric parameters of these stars will be presented in a separate paper (Catanzaro et. al, in preparation).

The results of the spectroscopic analysis of stars observed with the other telescopes will be presented in a forthcoming paper (Niemczura et al., in preparation).

Observation Primary Observers
dates Investigator
2010 May-Sept. M. Briquet KU Leuven observers
2011 July E. Niemczura D. Drobek
2011 Sept. E. Niemczura J. F. S. Gameiro
2012 July K. Uytterhoeven E. Niemczura
2012 August P. De Cat P. I. Pápics
2012 Sept. P. De Cat P. M. Arenal

E. Van Aarle, S. Bloemen, N. Gorlova, P. Lampens, H. Van Winckel, B. Vandenbussche

Table 1: HERMES spectroscopic surveys of Kepler stars.
KIC HERMES N SpT1 SpT2 Notes
Number observations [mag] [this paper]
1294756 2010 May-September 1 8.96 A2 A3 IV mult. + rot. +  Sct
2571868 2011 July 1 8.67 A0 A3 IVn hybrid (low. ampl.  Dor)
2694337 2010 May-September 1 10.35 F0.5 IVs hybrid
2695344 2012 July 1 9.62 F2 F2 V SB2;  Dor
2012 September 1
3230227 2010 May-September 1 8.94 A5 A5 IV: SB2; ecl. bin. +  Sct
3231985 2011 September 1 9.22 A2 A4 IV Ca weak (A3) rot.
3347643 2010 May-September 1 7.69 A2 A3 Van hybrid (low. ampl.  Dor)
3441230 2012 July 1 9.79 A2 kA2hA5mF0 V Am rot.
3656913 2012 July 1 9.81 A0 kA2hA6mF0 IV Am SB2; bin.
3850810 2010 May-September 1 10.10 F1 Vs wk met (F0) hybrid + rot.

KIC 3230227: Incorrectly normalised H lines; wings are unusable but core matches A5 IV. Metal lines indicate A5 V.

KIC 3231985: Weak Ca K and Ca 4226 line.

KIC 3850810: Considerably metal weak or poorly normalised.

Heckmann & Dieckvoss (1975)

Table 2: Journal of spectroscopic observations and the derived spectral classes. is the number of available spectra; SpT1 indicates the spectral classification from the literature, whereas SpT2 means spectral classifications obtained in this work. Superscript in the first column denotes the stars with additional information below the table. Additional notations for luminosity class: “b” – lower luminosity main-sequence star; “a” or “a+” – higher luminosity main-sequence star; “IV-V” – between IV and V; “IV/V” – either IV or V; “n” – nebulous (i.e. broad-lined); “s” – sharp lines; “nn” – very rapid rotators; “wk” met – weak metal lines; “met str” – strong metal lines. In the last column additional information taken from spectra and Kepler data are presented. “SB2” means double-lined spectroscopic binary; “rot.” – rotational frequency detected; “bin.” – orbital frequency detected; “mult.” – orbital frequencies of both components detected; “hybrid” – frequencies of  Dor and  Sct visible in the spectra; “contamination” – Kepler data contaminated by a nearby star; “?” – our classification of variability type is uncertain. The full table is available in the electronic form.

3 Spectral Classification

Stellar spectral classification in the Morgan-Keenan (MK) system involves determining the luminosity class and spectral type for a given star by comparing its spectrum with spectra of standards (Gray & Corbally, 2009). The spectral type in A stars is determined mainly from three spectral features: (i) hydrogen lines, (ii) the Ca ii K line, and (iii) lines of metals. In a non-peculiar star, all three features should consistently provide the same spectral type, if the correct luminosity class has been adopted. Determination of luminosity class in A stars is difficult. While the hydrogen lines are very sensitive to luminosity in early A stars, the A5 and later spectral types lose much of this sensitivity and at F0 the hydrogen lines become almost insensitive to the luminosity class. For F0, ionised lines of Fe and Ti become useful in determination of the luminosity class.

At late A spectral type, metallic (Am) stars are common. For these stars spectral types determined from the Ca ii K-line and metal lines differ by five or more subclasses. The Ca ii K line is weak in their spectra and thus yields an earlier spectral type than hydrogen lines, while the metal lines are enhanced and typical for a later subtype. They are denoted as, for example, kA2hA5mF0 (IV), where k, h and m precedes the spectral type obtained from the Ca ii K line, hydrogen lines and metal lines, respectively. These symbols are followed by the luminosity class. If the luminosity class is in parentheses, as it is in this example, then the star exhibits the anomalous luminosity effect (ALE, see Gray & Corbally 2009) where the violet part of the spectrum indicates a more luminous class than the blue part. Stars showing mild Am features, where there are fewer than five spectral subtypes between the Ca ii K line and metal line types, are called marginal Am stars, denoted as “Am:”.

Most of the objects were classified as A stars, including slow, fast and very fast rotators (“s”, “n” and “nn” added to the spectral type, respectively). The classified CP stars are of Am (17), Ap CrSrEu (1) and  Bootis types (2). Additional, mild peculiarities were discovered in some of the stars, such as weak Ca or metal lines. Only three stars were classified as late B, with one of them appearing to be He strong (He-s). The assigned spectral classes are presented in Table 2 and some notes are provided in the comments below the table.

4 Atmospheric parameters

To perform an abundance analysis, one needs an appropriate atmospheric model of the star, which requires the effective temperature, surface gravity and metallicity to be known. To determine these parameters, we applied both photometric and spectroscopic methods.

KIC  KIC  Na D & SED  Metal & Balmer lines
Number
[K] [mag] [mag] [K] [K] [km s] [km s]
10263800
10264728
10355055
10533616
10549371
10555142
10590857
10721930
10977859
11013201
Table 3: Atmospheric parameters of the investigated Kepler A and F stars. We present parameters derived from the KIC photometry (KIC), those obtained from spectral energy distributions (Na D & SED) and from analysis of metal and Balmer lines. The uncertainties are given for effective temperatures derived from SED fitting, and for  and values derived from metal lines. The full table is available in electronic form.

4.1 Effective temperatures from Spectral Energy Distributions

Effective temperature can be determined from the stellar spectral energy distribution (SED). For our target stars these were constructed from the available photometry, using 2MASS (Skrutskie et al., 2006), Tycho and magnitudes (Høg et al., 1997), USNO-B1 magnitudes (Monet et al., 2003), and TASS magnitudes (Droege et al., 2006), supplemented, if available, by Geneva photometry (Rufener, 1999), CMC14 magnitudes (Evans et al., 2002) and TD-1 ultraviolet flux measurements (Carnochan, 1979).

SEDs can be significantly affected by interstellar reddening. We have, therefore, estimated this parameter from the interstellar Na D (5889.95 Å) lines. The equivalent widths of the D lines were measured and values were obtained using the relation given by Munari & Zwitter (1997). For resolved multi-component interstellar Na D lines, the equivalent widths of the individual components were measured. The total in these cases is the sum of the reddening per component, since interstellar reddening is additive (Munari & Zwitter, 1997). The SEDs were de-reddened using the analytical extinction fits of Seaton (1979) for the ultraviolet and Howarth (1983) for the optical and infrared.

Effective temperatures were determined by fitting Kurucz (1993a) model fluxes to the de-reddened SEDs. The model fluxes were convolved with photometric filter response functions. A weighted Levenberg-Marquardt non-linear least-squares fitting procedure was used to find the solution that minimized the difference between the observed and model fluxes. Since  and metallicity [M/H] are poorly constrained by our SEDs, we fixed and [M/H]  dex for all of the fits. The results are given in Table 3. The uncertainties in  include the formal least-squares error and adopted uncertainties in of 0.02 mag, of 0.5 dex and [M/H] of  dex added in quadrature.

4.2 Atmospheric parameters from spectroscopy

We used atmospheric models and synthetic spectra to obtain atmospheric parameters, abundances of chemical elements and rotational velocities from the high-resolution HERMES spectra. All the necessary atmospheric models were computed with the line-blanketed, local thermodynamical equilibrium (LTE) ATLAS9 code (Kurucz, 1993a). The physics of the models includes plane-parallel geometry, and hydrostatic and radiative equilibrium. Models were pre-calculated in a grid for effective temperatures between and  K with a step of  K, surface gravities from to  dex with a step of  dex, microturbulence velocities between and  km s with a step of  km s, and for solar metallicity. Additional models of non-solar metallicity were calculated for smaller ranges of effective temperatures and surface gravities. The synthetic spectra were computed with the SYNTHE code (Kurucz, 1993b). Both codes, ATLAS9 and SYNTHE, were ported to GNU/Linux by Sbordone (2005) and are available online13. The stellar line identification and the abundance analysis in the entire observed spectral range were performed on the basis of the line list from Castelli & Hubrig (2004).

We analysed the metal lines following the spectrum synthesis methodology established by Niemczura & Połubek (2006) and Niemczura et al. (2009). It relies on an efficient least-squares optimisation algorithm (see Bevington, 1969; Takeda, 1995), and is appropriate over a wide range of  values. Our observed range is to  km s, which includes stars rotating too rapidly for equivalent width analysis to be possible due to line blending. The spectrum synthesis method allows for a simultaneous determination of various parameters influencing stellar spectra and consists of the minimisation of the deviation between the theoretical and observed spectra. The synthetic spectrum depends on stellar parameters such as , , , , radial velocity , and the relative abundances of the elements , where denotes the individual element. All atmospheric parameters are correlated. In our method, the , and  parameters were obtained prior to the determination of abundances of chemical elements and were considered as input parameters. The remaining parameters (, and ) were determined simultaneously because they produce detectable and different spectral signatures. The  values were determined by comparing the shapes of observed metal line profiles with the computed profiles, as shown by Gray (2005).

The starting values of  were taken from the SED fits, and values  came from the KIC. These parameters were then improved by using the sensitivity of Balmer lines to effective temperature and gravity, following the method proposed by Catanzaro et al. (2004). However, hydrogen lines cease to be good indicators of  for effective temperatures cooler than  K. Therefore, for stars with lower effective temperatures, the surface gravity was initially assumed to be equal to  dex. We used an iterative approach to minimize the differences between observed and synthetic H, H and H profiles (see Catanzaro et al., 2004).

The effective temperatures and surface gravities determined from Balmer lines were improved through the analysis of lines of neutral and ionised iron. In this method, we adjusted , and  by comparing the abundances determined from Fe i and Fe ii lines. The analysis was based on iron lines because they are the most numerous in our spectra. In general, we required that the abundances measured from Fe i and Fe ii lines yield the same result. In the spectra of stars with lower effective temperatures, the number of Fe i lines outnumbers those of the Fe ii lines. The strength of Fe i lines depends on , and metallicity [M/H], but is practically independent of . On the other hand, the lines of Fe ii are weakly sensitive to effective temperature and metallicity, but highly sensitive to gravity (Gray, 2005). First, we adjusted microturbulence until we saw no correlation between iron abundances and line depths for the Fe i lines. Next,  was changed until we saw no trend in the abundance versus excitation potential for the Fe i lines. Determinations of microturbulence and effective temperature are not independent, but since the influence of microturbulence on spectral lines is stronger, we adjusted this parameter first. Then, surface gravity is obtained by fitting the Fe ii and Fe i lines, and by requiring the same abundances from the lines of both ions. For hotter stars, whose spectra show more Fe ii than Fe i lines, the Fe ii lines were used to determine  and , whereas  was derived from Balmer lines (Gray, 2005). This approach is limited by stellar rotation, since line blending grows rapidly with . In the spectra of rapid rotators all the lines are blended and it is impossible to use individual iron lines for ,  and  determination. In such cases we assumed that the abundances obtained from Fe i and Fe ii lines are the same for the correct , and . The analysis relies on the same physics, regardless of rotation. The only difference is that for slow and moderate rotators we can use separate iron lines, including weak lines, to calculate the atmospheric parameters. For rapid rotators we are comparing iron abundances derived from broader spectral regions. Typically such regions contain many iron lines, both neutral and ionised. Of course, in such cases, stronger lines have more influence on the result.

Having derived ,  and , the determination of abundances could be performed. The analysis of chemical abundances began with the selection of spectral regions suitable for abundance analysis. The length of a chosen part depended mainly on . For slowly rotating stars ( km s), short sections covering only one or a few blended spectral features were used. For stars with  km s we used broader spectral ranges re-normalised by comparison with theoretical spectra, if necessary. The next step was line identification in the chosen region. In the case of rapidly rotating stars, where all lines are blended, only the most important elements producing spectral features were considered. In our method we can take into account all elements which have lines in the chosen spectral region. Elements that have little or no influence on the spectra are included in the analysis, but no attempt is made to calculate their abundances. The stellar rotation and signal-to-noise ratio of the spectrum dictate for which elements the abundance analysis is performed.

The line list used is that of Castelli & Hubrig (2004), which is based on Kurucz’s line list14 and supplemented by additional sources. It was prepared for the analysis of chemically peculiar Ap and HgMn stars, hence is suitable for the analysis of our CP stars. Each selected spectral region was analysed by the spectral synthesis method described above. We adjusted ,  and chemical abundances to minimize the difference between the calculated and observed spectra. The final solution was typically reached after ten iterations. The iterations stopped if the values of the determined parameters remained the same within 2% for three consecutive steps. In our method, chemical abundances and  values were determined from different parts of the spectrum. In the final step of the analysis, we derived the average values of  and abundances of all the chemical elements considered for a given star and estimated the uncertainty on the derived abundances.

In our analysis we do not take into account macroturbulence, whose influence on the line profile is degenerate with . The effect of macroturbulence on F stars was discussed in detail by Fossati et al. (2011), who found that stronger lines required higher macroturbulent velocities and interpreted this as the effect of unmodelled depth-dependent velocities in the atmosphere. Since these are not observable in A stars, the required macroturbulent velocities remain uncertain, but the effect of macroturbulence is mainly limited to . The other parameters can be slightly affected, especially in case of slowly rotating stars (Fossati et al., 2011). The  values of slow rotators may be overestimated due to non-inclusion of macroturbulence. However, for  km s the effect of modest (e.g.  km s) macroturbulence is negligible at our resolution and signal-to-noise.

The atmospheric parameters obtained from our analysis are given in Table 3, and the abundances of analysed elements are given in Table 4.

KIC 10263800 KIC 10264728 KIC 10355055 KIC 10533616 KIC 10549371 KIC 10555142 KIC 10590857
C 8.850.14 (1) 8.480.12 (2) 8.900.12 (2) 8.380.07 (7) 8.310.14 (14) 8.780.16 (5) 8.360.19 (5)
N 8.330.14 (1)
O 8.740.14 (2) 9.040.12 (1) 8.440.12 (1) 8.920.14 (2) 8.830.14 (1) 8.910.17 (1)
Ne
Na 6.600.12 (1) 7.100.12 (1) 6.890.14 (2) 6.460.14 (2) 6.630.18 (1) 6.220.17 (1)
Mg 7.680.14 (1) 8.010.12 (4) 7.280.19 (3) 7.980.10 (6) 7.560.13 (8) 7.640.18 (2) 7.570.10 (4)
Al
Si 6.800.14 (2) 7.940.12 (2) 7.500.12 (2) 7.600.18 (5) 7.540.20 (21) 7.520.29 (5) 7.560.25 (8)
P
S 7.320.01 (3) 7.290.17 (1)
Cl
K
Ca 6.620.14 (2) 6.640.12 (4) 6.610.16 (8) 6.650.16 (15) 6.540.16 (25) 6.300.23 (7) 6.670.23 (4)
Sc 4.030.14 (1) 3.090.12 (3) 2.960.11 (3) 3.100.16 (5) 3.240.10 (9) 2.910.21 (4) 3.280.13 (6)
Ti 4.430.14 (3) 5.070.16 (8) 5.050.18 (9) 5.220.17 (16) 5.060.15 (47) 5.050.15 (8) 5.270.20 (14)
V 2.790.12 (1) 4.430.11 (9) 4.670.18 (2)
Cr 5.460.14 (2) 5.690.19 (5) 5.600.04 (6) 5.860.09 (8) 5.530.16 (26) 5.690.18 (2) 5.670.18 (16)
Mn 5.600.12 (2) 5.290.12 (1) 5.570.14 (2) 5.350.16 (10) 5.430.04 (3) 5.590.17 (2)
Fe 6.930.15 (5) 7.540.12 (11) 7.290.11 (17) 7.700.11 (36) 7.320.09 (101) 7.450.14 (20) 7.540.09 (36)
Co 5.880.14 (2)
Ni 6.360.07 (4) 6.110.11 (6) 6.350.23 (11) 6.190.18 (41) 6.270.13 (6) 6.410.16 (13)
Cu 4.410.14 (1)
Zn 4.780.14 (2)
Ga
Sr 2.550.12 (1) 1.670.12 (1) 2.960.14 (2) 3.230.14 (1) 1.360.18 (1)
Y 3.140.12 (1) 2.240.12 (1) 2.820.14 (1) 2.510.17 (10) 3.030.29 (3) 2.850.17 (1)
Zr 3.180.20 (10) 2.980.18 (1)
Ba 3.780.12 (2) 3.110.12 (1) 2.480.20 (3) 2.750.13 (3) 2.690.18 (2) 2.430.17 (2)
La 1.600.14 (2)
Ce 1.690.14 (2)
Pr
Nd 1.720.20 (5)
Table 4: Abundances of chemical elements for a sample of analysed stars (on the scale in which (H) = 12). The median value is given; standard deviations were calculated only if the number of analysed parts of a spectrum was greater than three. In other cases the average value calculated from standard deviations of other elements are given. Number of analysed parts is given in brackets. The full table is available in the electronic version.

4.3 Discussion of uncertainties

The derived atmospheric parameters and chemical abundances are influenced by errors from a number of sources. First of all, the adopted atmosphere models and synthetic spectra were calculated taking into account several simplifications, e.g. LTE, plane parallel geometry (1-D), and hydrostatic equilibrium (Kurucz, 1993a), instead of a 3-D and non-LTE approach (e.g. Bergemann et al., 2012).

Recently, the non-LTE line formation for Fe i and Fe ii for stars with  between  K and  K and   between and  dex was discussed by Mashonkina (2011). In this paper the non-LTE calculations were performed with a complete model atom, as described therein. The departures from LTE caused the depletion of total absorption in the Fe i lines and positive abundance corrections, leading to an underestimation of Fe abundance in LTE calculations. However, the corrections to LTE abundances were found to be smaller than  dex for main-sequence stars. For Fe ii, the non-LTE effect was again found to be negligible, with corrections of to  dex over the whole range of considered stellar parameters. Most stars in our sample have effective temperatures lower than  K. The non-LTE corrections of Fe abundances for these stars less than  dex. Application of this correction would change  determinations from the ionisation balance of Fe i and Fe ii lines. If we add non-LTE corrections to the obtained abundances of Fe i, the determined surface gravity would be higher by about  dex. However, in most of our investigated stars, high projected rotational velocities prevents us from obtaining abundances of Fe i and Fe ii separately. Therefore, we did not include the non-LTE corrections in the present analysis.

Other important factors affecting spectral analysis are atomic data (Kurucz, 2011), quality (resolution, S/N), wavelength range of the observed spectra and their normalisation. It is quite easy to normalise the optical spectrum of an A star with   km s. The difficulty increases for heavily blended spectra of stars with  much higher than this value, even for high-resolution and high S/N data.

All the above-mentioned factors influenced our derived atmospheric parameters and chemical abundances. The chemical abundances are also affected by inaccurate atmospheric parameters. The influence of errors in atmospheric parameters on the derived chemical abundances can be checked by perturbing the synthetic spectrum by the adopted step size in the grid of atmospheric models. For example,  K leads to changes in abundances smaller than  dex in most cases. Similarly, changes chemical abundance by about  dex, and  km s changes abundances by  dex at most. It should be noted here that the uncertainty on the microturbulence velocity increases with , from  km s for small and moderate velocities, to  km s for rapid rotators. Chemical abundances determined for rapidly rotating stars are affected by the greatest uncertainty. The combined errors of chemical abundances were calculated as:

In most cases was less than  dex.

4.4 Comparison of results from different methods

Figure 1: Distributions of (left panel) and (right panel). Gaussian fits are shown by solid lines. Central values, arithmetic means and standard deviations of both distributions are indicated in figures.
Figure 2: Left panel: differences as a function of (spec). A linear function is fitted. The dotted horizontal line indicates a difference of mag. Right panel: the distribution of differences . Arithmetic mean and standard deviation of the distribution are indicated in the figure.

Interstellar reddening

The values from the KIC, , and those determined here from the equivalent widths of sodium D lines, , are given in Table 3. The distributions of these parameters are presented in Fig. 1. It is clear that the photometric interstellar reddening given in the KIC is higher than . The average value of is equal to  mag, whereas from sodium lines it is  mag. The reddenings obtained from Na D lines are lower than  mag for all stars, except for KIC 9552758 (HD 190293) and KIC 8915335 (HD 190566), for which they equal  mag. We can see from Fig. 2 that the differences between interstellar reddenings obtained from Na D lines and those listed in KIC slightly increase with effective temperature. As explained by Brown et al. (2011), the values of given in the KIC were approximated using a simple model of dust distribution and with an adopted disk scale-height of  pc – larger than suggested by earlier estimates (e.g. Koppen & Vergely, 1998; Marshall et al., 2006). Brown et al. (2011) assumed a wavelength dependence of reddening taken from Cardelli et al. (1989) and used in all cases. As described by Brown et al. (2011), reddening was precomputed for all filters and for a range of stellar parameters. Then, the authors calculated the reddening vectors for typical stars in the Kepler field ( K, , [M/H]) and applied them to all stars. Values of interstellar reddening influenced the other parameters determined from KIC photometry, especially the effective temperatures. For this reason, it is necessary to check the available KIC values by determining them with an independent method.

Figure 3: A comparison of effective temperatures derived by different methods. In the panel (a) the differences (spec)(KIC) are plotted as filled circles, whereas (spec)(H2014) are shown as stars. In the panel (b) (spec)(SED) are plotted. The stars described in the text are indicated by names and open circles. The distributions of differences  and their characteristics are given in the right panels.

Effective temperatures

The effective temperatures determined using the three methods described in Sect. 4 are presented in Table 3 and compared in Fig. 3, where the differences between (spec) obtained from the analysis of the high-resolution spectra and values taken from the KIC, (KIC), and from SED fitting, (SED), are shown as a function of (spec).

For most of the stars, the values of (spec)(KIC) agree to within  K. Differences larger than  K occur for stars with peculiarities identified in their spectra (KIC 3231985, 4150611, and 8703413) and for some rapidly rotating stars with  exceeding  km s (KIC 11506607, 3347643, 4681323, 5471091, and 6450107). The Spearman’s and Kendall’s rank correlation coefficients between (spec)(KIC) and (spec) are equal to and , respectively, indicating very low correlation. The (KIC) values are underestimated for stars hotter than about  K. For stars with lower effective temperatures the agreement is better and the outliers can be explained by rapid rotation or chemical peculiarities in their spectra, which can influence the determination of atmospheric parameters. This result is in agreement with the previous spectroscopic determinations of atmospheric parameters given by Molenda-Żakowicz et al. (2011), Lehmann et al. (2011) and Tkachenko et al. (2013b), where the systematic underestimations of the KIC  values for stars hotter than about  K were pointed out. These deviations probably originate from incorrect values of interstellar reddening.

Recently, Huber et al. (2014) (hereafter H2014) published a new catalogue of revised effective temperatures for Kepler stars. The catalogue is based on a compilation of literature values for atmospheric properties (, and [Fe/H]) derived from different methods, including photometry, spectroscopy, asteroseismology, and exoplanet transits. Effective temperature is available from the H2014 catalogue for all targets in our sample. In Fig. 3a we show the differences between (H2014) and (spec). As one can see, there is a similar scatter in differences [(spec)(H2014)], as was the case for (spec)(KIC). For stars with  K, the differences [(spec)(H2014)] are lower than in the case of (KIC), which means that the (H2014) values are higher than (KIC) for such stars. Neither (KIC) nor (H2014) are systematically underestimated between and  K.

Figure 3 b compares (spec) with that determined from the spectral energy distribution, (SED). The consistency between these two methods is good: there is no significant trend of (spec)(SED) with effective temperature. The Spearman’s and Kendall’s rank correlation coefficients between (spec)(SED) and (spec) are equal to and , respectively. This shows that the correlation is smaller than between (spec)(KIC) and (spec). The differences exceed  K only for three rapidly rotating stars: KIC 5294231, 5786771, and 8915335. The (SED) analysis supports the claimed accuracy of our spectroscopic temperatures.

Figure 4: Left panel: the differences (spec)(KIC) (filled circles) and (spec)(H2014) (stars) as a function of (spec). Right panel: the distributions of differences  and their characteristics.

Surface gravities

Figure 4 compares the (spec) values obtained from high-resolution spectroscopy with those taken from the KIC, (KIC), and the H2014 catalogue, (H2014). For most stars the surface gravities are consistent to within  dex. The distribution of differences (spec)(KIC) and (spec)(H2014) does not indicate any trend with surface gravity or effective temperature. The higher differences were determined mostly for rapidly rotating stars with   km s. Irrespective of rotational velocity, spectroscopic surface gravity values determined from iron lines are adopted for subsequent analysis of these stars. The high-resolution spectra are much more sensitive to  than the photometric indices upon which the KIC is based.

Figure 5: Microturbulent velocities as a function of effective temperature. Non-CP stars are shown as open circles, Am and Ap stars as filled circles and star symbol, respectively. Sizes of symbols for non-peculiar stars depend on rotational velocity, see the legend.

Microturbulent velocities

In Fig. 5 the derived values of microturbulence are shown as a function of effective temperature, with different symbols for CP and non-CP stars. The values of  obtained decrease towards higher effective temperature. For stars with  between and  K,  is within the range  km s, for those with  between and  K,  is approximately  km s, while for even higher values of  microturbulence decreases to about  km s. This behaviour is in agreement with the results of Gebran et al. (2014), based on the analysis of high-resolution spectroscopic data of 61 A field stars, 55 A and 58 F in open clusters (Pleiades, Coma Berenices, Hyades and Ursa Major). The results are also in agreement with Smalley (2004). The determined microturbulent velocities of CP stars are similar to those of non-CP stars. According to Landstreet et al. (2009) the microturbulent velocity in Am stars is generally higher than in normal stars. This tendency is not seen in our sample. However, in Landstreet et al. (2009) it is not possible to directly compare the microturbulence velocities of the Am stars with normal stars of the same temperature, because, as Murphy (2014) notes, the temperature distributions do not overlap. We conclude that in our sample the distribution of microturblent velocities among the Am stars appears to be the same as that of normal stars. As shown by Landstreet et al. (2009), the line profiles of slowly rotating Am stars show different asymmetries and anomalies caused by velocity fields. These asymmetries affect line shapes and equivalent widths and lead to a higher microturbulence velocities in comparison with chemically normal A stars. However, since the analysis of Landstreet et al. (2009) is based on a sample of stars with low projected rotation velocity (  km s), it cannot be ruled out that this phenomenon is diluted for faster rotators. It is beyond the scope of this paper to examine this effect in more detail.

As expected, the microturbulent velocity determined for the Ap star is small,  km s (e.g. Folsom et al., 2007).

Figure 6: Left panel: the differences (spec) (KIC) as a function of (spec). Non-CP and CP stars are plotted by different symbols as explained in the legend. The four most deviant objects are encircled (see the description in the text). The line is the linear fit to all points. Right panel: the distribution of differences and their characteristics.

Iron abundances

Figure 6 compares the iron abundances derived from high-resolution spectra, (spec), with those from the KIC, (KIC). As can be seen, the differences (spec)(KIC) are strongly correlated with the iron abundance. The same situation occurs for differences between (spec) and (H2014), taken from the H2014 catalogue. The Spearman’s and Kendall’s rank correlation coefficients are equal to and , respectively. Only four stars (KIC 4150611, 6951642, 7119530, 7299869; the encircled symbols in Fig. 6) do not follow this general tendency. KIC 4150611 is a possible triple system classified as F1 V with weak metal lines suggesting spectral type A9. This classification was not confirmed by detailed analysis of chemical abundances, as iron abundance is only  dex lower than the solar abundance from Asplund et al. (2009). The other elements do not show significant deviations from the solar values either. This effect can be caused by analysing a triple object as a single one. Additionally, the rotational velocity of this star is  km s. For a star with  K, this causes strong blending of all metal lines. It is worth mentioning that the differences between effective temperatures derived from the KIC photometry, SED, and high-resolution spectrum are also high for this star. Rapid rotation ( km s) is also likely to be responsible for the deviations in chemical abundances for the other three stars. Additionally, the analysis of high-resolution spectra resulted in higher iron abundances for KIC 4768731, an Ap (SrCrEu) star for which a high iron abundance is expected. The same trend of (spec)(KIC) with iron abundance was noticed by Tkachenko et al. (2013b). We found no correlation of (spec)(KIC) differences with effective temperature or surface gravity.

5 Projected rotational velocities

Figure 7: The distribution of rotational velocities for all analysed stars (panel a) and for all non-CP stars (panel b). Panels c, d, e, and f show distributions of  for stars with  K,  K,  K, and  K, respectively.

Rotational velocity is one of the fundamental parameters that influences stellar evolution. It has been frequently demonstrated that stellar rotation strongly depends on spectral type, with most A stars being rapid rotators (e.g. Royer, 2009; Abt, 2009).

Figure 7 a shows the distribution of derived values of . The stars analysed have projected rotational velocities ranging from to about  km s. Stars with  lower than  km s are considered here as slowly rotating stars. There are 21 such objects in our sample. The obvious members of this group are all either Am or Ap (CrSrEu) stars. These objects will be thoroughly discussed in the next section. Another slowly-rotating star, KIC 3231985, also shows small peculiarities its spectrum and was classified as A4 IV Ca weak (A3). Figure 7 b shows the distribution of  for all non-CP stars. One can see that the peak for slow rotation has disappeared and the mean value of  for this sample is equal to  km s. The maximum of the distribution is located around  km s. This result is consistent with the rotational velocities of A stars given in the literature (Royer, 2009). Abt (2000) pointed out that most of the normal A0–F0 main-sequence stars have equatorial rotational velocities greater than  km s. Panels c-f in Fig. 7 show distributions of  for stars with different ranges of . Most of the stars with  K have projected rotational velocities below  km s, consistent with the predictions for the late A and early F stars (Royer, 2009).

There are nine slowly rotating stars in our sample that were not classified as CP objects. Three of them were classified as F stars: KIC 6519869 (F1 IV-Vs), KIC 7661054 (F2.5 V), and KIC 4077032 (F2 IIIs). On average, F stars have smaller rotation velocities than A stars. The rotation velocity drops rapidly from high values ( km s) for late A stars to low values ( km s) for late F stars (Noci, Pomilia, & Ortolani, 1986; Royer, 2009). All slowly rotating non-CP F stars have abundance patterns close to solar, except for some elements. For KIC 4077032 heavy and rare-earth elements are overabundant, while for KIC 6519869 only abundances of elements derived from weak blends are significantly different from solar. Four of the slow rotators were classified as mid-A stars: KIC 10721930 (A5 IV), KIC 3231985 (A4 IV Ca weak (A3)), KIC 5355850 (A5 IVs), and KIC 8386982 (A4 IV/V:). For KIC 10721930 heavy and rare-earth elements are overabundant. For KIC 3231985 almost all elements are overabundant. This star has been identified as a binary (e.g. Dommanget & Nys, 2002; Fabricius et al., 2002). Two objects were classified as early-A stars: KIC 12153021 (A2 IV-V) and KIC 12736056 (A0.5 IIIs).

The occurrence of non-CP A stars with low  has been discussed in many papers (e.g. Takeda et al., 2008; Abt, 2009; Royer et al., 2014). Abt (2009) presented this problem in detail, considereding binarity and the amount of time necessary for slow rotators to become Ap or Am stars by a diffusion mechanism (see also Talon, Richard, & Michaud (2006) and Michaud & Richer (2013)). The frequency of binaries among A stars turned out to be normal. On the other hand, the time-scale for developing peculiarities is rapid for Ap(Si), Ap(HgMn), and Am stars, but slow for Ap(SrCrEu) stars, for which it takes about half of their main-sequence lifetime. According to Abt (2009), if a constant formation rate of field A stars is assumed, about half of the potential Ap(SrCrEu) stars will appear as normal slow rotators and that is why there are normal slow rotators in the A0-A3 group. Therefore, slowly rotating, non-CP stars could be Ap stars that already underwent magnetic braking, but do not yet show chemical peculiarity. Accordingly, the five slowly rotating stars from our sample (KIC 10721930, 12153021, 12736056, 5355850, and 8386982) might become Ap(CrSrEu) stars in the future.

A more probable reason for the low observed  in some A stars is the projection effect. The influence of this effect on the spectral lines of Vega has been investigated by Gulliver et al. (1994). For many years Vega was considered as a standard star with  equal to  km s (e.g. Royer et al., 2007). Gray (1985) first suggested that Vega is a rapid rotator seen pole-on. Gulliver et al. (1994) analysed high-resolution, high signal-to-noise spectra of Vega and found out two distinct types of profiles. The strong lines exhibit classical rotational profiles with enhanced wings, but the weak lines have different, flat-bottomed profiles. Gulliver et al. (1994) used ATLAS9 model atmospheres and SYNTHE synthetic spectra to model the spectral lines and confirmed that Vega is in fact a rapidly rotating, nearly pole-on star with a gradient in temperature and gravity in the photosphere. The equatorial rotational velocity of Vega ( km s) and inclination () were derived by fitting the flat-bottomed line profiles of Fe i4528Å and Ti ii4529Å. To investigate the effect that low inclination has on the observed spectrum of a rapid rotator, high-resolution (R ), and very high signal-to-noise (S/N ) spectra are necessary. Takeda et al. (2008) found that lines of ionized elements in the spectrum of Vega show classical profiles, while the lines of neutral elements display very characteristic shapes: from flat-bottomed to reversed profiles, due to the temperature gradient from the pole to the equator.

Royer et al. (2012) analysed a sample of A stars to find possible rapid rotators seen pole-on. The line profiles in each spectrum were analysed globally, using a Least-Square Deconvolution (LSD) method (Donati et al., 1997; Kochukhov et al., 2010) to enhance the signal and recover the broadening function. To investigate the possibility that the slow rotators in our sample are in fact rapid rotators seen pole-on, we followed the methods presented by Gulliver et al. (1994) and Takeda et al. (2008). We used the revised LSD code (Tkachenko et al., 2013a) to calculate the LSD profiles for intrinsic line depths from to and from to , which enables the distinction between the line profiles of strong and faint lines in the spectra of the slow rotators. All line lists were computed with the individual parameters of , , , and [M/H] determined in this paper. The boundary for the separation of line depths was chosen so that the S/N in both LSD profiles, computed from weak and strong lines, were approximately equal. The wavelength range used (4890 to 5670 Å) included the region between H and the first strong telluric lines. We found no significant difference in profile shapes for any of the investigated stars, which can be attributed to the insufficient quality of our spectra.

6 Chemical abundances

Figure 8: The difference between the obtained chemical abundances and solar values (Asplund et al., 2009) as a function of the atomic number. The results for non-CP stars are given as filled circles; for CP stars as filled stars.
Figure 9: The distributions of the abundances of the elements: Fe, Sc, Ti, Cr, Mn, and Ca in the analysed sample of stars. The solar abundances taken from Asplund et al. (2009) are shown as vertical lines.
Figure 10: Abundances of different elements as a function of rotational velocity for CP stars (triangles), non-CP stars with  K (filled circles), and non-CP stars with  K (open circles).
Figure 11: Abundances of different elements as a function of iron abundance for CP stars (triangles), non-CP stars with  K (filled circles), and non-CP stars with  K (open circles).
El. No. of Normal All Solar
all stars stars values
stars
C
N
O
Na
Mg
Al
Si
P
S
K
Ca
Sc
Ti
V
Cr
Mn
Fe
Co
Ni
Cu
Zn
Ga
Sr
Y
Zr
Ba
La
Ce
Pr
Nd
Sm
Eu
Gd
Dy
Er
Table 5: Average abundances for the entire sample and the non-CP stars compared to the solar values of Asplund et al. (2009).

In Table 5 the average abundances for the stars in our sample are compared with the photospheric solar values of Asplund et al. (2009). Only the abundances of Mg, Ti, and Fe were derived for all stars. We excluded the O i 7771-5 Å triplet from our analysis and used other O i lines (e.g., 5331, 6157, 6158, and 6456 Å) that are much less affected by non-LTE effects (e.g. Przybilla et al., 2000). Lines of heavy elements except for Sr, Y, Ba, and rare-earth elements are very sparse in spectra of A stars. Abundances of these elements were investigated only for slowly and moderately rotating stars and in most cases, only one or two blends were available.

In Fig. 8, the determined average abundances for all stars are compared with the solar abundances of Asplund et al. (2009). In Fig. 9, histograms show the distribution of abundances of the iron-peak elements, Ti, Cr, Mn, Fe, Sc, and Ca. The average abundances of most of the light and iron-peak elements are close to the solar values. The largest differences occur for the elements represented by weak blends only. Naturally, the abundances of some elements determined for CP stars are significantly different from the solar values. For instance, the highest Cr abundance, , was determined for the Ap CrSrEu star KIC4768731, and the lowest abundances of Sc were derived for Am stars.

We found that apart from 14 CP stars, there are two low-metallicity objects in our sample, KIC 9828226 and KIC 8351193. KIC 9828226 was classified as an A1.5 V star for which the metal line spectrum is weak and matches A0 type. The moderate rotational velocity,  km s, allowed us to obtain the abundances of Mg, Si, Ca, Sc, Ti, Cr, Fe, Ni, Sr, C, and N. All these elements are underabundant. KIC 8351193 (HD 177152) was classified as a B9.5: star. With a rotational velocity of  km s the spectrum of this star is almost featureless. It was possible to determine only the abundances of Mg, Ti, and Fe on the basis of single lines. All these elements are underabundant. However, because of the high , a small error in continuum location significantly changes the results. Tkachenko et al. (2013b) analysed the high-resolution TLS spectrum of this star and obtained  K and  dex from the analysis of Balmer lines. The effective temperature and surface gravity obtained in the present paper, (spec)  K and (spec)  are higher than given by Tkachenko et al. (2013b). The high rotational velocities obtained in both papers are consistent within the error bars. Although Tkachenko et al. (2013b) also found a low metallicity for this star, their iron abundance of  dex is much lower than the value obtained in this paper ( dex). The difference is caused mostly by differences in effective temperatures. The effective temperatures (SED) and (spec) of KIC 8351193 obtained in this paper are consistent with each other and with the values  K and  K, determined from the SED and Strömgren photometry by Balona et al. (2011). On the basis of Kepler data, Balona et al. (2011) classified KIC 8351193 as a “rotationally modulated” star with period of  d. We confirm this classification (see Sect. 7).

6.1 Correlation of chemical abundances with other parameters

The analysis of a large sample of stars allowed us to search for correlations between the abundances of chemical elements and stellar parameters. We analysed separately three groups: the CP stars formed one group, and the normal stars were separated into hotter ( K) and cooler ( K) populations. For some elements, large scatter or a small number of stars with determined abundances prevents the correlation analysis. This is the case for most of the heavy and rare-earth elements. Comparing samples with hotter and cooler stars, we found a larger scatter of abundances for hotter stars, especially for  higher than  km s. The same effect was observed for the abundances of A and F stars in the Hyades (Gebran et al., 2010), Pleiades (Gebran & Monier, 2008) and Coma Berenices open clusters (Gebran, Monier, & Richard, 2008).

No significant correlations of element abundances with effective temperature and surface gravity were found in the normal stars, with absolute values of Spearman’s and Kendall’s correlation coefficients less than . No significant correlations between the abundances of most elements and rotational velocity were found either. Typical relations for some elements are presented in Fig. 10. A small positive correlation was found between the abundances of C and  for hotter stars, with correlation coefficients of about 0.3–0.4 according to both considered statistics. In an investigation of 23 normal A and F stars, Fossati et al. (2008) found no correlation between abundances and rotational velocity. A larger sample was examined by Takeda et al. (2008), who found negative correlations for C, O, Ca and strong positive correlations for the Fe-peak elements plus Y and Ba, with rotational velocity.

Significant positive correlations with iron abundance in normal stars were found for almost all elements (see Fig. 11). The highest correlation coefficients were obtained for light elements Na, Mg and Si, iron-peak elements Ca, Ti, Cr, Mn, Ni, and heavy elements Sr and Ba, all of which have correlation coefficients larger than . Less significant positive correlations were obtained for Sr and Y, for which both correlation coefficients were smaller than . Takeda et al. (2008) also found strong positive correlations between the abundances of Si, Ti and Ba, and that of Fe.

We classified 21 CP stars in our sample, including 17 Am stars, one Ap CrSrEu, two  Boo stars and one He-strong B star. We checked for correlations between abundances in Am stars and their atmospheric parameters, iron abundances, and rotational velocities. The abundances of most chemical elements appeared to be independent of the effective temperature. The negative correlation determined for O and the positive correlations found for Cu, Zn, and Zr could be caused by the small sample size. We found no correlations between abundances and surface gravity. On the other hand, the abundances of Na, Si, Ti, V, Cr, Mn, Ni, Cu, and rare-earth element La were all positively correlated with the iron abundance (Fig. 11). The only element found to be negatively correlated with iron was Zr, with a correlation coefficient of , according to Spearman’s statistic. The abundances of many elements were also found to be strongly correlated with that of iron for non-CP stars.

The correlation between Na and Fe abundances was discussed by Takeda et al. (2012), who performed a non-LTE analysis of the abundances of alkali elements (Li, Na, and K) determined from the Li i  Å, Na i / Å, and K i  Å lines for 24 slowly rotating A stars ( km s) in the broad temperature range  K  K. Many of these stars showed Am peculiarities to different degrees. They found a significant trend of Na abundances with Fe, independent of . This means that Na becomes enriched along with Fe in accordance with the degree of Am peculiarity. The same trend was derived for our sample of stars with a Spearman’s correlation coefficient of . The positive correlation of Na abundance with Fe disagrees with the predictions of atomic diffusion theory (e.g. Talon, Richard, & Michaud, 2006). The same is true of Si abundances in Am stars, where we see a positive correlation of , according to Spearman’s statistic, that is not theoretically expected.

No significant correlations with  were found for the abundances of any element, albeit in a small analysed range of . This contradicts the finding of Fossati et al. (2008) who analysed Am stars in the Praesepe cluster and reported strong correlations between abundances and  of Am stars for all peculiar elements, except for Sc and Ti. One reason for the contradiction may be that the field stars analysed here have different ages and initial chemical composition.

6.2 Chemically peculiar stars

Figure 12: Chemical abundances of Am stars compared with the solar values (Asplund et al., 2009) as a function of the atomic number.
KIC number HD number SpC (new)
New Am stars
3441230 kA2hA5mF0 V Am
5199439 225365 kA3hA5mA5 (IV)s Am:
5200084 225410 kA3hA6mF1 (IV) Am
5633448 225463 kA3hA4mA7 V Am
6292398 kA3hA7mA7 V Am:
6292925 kA2.5hA3mA7 (IV) Am:
7767565 196995 kA5hA7mF1 IV Am
Known Am stars
8323104 188911 kA2.5hA6mA7 (IV) Am
8703413 187254 kA3hA5mA9 (IV)s Am
9117875 190165 kA3hF0.5mF3 (III) Am
9147002 180237 kA3hA5mF3 V Am
9204718 176843 kA3hA9mF1 (IV) Am
Am stars in SB2 systems
8112039 187091 kA3hA5mA5 (IV)s Am:
4660665 kA3hA5mA5 IV Am:
6865077 kA3hA8mA6 (IV) Am
8692626 kA2hA4mA6 (IV) Am:
10616138 kA3hA6mA7 (IV)s Am:
Ap star
4768731 225914 A5 Vp (SrCrEu)
 Boo stars
9828226 hA2kA0mB9 Vb.
5724440 187234 F0 Vnn weak met (A5)
Table 6: CP stars investigated in this work.

There are two main classes of CP stars among the main-sequence Population I stars: the metallic line (Am) and the peculiar A stars (Ap). The Am stars manifest overabundances of most iron-peak elements and some heavy elements like Zn, Sr, Zr and Ba, and underabundances of Ca and Sc (Gray & Corbally, 2009). These anomalies are explained by radiative diffusion in a non-magnetic star (e.g. Richard, 2005; Talon, Richard, & Michaud, 2006). Their projected rotational velocities are generally smaller than those of normal A stars (Abt & Hudson, 1971), which permits the segregation of elements by diffusion (Abt & Moyd, 1973). The majority of Am stars are members of close binary systems, where their slow rotation results from tidal braking.

The Ap stars form the second group of CP stars with overabundances of some elements, such as Si, Sr, Cr, and rare-earth elements (e.g. Kurtz et al., 2007a, b). Most of these stars have strong magnetic fields with polar field strengths typically of several kG (Hubrig et al., 2005; Kurtz et al., 2006). These stars have much slower rotational velocities than normal A and B stars. The distribution of chemical abundances changes both horizontally, with the formation of spots, and vertically, where they are stratified (e.g. Ryabchikova, Leone, & Kochukhov, 2005; Ryabchikova, 2014; Kochukhov et al., 2015). The abundance anomalies and spots are caused by diffusion in the presence of a magnetic field (e.g. Babel, 1993; Turcotte, 2003).

There are other classes of CP stars, such as the  Bootis stars, which are a class of metal-weak population I A stars (Gray & Corbally, 2009). They are characterised by a weak Mg II4481 Å line and solar abundances of C, N, O, and S (Gray, 2005). Their projected rotational velocity distribution matches that of normal A stars (Abt & Morrell, 1995; Gray & Corbally, 2009). The metal deficiencies of  Boo can be distinguished from those of Population II stars on the basis of detailed abundance analyses (Paunzen et al., 2014). This suggests that the abundances of  Boo stars are restricted to the outer layers of their atmospheres and can be explained by accretion of dust-depleted circumstellar material (Venn & Lambert, 1990; King, 1994).

There are 20 peculiar A stars in our sample (see Table 6), which are comprised of 17 Am stars, one Ap star (KIC 4768731), and two  Bootis candidates (KIC 9828226 and KIC 5724440).

Four of the Am stars were previously known from multicolour photometry (Mendoza, 1974) and spectroscopy (Floquet, 1975). These are KIC 8323104, 8703413, 9117875, and 9204718. KIC 8323104 has a very low   of  km s, which enabled us to derive the abundances of alkali elements Na and K. Determinations of Na and K abundances in the atmospheres of Am stars are rare (Takeda et al., 2012), but important for establishing the theory explaining the chemical abundance peculiarities of A stars. The effective temperature of KIC 9117875, at  K, is low in comparison with other Am stars in our sample. The star falls among the  Dor stars in the  – diagram, but it is not a  Dor star (see Sect. 7). The fifth known Am star, KIC 9147002 was previously classified as an A2p star by Bertaud (1960). Its abundance pattern is typical of Am stars, with a significant Sc underabundance and strong overabundances of heavy elements. The Am stars KIC 8703413, 9117875 and 9204718, were also studied by Catanzaro & Ripepi (2014). Their results are in agreement with the values obtained here, within the errors.

Six newly discovered Am stars are KIC 3441230, 5199439, 5200084, 5633448, 6292398, and 6292925. These have low and abundance patterns typical of Am stars. A seventh Am star is KIC 7767565, of spectral type kA5hA7mF1 IV, whose chemical abundances are not typical, because Ca and Sc are overabundant. This star has a moderate  of  km s.

KIC 8112039 is a highly eccentric binary, also known as KOI-54 (Welsh et al., 2011; Burkart et al., 2012). Welsh et al. (2011) analysed the high-resolution spectrum taken in quadrature and determined atmospheric parameters of both components, finding and equal to and  K, and surface gravities and , respectively. Both stars appeared to be metal-rich, with [Fe/H]  dex. They also determined low  values for the components, at and  km s. The spectrum investigated in this paper was taken in conjunction, so that the lines of both stars overlap. From the analysis of our spectrum we derived  K,  dex, [Fe/H]  dex, and  km s.

Four other stars were classified as Am stars and belong to SB2 binary systems: KIC 4660665, KIC 6865077, KIC 8692626, and KIC 10616138. These stars, together with KIC 8112039, will be analysed in a forthcoming paper (Catanzaro et al., in preparation).

The only Ap star in our sample, KIC 4768731, has the spectral type A5 Vp (SrCrEu). According to Mason et al. (2001), the star is the primary of a visual binary with a separation of about  arcsec with  mag (primary) and  mag (secondary). Previously, KIC 4768731 was classified as an A7 star (see Table 2), without indication of chemical peculiarity. The chemical pattern of KIC 4768731 is typical for Ap (SrCrEu) stars, with significant overabundances of Cr and Sr (see Fig. 13). All Fe-peak elements, except for Sc, V, and Zn, are overabundant. Additionally, KIC 4768731 was discovered as a rapidly oscillating (roAp) star (Smalley et al., in preparation). Detailed spectral analysis of Ap stars requires taking into account stratification of the elements in their atmospheres. This effect has a substantial impact on spectral lines, so atmospheric parameters determined here from the ionisation and excitation equilibria and wings of hydrogen lines can be incorrect. Results of spectral analysis with stratification of chemical elements for KIC 4768731 will be presented by Niemczura, Shulyak et al. (in preparation).

We classified two stars as probable  Boo stars. The spectral type obtained for KIC 9828226 is hA2kA0mB9 Vb. Its hydrogen lines are very broad, the metal lines are weak, and the magnesium Mg ii4481 line is especially weak when compared with the metal line type. The star has  km s. KIC 5724440 was given the spectral type F0 Vnn (weak met (A5)). Its hydrogen lines are broad and shallow – typical for a rapid rotator – and are consistent with an F0 spectral type, whereas its metal lines are weak and are similar in strength to the A5 standard HD 23886. The Ca ii K line is broad but shallow, having the depth of the A3 V standard but the width of the A7 V standard. The moderate resolution spectrum of this star was analysed by Catanzaro et al. (2011), who indicated a problem with the effective temperature determination from Balmer lines. This problem was also noticed in this paper. The Mg ii4481 in the spectrum of KIC 5724440 is weaker than in the spectrum of a normal star with the same atmospheric parameters. The iron-peak elements and especially Sr are slightly underabundant. This supports the  Bootis classification of this star.

Figure 13: Chemical abundances of the Ap star KIC 4768731 compared with the solar values (Asplund et al., 2009) as a function of the atomic number.
Figure 14: Obtained chemical abundances of the suspected  Boo stars KIC 5724440 (filled circles) and KIC 9828226 (filled squares) compared with the solar values (Asplund et al., 2009) as a function of the atomic number.

7 Discussion and Conclusions

We analysed high-resolution HERMES/Mercator spectra of 117 A and F stars collected during the spectroscopic campaign of Kepler stars. Spectral classification showed that our sample consists mostly of non-evolved stars of luminosity types V to IV and revealed 20 chemically peculiar (CP) stars of Am, Ap (CrSrEu), and  Bootis types. These peculiarities were confirmed by detailed investigation of the stellar spectra. Nine new CP stars were found, of which six were Am, one was Ap, and two were  Boo stars.

The atmospheric parameters, i.e. , , , abundances of chemical elements, and  for all stars were determined by using spectral synthesis. The atmospheric parameters were compared with the values taken from the KIC and those determined from SED fitting. Similar to the results reported by Tkachenko et al. (2013b) and Lehmann et al. (2011), the KIC values of effective temperatures were found to be lower than the spectroscopic ones for stars with   K. One of the reasons for these differences is probably the use of incorrect interstellar reddening when estimating the KIC temperatures. The E(BV) values from the KIC were compared with the spectroscopic values determined from interstellar Na D lines. For  higher than  K the differences E(BV)(Na)–E(BV)(KIC) increase. The spectroscopic effective temperatures derived from iron and Balmer lines are in good agreement with values obtained from SED fitting. For most stars, the surface gravities determined from iron lines and taken from the KIC are consistent to  dex. There is a positive correlation between the differences of iron abundances determined from spectroscopy and taken from KIC and . The same trend was noticed by Tkachenko et al. (2013b). This suggests that the KIC iron abundances of metal-weak stars are systematically overestimated, whereas iron abundances of metal-strong stars are underestimated. The obtained microturbulence velocities follow the trend presented by Landstreet et al. (2009) and are smaller towards higher effective temperatures. The obtained values of  range from to about  km s with the average equal to  km s, which is typical for A and early F stars.

Figure 15:  –  diagram for the stars analysed in this paper. Stars identified as single  Sct are indicated by filled circle,  Sct in a binary systems as open circle, single  Dor as filled squares,  Dor in binary systems as light blue squares. Single stars with hybrid pulsations are plotted as filled diamonds, whereas hybrid stars in binary systems are open diamonds. Non-pulsating stars (including Am objects) are shown as stars.

The atmospheric parameters determined in this work are very important as input ingredients for further seismic modelling. We analysed A and F stars, so our sample contains pulsating  Sct,  Dor, and hybrid stars. To find out if the obtained atmospheric parameters of the stars are consistent with theoretical instability strips in this part of the H-R diagram, we performed an analysis of the Kepler data for all stars investigated here. We used the corrected light curves available through the MAST 15 database. We considered all the long cadence data available for a given star. If necessary, short cadence data were used to verify our classification. The Kepler data from each quarter were shifted to have the same mean level and all quarters were merged. For some stars, it was necessary to remove the trends from each quarter separately. Then, we visually examined all light curves, the frequency spectra, and detected frequencies to identify candidate  Sct,  Dor, hybrid, and other types of variability. The results of our initial classification are presented in the last column of Table 2. We classified 15 stars as  Sct, 13 as  Dor, and 71 as  Sct/ Dor hybrid stars. In this paper, all objects for which we can find frequencies typical for  Sct and  Dor stars are termed “hybrid”. Moreover, we found 14 pulsating stars in binary systems, including eclipsing binaries, and 30 stars which do not show pulsations (these stars were discussed in detail by Murphy et al. (2015). The spectroscopic parameters provided here were critical for assigning locations of the non-pulsators with respect to the  Sct instability strip. Here we only mention that light curves of the stars belonging to the latter group show variability resulting from rotation and spots on the surfaces of these stars and/or from binarity (see Table 2).

All of the stars investigated are shown in the    diagram (Fig. 15). We add evolutionary tracks, calculated with Time Dependent Convection (TDC, Grigahcène et al., 2005), for stars having [Fe/H]  and covering a range of masses from to  M. The stars are plotted with different symbols depending on the type of pulsations, type of chemical peculiarity and other characteristics. As we can see in Fig. 15, non-pulsating stars are situated inside and on the blue edge of the  Sct instability strip. Most Am stars are located inside the instability strip. Most hybrid stars are located inside or close to the blue edge of the  Sct instability strip. The exceptions are KIC 12736056 and 8489712, whose effective temperatures are too high for hybrid  Sct/ Dor stars. The same applies for KIC 4056136, classified as a  Sct star, KIC 11180361, showing variability typical for hybrid stars and binarity, and for KIC 6292398, for which we found  Dor frequencies and evidence for binarity. All these stars need further investigation to verify the possibility of contamination by another star. The frequency analysis of the Kepler data revealed 31 stars with low-frequency peaks which we attribute to binarity. For these stars lines from only one star are present in the spectrum and we see  Sct and/or  Dor frequencies in the Kepler data.

We investigated 13 Am type stars, of which 7 are newly discovered. Most of them show only the variability typical for stars in binary systems or due to rotation and spots. Only three stars are possible pulsators. KIC 6292398 has a frequency pattern found in  Dor stars and binarity, KIC 7767565 is probably also a  Dor star, and KIC 9204718 shows frequencies typical of binarity, rotational changes and  Sct pulsations. The Kepler photometric data for these stars need to be analysed in detail.

We have also investigated a possible relationship between metallicities or chemical peculiarities and the hybrid pulsations. The first  Sct/ Dor hybrid, HD 8801, was discovered by Henry & Fekel (2005) and classified as an Am star. This classification was confirmed by Neuteufel et al. (2013). Hareter et al. (2011) analysed the spectra of two other hybrid stars, HD 114839 and BD +18 4914, discovered with MOST photometry by Rowe et al. (2006) and King et al. (2006), respectively. The evolved star BD +18 4914 has an abundance pattern typical for an Am star, but the less evolved star HD 114839, previously classified in the literature as a mild Am star, does not show any peculiarity. For 58 hybrid candidates investigated here the spectral analyses do not indicate any peculiarity.

Acknowledgments

We are grateful to the anonymous referee for their thorough examination of this paper. EN acknowledges the support from Fundaç ao para a Ciência e a Tecnologia (FCT) through the grant ‘Coopera¸c ao Cientica e Teccnologica FCT/Polonia 2011/2012 (Proc. 441.00 Polonia)’ funded by FCT/MCTES, Portugal and the support from 1007/S/IAs/14 funds. This work was supported by NCN grant 2011/01/B/ST9/05448. DD would like to acknowledge the NSC grant No. 2011/01/N/ ST9/00400. Calculations have been carried out in Wrocław Centre for Networking and Supercomputing (http://www.wcss.pl), grant No. 214. This research was supported by the Australian Research Council. Funding for the Stellar Astrophysics Centre is provided by the Danish National Research Foundation (grant agreement no.: DNRF106). The research is supported by the ASTERISK project (ASTERoseismic Investigations with SONG and Kepler) funded by the European Research Council (grant agreement no.: 267864). KU acknowledges support by the Spanish National Plan of R&D for 2010, project AYA2010-17803. The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007–2013)/ERC grant agreement No. 227224 (PROSPERITY). This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. MB is F.R.S.-FNRS Post-doctoral Researcher, Belgium. AP acknowledges support from the NCN grant No. 2011/03/B/ST9/02667. PIP is a Postdoctoral Fellow of the Fund for Scientific Research (FWO), Flanders, Belgium. SB is supported by the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO).

KIC HERMES N V SpT1 SpT2 Notes
Number observations [mag] [this paper]
1294756 2010 May-September 1 8.96 A2 A3 IV mult. + rot. +  Sct
2571868 2011 July 1 8.67 A0 A3 IVn hybrid (low. ampl.  Dor)
2694337 2010 May-September 1 10.35 F0.5 IVs hybrid
2695344 2012 July 1 9.62 F2 F2 V SB2;  Dor
2012 September 1
3230227 2010 May-September 1 8.94 A5 A5 IV: SB2; ecl. bin. +  Sct
3231985 2011 September 1 9.22 A2 A4 IV Ca weak (A3) rot.
3347643 2010 May-September 1 7.69 A2 A3 Van hybrid (low. ampl.  Dor)
3441230 2012 July 1 9.79 A2 kA2hA5mF0 V Am rot.
3656913 2012 July 1 9.81 A0 kA2hA6mF0 IV Am SB2; bin.
3850810 2010 May-September 1 10.10 F1 Vs wk met (F0) hybrid + rot.
3851151 2010 May-September 1 9.82 A2 A3 V hybrid + rot.
3942392 2012 July 1 9.87 F2 F2 Vs rot. +  Dor
4035667 2010 May-September 1 10.05 A3 IVn hybrid + bin. + rot.
4044353 2010 May-September 1 9.76 A2 A5 IVs hybrid
4048494 2011 July 1 9.55 A6 V SB2; bin. +  Sct
4056136 2011 September 1 9.55 A0 B9 IV-Vnn rot.;  Sct; contamination
4077032 2011 July 1 9.71 F0 F2 IIIs hybrid (low ampl.  Dor )
4150611 2010 May-September 1 8.00 A2 F1 V wk met (A9)  Sct + ecl. bin.
4281581 2011 July 1 9.32 A2 A3 IV-Vs hybrid
4480321 2011 September 1 10.27 A9 V wk met (A5) SB2; hybrid + bin. + rot.
4572373 2012 July 1 9.84 A3 Van (wk met A2) rot.
2012 September 1
4660665 2011 September 1 7.86 A2 kA3hA5mA5 IV Am: SB2; bin. + rot. +  Dor
4671225 2011 September 1 9.88 A7 A3 IV-Vn hybrid
4681323 2011 July 2 9.06 A0 A1 IV-s bin. + rot.
4768731 2011 July 1 9.17 A7 A5 Vp SrCrEu Ap; bin. + rot.
4831769 2011 September 1 9.55 A2 A3 Va bin.
4832225 2011 September 1 9.08 A0 B9.5 V SB2; bin. +  Dor
4840675 2012 July 1 9.66 A5 Vnn SB2; hybrid (low ampl.  Dor )
5018590 2012 July 1 9.87 F0.5 V  Dor
2012 September 1
5113797 2011 July 1 9.17 A3 A3 IV-V hybrid (low ampl.  Sct)
5199439 2012 July 1 9.51 A0 kA3hA5mA5 (IV)s Am: bin. + rot.
2012 September 1
5200084 2011 September 1 9.16 A0 kA3hA6mF1 (IV) Am bin.
5294231 2012 July 1 9.74 A2 A0 IVn rot.
2012 September 1
5355850 2011 September 1 9.44 A5 A5 IVs bin. + hybrid
5429163 2011 July 1 9.80 A5 V A4 V hybrid? + rot.
5471091 2012 July 1 9.97 F2 F2 IV bin. + rot.
2012 September 1
5473171 2010 May-September 1 8.97 A2 A5 IV/Vn hybrid
5524045 2011 September 1 9.36 A0 A0.5 Va+ rot.
5525210 2012 July 1 9.90 F1 V SB2; bin. + rot.
2012 September 1
5608334 2011 July 1 9.88 F2 Vs wk met (A7) hybrid?
5633448 2011 September 1 8.96 A2 kA3hA4mA7 V Am rot.
5724440 2011 July 1 7.89 A5 F0 Vnn wk met (A5)  Boo?; hybrid
5786771 2011 July 2 9.05 A2 A0.5nn: bin. + rot.
5954264 2010 May-September 1 8.16 F0 F1 IV-V  Dor
5988140 2010 May-September 1 8.81 A7.5 IV-III A8 IIIs  Sct + bin.
Table 2: Journal of spectroscopic observations and the derived spectral classes. is the number of available spectra; SpT1 indicates the spectral classification from the literature, whereas SpT2 means spectral classifications obtained in this work. Superscript in the first column denotes the stars with additional information below the table. Additional notations for luminosity class: “b” – lower luminosity main-sequence star; “a” or “a+” – higher luminosity main-sequence star; “IV-V” – between IV and V; “IV/V” – either IV or V; “n” – nebulous (i.e. broad-lined); “s” – sharp lines; “nn” – very rapid rotators; “wk” met – weak metal lines; “met str” – strong metal lines. In the last column additional information taken from spectra and Kepler data are presented. “SB2” means double-lined spectroscopic binary; “rot.” – rotational frequency detected; “bin.” – orbital frequency detected; “mult.” – orbital frequencies of both components detected; “hybrid” – frequencies of  Dor and  Sct visible in the spectra; “contamination” – Kepler data contaminated by a nearby star; “?” – our classification of variability type is uncertain.
KIC HERMES N V SpT1 SpT2 Notes
Number observations [mag] [this paper]
6032730 2011 July 1 8.61 A2 A3 IVn hybrid
6106152 2012 July 1 8.06 A0 A3 IV wk met (A1) bin.
6128236 2011 September 1 8.79 F2 F0 IV+ bin. + rot.
6192566 2011 September 1 9.31 A2 A6 IV-V SB2; bin. + rot. ?
6292398 2011 September 1 9.98 kA3hA7mA7 V Am: bin. +  Dor
6292925 2011 September 1 9.66 kA2.5hA3mA7 (IV) Am: bin.
2012 July 1
6352430 2011 September 1 7.91 B8 B8 Vs He str (B6) SB2;  Dor
6380579 2011 September 1 9.99 F2 Vs  Dor
6381306 2011 September 1 8.63 A0 SB2; hybrid
6450107 2011 September 1 7.53 A0 A1 IV-s rot.
6509175 2011 September 1 9.95 A2 A6 IVn: hybrid
6519869 2011 September 1 10.38 F1 IV-Vs  Dor
6587551 2010 May-September 1 9.72 A0 A3 IV hybrid
6756386 2010 May-September 1 8.62 A2 A3 IVn hybrid
6761539 2011 September 1 10.27 F0 V wk met (A5) hybrid
6844024 2012 September 1 9.92 F1 Vn: wk met (A6) bin. + hybrid
6865077 2011 July 1 9.76 kA3hA8mA6 (IV) SB2; hybrid
6937758 2010 May-September 1 9.72 A2 A5 V hybrid (low ampl.  Dor)
6951642 2011 September 1 9.70 A5 F1 V hybrid; contamination?
6967360 2011 September 1 9.42 A2 A1 Vann hybrid
7060333 2011 September 1 9.09 A2 A5 IV bin. +  Sct
7119530 2010 May-September 3 8.44 A3 A3 IVn hybrid
7299869 2012 July 1 9.70 A2 A3 IV-n hybrid
2012 September 1
7345479 2011 September 1 7.90 A0 A2 Vnn bin. + rot.
7530366 2011 September 1 8.31 A0 A0.5 IVnn bin. + rot.
7533694 2011 September 1 10.29 A2.5 IVn hybrid
7583939 2011 July 1 9.64 A0 A2 IV/V hybrid (low amp.  Dor)
7661054 2012 July