arxiv:astro-ph/ v1 13 Jan 2004

Transcription

1 Astronomy & Astrophysics manuscript no. mg3 November 13, 2018 (DOI: will be inserted by hand later) On the correlation of elemental abundances with kinematics among galactic disk stars. T.V. Mishenina 1, C. Soubiran 2, V.V. Kovtyukh 1, and S.A. Korotin 1 arxiv:astro-ph/ v1 13 Jan Astronomical Observatory of Odessa National University and Isaac Newton Institute of Chile, Odessa Branch, Ukraine 2 Observatoire Aquitain des Sciences de l Univers, CNRS UMR 5804, BP 89, Floirac, France Received accepted Abstract. We have performed the detailed analysis of 174 high-resolution spectra of FGK dwarfs obtained with the ELODIE echelle spectrograph at the Observatoire de Haute-Provence. Abundances of Fe, Si and Ni have been determined from equivalent widths under LTE approximation, whereas abundances of Mg have been determined under NLTE approximation using equivalent widths of 4 lines and profiles of 5 lines. Spatial velocities with an accuracy better than 1 kms 1, as well as orbits, have been computed for all stars. They have been used to define 2 subsamples kinematically representative of the thin disk and the thick disk in order to highlight their respective properties. A transition occurs at [Fe/H] = 0.3. Stars more metal-rich than this value have a flat distribution with Z max < 1 kpc and σ W < 20 kms 1, and a narrow distribution of [α/fe]. There exist stars in this metallicity regime which cannot belong to the thin disk because of their excentric orbits, neither to the thick disk because of their low scale height. Several thin disk stars are identified down to [Fe/H] = Their Mg enrichment is lower than thick disk stars with the same metallicity. We confirm from a larger sample the results of Feltzing et al (2003) and Bensby et al (2003) showing a decrease of [α/fe] with [Fe/H] in the thick disk interpreted as the signature of the SNIa which have progressively enriched the ISM with iron. However our data suggest that the star formation in the thick disk stopped when the enrichment was [Fe/H] = 0.30, [Mg/Fe] = +0.20, [Si/Fe] = A vertical gradient in [α/fe] may exist in the thick disk but should be confirmed with a larger sample. Finally we have identified 2 new candidates of the HR1614 moving group. Key words. Stars: fundamental parameters Stars: abundances Stars: kinematics Stars: atmospheres Galaxy: stellar content 1. Introduction The behaviour of stellar elemental abundances with age and kinematics, in various substructures of the Galaxy, puts important constraints to the construction of models of chemical and dynamical evolution of our galactic system. In this paper we put particular attention to the transition between the thin disk and the thick disk and to the abundances of the α-elements Mg and Si and the iron-peak element Ni. According to the current nucleosynthetic theory, α-elements are being produced as a result of α-capture reaction, taking place in the core of massive starsduringtheirexplosionassnii(burbidgeetal1957). Fe is produced by both massive SN II and less massive SN Ia. It is expected that 2/3 of the solarabundance of iron is produced by explosions of white dwarfs in double systems Send offprint requests to: tamar@deneb.odessa.ua Based on spectra collected with the ELODIE spectrograph at the 1.93-m telescope of the Observatoire de Haute Provence (France) (SN Ia), 1/3 by SN II (Timmes et al 1995). If the percentage of massive stars in the earlier Galaxy was higher than today, one can predict that the α/fe ratio is going to change over time. A well established observational fact shows that in old metal-poor stars of the Galaxy, [α/fe], in particular [Mg/Fe], is overabundant relative to Sun s value (Wallerstein 1961; Gratton & Sneden 1987; Magain 1989; Nissen et al 1994; Fuhrmann et al 1995 etc). If the α-element overabundance is a typical chemical pattern in halo stars in comparison with disk stars, there is a question whether there is a distinction in the Mg behaviour in other subsystems of the Galaxy, in particular in the thin and thick disks. Such a difference would have important consequences on the choice of the most probable scenario of formation of the thick disk (collapse, accretion etc.) and its timescale. Moreover [α/fe] could serve as an additional criterion for the deconvolution of thick and thin disk populations. It is especially important for the stars with metallicity in the interval of overlapping [Fe/H] values ( 0.7 <[Fe/H] < 0.3).

2 2 Mishenina et al.: [Mg/Fe] vs kinematics in the disk Several groups have recently attacked this problem by attempting to make consistent analyses of abundances of Mg and other elements in large (> 50) samples of stars of the thin and thick disks, using kinematics to distinguish the stellar populations. However, their results are ambiguous. Fuhrmann (1998) found the thick disk to be significantly older than the thin disk and its data exhibit an appreciable constant overabundance of MgI in the thick disk relative to the thin disk. He interprets the observed discontinuity in the chemical and kinematical behaviour of the thin disk and thick disk populations as a hiatus in star formation before the earliest stages of the thin disk. Mashonkina s group (Mashonkina & Gehren 2000, Mashonkina & Gehren 2001, Mashonkina et al 2003) is in agreement with Fuhrmann concerning the discontinuity between the thin disk and the thick disk. They observe a step-like decrease of Eu/Ba ratio at the thick disk to thin disk transition and estimate that the thick disk formed from well mixed gas during a short timescale of 0.5 Gyr. Gratton et al (2000) argue in favour of a constant and significant overabundance of Mg with respect to Fe in the thick disk, a fast formation of the thick disk, a sudden decrease in star formation during the transition between the thick and thin disk phases and a formation scenario mixing dissipational collapse and accretion on similar time scales. Contrary to these studies, Chen et al (2000) do not observe in their data any discontinuity between the thin disk and the thick disk. Their thick disk stars exhibit a slightly higher α/fe ratio than the thin disk, but not exceeding Very recently Bensby et al (2003) and Feltzing et al (2003) published a convincing study showing that the thick disk extents at solar metallicities with an abundance trend clearly separated from the thin disk despite a large overlap in metallicity. A knee and decrease in [Mg/Fe] vs [Fe/H] is interpreted asthe signatureofcontribution from SNIa to the enrichment of the interstellar gas out of which the thick disk stars formed. However these studies cannot be easily compared because they do not rely on the same definition of the thick disk.the thin disk andthe thick diskareknowntogreatly overlap in kinematics and metallicity, and the selection of a sample representative of each population is not obvious. The use of a single criterion (V rot, Z max or eccentricity of the orbit, [Fe/H],...) is usually not sufficient for this task. In the last decade, studies of the thick disk kinematics have been numerous and its velocity distribution is now well constrained. This allows to make a rigourous classification of the stellar populations by computing the probability of each star of a sample to belong to the thin disk and the thick disk according to its velocity (U,V,W), having first evaluated the selection biases which affect local samples. It is the way we have adopted in this paper to investigate the behaviour of Mg and Si (as alpha-element) and Ni (as iron-peak element) in the thin and thick disks. In the past years, Mg abundances in stars have been determined in a wide range of metallicities and temperatures, either in LTE approximation or through detailed NLTE calculations. In the visible range, there are about 10MgIlines, but mostofthem arestronglineswith equivalent width (EW) greater than 200 må, and for such lines NLTE effects can be significant. As shown by several authors (Zhao & Gehren 2000; Idiart & Thévenin 2000; Shymanskaya et al 2000), departures from LTE are observed in MgI lines in stars of various spectral types. We have therefore used the NLTE approach to determine Mg abundances more exactly. The lines that we have used for the other elements have small or moderate EW and considering our temperature and metallicity interval, we have negleted NLTE corrections. As a matter of fact, in the Sun the abundance deviations due to NLTE effects are generally small and do not exceed 0.1 dex, as established in former calculations: 0.07 dex according to Shchukina & Bueno(2001), <0.10 dex according to Gehren et al(2001a, 2001b)forFeIlines,and-0.01dexforSiIlines(Wedemeyer 2001). Gratton et al (1999) and Thévenin & Idiart (1999) have investigated NLTE effects for iron lines in dwarfs and giants of different metallicity. In both papers, the dominating NLTE effect for Fe is the overionization by ultraviolet radiation (UV photons). Gratton et al (1999) found that NLTE corrections may be neglected in most cases, including the stars on the main-sequence and red giant branch. Thévenin & Idiart (1999) derived metallicity corrections of about 0.3 dex for stars with [Fe/H] < 3.0, however for stars with [Fe/H] > 1 this value does not exceed 0.1 dex. Abundances of Fe, Si and Ni were therefore determined in LTE approximation. This paper presents the determination and analysis of kinematical parameters and abundances of Mg, Si and Ni for174dwarfstarsin the domain of 1.0 <[Fe/H] < This study was carried out using a homogeneous spectral material, uniform methods of treatment and NLTE calculations for 9 MgI lines. Two subsamples have been defined on the basis of kinematics to be representative of the thin disk andthe thick diskand areused tohighlightthe chemical behaviour of the two stellar populations. 2. Observations and reduction All the spectra used in this paper are extracted from the most recent version of the library of stellar spectra collected with the ELODIE echelle spectrograph at the Observatoire de Haute-Provence by Soubiran et al (1998) and Prugniel & Soubiran (2001). The performances of the instrument mounted on the 193cm telescope, are described in Baranne et al (1996). ELODIE is a very stable instrument, build to monitor radial velocity oscillations of stars with exoplanets, at a resolving power of in the wavelength range λ λ ÅÅ. Spectrum extraction, wavelength calibration and radial velocity measurement have been performed at the telescope with the on-line data reduction software while straightening of the orders, removing of cosmic ray hits, bad pixels and telluric lines were performed as described in Katz et al (1998). All the spectra of the library have been parametrized in terms of (T eff, logg, [Fe/H] ), either from the literature or internally with the TGMET code (Katz et al 1998).

3 Mishenina et al.: [Mg/Fe] vs kinematics in the disk 3 This allowed us to select a set of 174 FGK dwarfs and subgiants (3.7<logg<4.6) for this study, with metallicities in the range 1 <[Fe/H] < Several solar spectra taken on the Moon and asteroids were also considered as references. The selected spectra have a signal to noise ratio (S/N) at 5500 Å ranging from 100 to 300. The continuum level drawing and equivalent width measurements were carried out by us using DECH20 code (Galazutdinov, 1992). Equivalent widths of lines were measured by Gaussian function fitting. Their accuracy was estimated by comparing our measurements on solar spectra to those obtained by other authors. The mean difference with Edvardsson et al (1993) is : EW Edv EW our = 2.3(σ = 2.3) må for 27 lines of FeI, FeII, SiI and NiI in common; and with Reddy et al (2002) : EW Red EW our = 1.2(σ = 2.1) må for47 lines in common. The comparison is shown in Fig. 1. EW(lit) EW(our) Fig.1. Comparison of solar EW measured in this study with those from the literature (Edvardsson et al 1993 filled circles, Reddy et al 2002 open circles) excitation potentials, calibrated from 45 reference stars in common with Alonso et al (1996), Blackwell & Lynas- Gray (1998) and di Benedetto (1998). The zero-point of the temperature scale was directly adjusted to the Sun, based on 11 solar reflection spectra taken with ELODIE, leading to the uncertainty in the zero-point of about 1K. Judging by the small scatter in our final calibrations and T eff, the selected combinations are only weakly sensitive to effects like rotation, micro- and macroturbulence, non- LTE and other. The application range of the line-depth method is 0.5 <[Fe/H] < For the most metal-poor stars, T eff was determined earlier (Mishenina & Kovthyukh 2001). The H α line-wing fitting was used for stars studed in this work. None of the stars from Mishenina & Kovthyukh (2001) had their temperature measured by line depth ratio because they are too metal poor. In order to testify that the temperature scales adopted in Mishenina & Kovthyukh (2001) and Kovtyukh et al (2003) are consistent, we show in Fig. 2 our adopted temperatures versus those estimated by Alonso et al (1996), Blackwell & Lynas-Gray (1998) and di Benedetto(1998). Our temperature scale is slightly hotter than their by 20 30K, as mentioned in Kovtyukh et al (2003), but the dispersion is satisfactory ( 80K, Tab. 1). We have checked the 2 outliers HD and HD We are certain that Alonso et al s temperature for HD (5483K) is in error because recently Heiter & Luck (2003) determined a temperature of 6000K similar to ours (5932K), a value confirmed by the H α profile. HD is a variable, active, rotating spectroscopic binary rending its photometric measurements suspicious. On the basis of IR photometry, Alonso et al (1996), Blackwell & Lynas-Gray (1998) and di Benedetto (1998) determined respectively 4978K, 4983K, 4937K. Our determination of 5314K agrees well with that of Zboril & Byrne (1998) who find T eff =5260K. 3. Atmospheric parameters In order to perform reliable abundance determinations, it is crucial to derive accurate stellar parameters, especially effective temperatures T eff. V-K or IR flux methods are often advocated as the best (Asplund 2003) but require homogeneous infrared observations which are not available for all stars. Recently our group has improved the line-depth ratio technique pioneered by Gray (1994), leading to high precision T eff for most of our program stars (Kovtyukh et al 2003). This method, relying on ratios of the measured central depths of lines having very different functional dependences on T eff, is independent of interstellar reddening and takes into account the individual characteristics of the star s atmosphere. Briefly, a set of 105 relations was obtained, the mean random error of a single calibration being 60 70K(40 45Kin the most and K in the least accurate cases). The use of calibrations per spectrum reduced the uncertainty to 5 7 K. These 105 relations were obtained from 92 lines, 45 with low (χ <2.77 ev) and 47 with high (χ >4.08 ev) Fig. 2. Comparison of our adopted temperatures with those estimated by Alonso et al (1996), Blackwell & Lynas-Gray (1998) and di Benedetto (1998). Crosses indicate stars with [Fe/H] <-0.5. The two outliers HD and HD are discussed in the text. HD was measured by Alonso, Blackwell and Di Benedetto, that is why there are 3 points for the same star. Finally, Fig. 3 shows a general lack of correlation between the Fe and Mg abundances and T eff that testifies to a correct choice of the effective temperatures.

4 4 Mishenina et al.: [Mg/Fe] vs kinematics in the disk Table 1. Mean difference and standard deviation between our temperatures and those estimated by Alonso et al (1996), Blackwell & Lynas-Gray (1998) and di Benedetto (1998) for 40 stars in common. Reference N T eff (K) σ(k) Alonso et al (1996) Blackwell & Lynas-Gray (1998) di Benedetto (1998) Table 2. Number of stars in common, mean difference and standard deviation of parameter comparison with 5 other studies. N T eff σ log g σ [Fe/H] σ Ref * 1 Edvardsson et al (1993); 2 Gratton et al (1996); 3 Fuhrmann (1998); 4 Chen et al (2000); 5 Fulbright (2000); 0.0 [Fe/H] [Mg/H] Teff Teff Fig.3. Abundances of Fe and Mg vs T eff the wings-fitting for the MgIb lines. For all stars the difference does not exceed 0.2 dex (such difference was detected only for a few program stars). Microturbulent velocities V t were determined by forcing the abundances determined from individual FeI lines to be independent of equivalent width. Starting with V t =1.1 kms 1, we varied it until abundances computed for FeI lines (20 må< EW < 150 må) and plotted as a function of EW-s showed a zero slope. The precision of the microturbulent velocity determination is 0.2 kms 1. The abundance of iron relative to solar one [Fe/H] =log(fe/h) log(fe/h) was used as the metallicity parameter of a star and was obtained from FeI lines. The adopted parameters of the target stars are given in Table 8. To check the adopted model atmosphere parameters, we compared our values with those derived by Edvardsson et al (1993), Gratton et al (1996) Fuhrmann (1998), Chen et al (2000) and Fulbright (2000). The mean differences (other us) and standard deviations are given in Table 2. On the whole, the agreement is rougly good. The two most commonly used techniques to determine surface gravities are the ionization balance of neutral and ionized species and the fundamental relation which expresses the gravity as a function of mass, temperature and bolometric absolute magnitude deduced from parallaxe. For the latter method, stellar masses can be estimated from evolutionary tracks, but metallicities and eventual α- enhancement have to be known first. We therefore choosed to derive log g from ionization balance, a method which might be affected by NLTE effects. However, a detailed study of surface gravities derived by different procedures was preformed by Allende Prieto et al (1999) who concluded that astrometric and spectroscopic (iron ionization balance) gravities were in good agreement in the metallicityrange 1.0 <[Fe/H]< +0.3.Wecomparedouradopted surface gravities to those determined astrometrically by Allende Prieto et al (1999) and obtain a mean difference and standard deviation of and 0.15 respectively for 39 stars in common. This is consistent with an accuracy of 0.1 dex on our estimated spectroscopic gravities. As an additional criterion of the reliable choice of logg we used 4. Abundance determination 4.1. Fe, Si, Ni abundances Fe, Si, Ni abundances were determined from an LTE analysis of equivalent widths using the WIDTH9 code and the atmosphere models by Kurucz(1993). Appropriate models for each star were derived by means of standard interpolation through T eff and logg. The model metallicities were taken with an accuracy of ±0.25 dex around the TGMET first guess. Abundances of the investigated elements Fe, Si, Ni were carried out from a differential analysis relative to the Sun. Solar abundances of Fe, Si and Ni were calculated with our solar EWs and the oscillator strengths log gf from Kovtyukh & Andrievsky (1999), who derived loggf for 565 lines of the 27 chemical elements using solar spectra. We determined : log A(Fe) = 7.55, log A(Si) = 7.55, log A(Ni) = 6.25, where log(h) = 12 (according to Grevesse & Sauval (1998) log A(Fe) = 7.50, log A(Si) = 7.55, log A(Mg) = 7.58, log A(Ni) = 6.25). Depending on authors, the solar iron abundance varies from7.44to7.64(asplund etal2000;blackwelletal1995;

5 Mishenina et al.: [Mg/Fe] vs kinematics in the disk 5 Table 3. List of Fe, Si and Ni lines, atomic parameters and solar EW. Only available at the CDS Gehren et al 2000a; Grevesse & Sauval 1998; Shchukina & Bueno 2001, etc). This disagreement is caused by several factors: the systems of oscillators strengths, the models of solar atmosphere (empirical, hydrodynamical, LTE or NLTE assumptions), the solar spectrum LTE or NLTE synthesis etc. Using laboratory log gf, Blackwell et al (1995) have obtained a solar iron abundance of From calculated and laboratory log gf, Gehren (2001a) found an abundance 7.50 and 7.56 from Fe II lines, and 7.47 and 7.56 from Fe I lines. The detailed hydrodynamical models of the Sun result in the solar value of iron abundance of 7.44 (Fe I lines) and 7.45 (Fe II lines) (Asplund et al ). In the same time of the meteoretic value of the iron abundance is 7.50 (Grevesse & Sauval 1998). The same value was found for the solar atmosphere by Shchukina & Bueno (2001) using the NLTE approximation and hydrodynamical solar model. The solar silicon abundance obtained within the framework of the hydrodynamical models (Asplund 2000) is sligtly lower(7.51) than the meteoric silicon abundance (7.55) (Grevesse & Sauval 1998). To reduce the influence of the used systems of oscillators strengths, atmosphere models and instrumental errors we used a differential method of abundance determination. The list of Fe, Si and Ni lines, atomic parameters and solar EW are given in Table 3. Relative abundances of Fe, Si and Ni are given in Table NLTE calculations for Mg The determination of Mg abundance was carried out through detailed NLTE calculations using equivalent widths of 4 lines (λ λ 4730, 5711, 6318, 6319 Å) and profiles of 5 lines (λ λ 4571, 4703, 5172, 5183, 5528 Å). NLTE abundances of Mg were determined with the help of a modified version of the MULTI code of Carlsson (1986) described in Korotin et al. (1999a 1999b). We have included in the code opacity sources from the ATLAS9 program (Kurucz 1992). This enables a much more accurate calculation of the continuum opacity and intensity distribution in the UV region which is extremely important in the correct determination of the radiative rates of b f transitions. For b b transitions include only continuum opacity and for b f transitions include continuum both opacities and line opacities from ATLAS9. A simultaneous solution of the radiative transfer and statistical equilibrium equations has been performed in the approximation of complete frequency redistribution for all lines. We employed the model of magnesium atom consisting of 97 levels: 84 levels of Mgi, 12 levels of Mgii and a ground state of Mgiii. They were selected from works of Martin & Zalubas(1980) and Biemont & Brault (1986). A detailed structure of the multiplets was ignored and each Fig. 4. Grotrian diagram for MgI. Displayed are the radiative transitions treated explicity in the non-lte calculations E(1/cm) Mg II 1S 1P' 1D 1F' 3S 3P 3P' 3D 3F' LS multiplet was considered as a single term. The fine structure was taken into account only for a level 3s3p 3 P 0, since it is closely connected to the most important transitions in magnesium atom ( λ λ ÅÅ, ÅÅ, 5167 Å, 5172 Å ). Within the described system of the magnesium atom levels, we considered the radiative transitions between the first59levelsof Mgiandgroundlevelof Mgii.Transitions between the rest levels were not taken into account and they were used only in the equations of particle number conservation. Only transitions with λ < Å were selected for the analysis. All 424 b b transitions were included in the linearization procedure. Grotrian diagram for the MgI are displayed in Fig. 4. Photoionization cross-sections were mainly taken from the Opacity Project (Yan et al 1987) keeping a detailed structure of their frequency dependence, including resonances. For some important b f transitions, the crosssection structure is extremely complicated making it difficult to describe it using only simple approximations like ν 3. Oscillator strengths were selected from the extensive compilative catalogue by Hirata & Horaguchi (1994). Some information was obtained through the Opacity Project. As we ignored a multiple structure of all the levels, the oscillator strengths for each averaged transition were calculated as f = gif i gi.

6 6 Mishenina et al.: [Mg/Fe] vs kinematics in the disk After the combined solution of radiative transfer and statistical equilibrium equations, the averaged levels have been split with respect to multiplet structure, then level populations were redistributed proportionally to the statistical weights of the corresponding sublevels and finally the lines of the interest were studied. Electron impact ionization was described using Seaton s formula (Seaton 1962). Collision rate between ground level Mgi and ground level of Mgii were approximated by use of the fits from Voronov (1997). For electron impact excitation for all allowed b b transitions we used the van Regemorter (1962) formula. Collisional rates for the forbidden transitions were calculated by using the semiempirical formula(allen 1973), with a collisional force of = 1. Inelastic collisions with hydrogen may play a significantroleinthe atmospheresofcoolstars.we tookinto account this effect with the help of Drawin s formulas(1968) offered Steenbock & Holweger (1984), with a correction factor of 1/3. Nevertheless, as it was shown in our test calculation, the influence of such collisions is negligible for the stars considered in our paper. In particular, the variation of the correction factor from 0 to 1 results in a change of the equivalent width less than 0.5%. For all the transitions we also took into account such broadening parameters of lines as radiative damping, Stark effect, van der Waals damping and microturbulent velocity. For temperatures, considered by us, the influence of Stark effect is small. The main influence on the profile was exerted by van der Waals broadening. The Unsöld s (1955) formula which is widely used to allow for this effect, is known to yield somewhat understimated coefficients C 6. To refine these coefficients, we conducted a special comparison of the observed and calculated line profiles in solar spectrum using method Shimanskaya, Mashonkina & Sakhibullin (2000). The Solar Flux Atlas of Kurucz et al. (1984) and Kurucz s model of the solar atmosphere (1992) were used. To take into account the chromospheric growth of temperature, this model was completed by a model of the solar chromosphere from work of Maltby et al.(1986). Nevertheless, the influence of a chromosphere on the equivalent width of considered lines was insignificant (less than 1%). The corrections log C 6 to the classic Unsöld (1955) collisional damping constant derived from solar line wing fitting are given in Table 4. These values are comparable for common lines with results obtained in the paper of Shimanskaya, Mashonkina & Sakhibullin (2000). Similar values were obtained by Barklem, Piskunov & O Mara (2000). For example, their correction factors for the lines λλ 5172, 5183 ÅÅ are log C 6 =0.85. Oscillator strengths for observed lines were selected from the compilative catalogue by Hirata & Horaguchi (1994). The estimates obtained from the NLTE profile analysis of the MgI lines in the solar spectrum give the solar magnesium content (Mg/H)=7.57±0.02 which agrees well with the determination of Grevesse & Sauval (1998) and Table 4. Atomic data for the Mg lines Lambda χ (ev) loggf C 6 C 6 (Barklem et al 2000) Shimanskaya, Mashonkina & Sakhibullin (2000), who derived (Mg/H)=7.58. The comparison the NLTE profiles with those observed are presented in Fig. 5. NLTE abundances ofmg are given in Table 8. The NLTE effects for the used lines appeared to be very small (less than 0.05 dex). The similar results were obtained in Shimanskaya, Mashonkina& Sakhibullin (2000). Relative flux Relative flux Relative flux λ λ Relative flux Relative flux Fig. 5. Comparison of synthetic NLTE profiles (continuous line) of the MgI lines to the observed spectrum of HD (dots) λ λ λ

7 Mishenina et al.: [Mg/Fe] vs kinematics in the disk 7 Table 5. Effects of uncertainties of model parameters and EW measurements on the derived abundances HD /4.4/1.1/ 0.24 El T eff logg Vt EW [Fe/H] Total FeI FeII MgI SiI NiI HD /4.5/1.1/ 0.77 FeI FeII MgI SiI NiI HD /4.0/1.0/-0.03 FeI FeII MgI SiI NiI (1) T eff = -100 K; (2) logg = +0.2; (3) V t=+0.2 kms 1 ; (4) [Fe/H] = 0.25 dex; (5) EW= ± 2mÅ ; (6) Total error. 5. Error analysis Average standard deviations of abundances obtained from lines of Fe I (the number of used lines differs from starto star),11-24lines ofsi Iand12-17linesofNiIare 0.10, 0.08 and 0.08 respectively. Final errors in abundances result mainly from errors in the choice of parameters of the model atmospheres and in equivalent width measurements (gaussian fitting, placement of the continuum) in the case of Fe, Si, and Ni or fitting a synthetic spectrum in the case of Mg. Table 5 lists the errors obtained when changing the atmospheric parameters by T eff = -100 K; logg =+0.2; V t =+0.2 kms 1 ; [Fe/H] =±0.25 and by assuming an uncertainty of ±2 må in the EW. These values were adopted taking into account the intrinsic accuracy of the atmospheric parameter determination, the processing of the spectra and the comparison of our parameter definition with those of other authors. This test has been performed for 3 stars with different characteristics. As seen in Table 5, the total uncertainty reaches 0.12 dex for iron abundance determined from II species and 0.10 dex in the case of I species. The standard deviation obtained by comparing our[fe/h] determinations to those obtained by other authors (Table 2) shows that we are consistent with them at a level lower than 0.10 dex. For Mg, Si and Ni, the total uncertainty due to parameter and EW errors is 0.08 dex, 0.05 dex and 0.09 dex respectively. We have also compared our Mg, Si and Ni Table 6. Abundance [Mg/H], [Si/Fe], [Ni/Fe] comparison with other authors N Mg σ Si σ Ni σ Ref * 1 Idiard & Thévenin (2000); 2 Gratton et al (2003); 3 Edvardsson et al (1993); 4 Chen et al (2000); 5 Reddy et al (2003) abundances with those determined by other authors. For [Mg/H], we have carried out the comparison with Idiart & Thévenin (2000) and Gratton et al (2003) who have also used the NLTE approach. Our [Si/Fe] and [Ni/Fe] determinations have been compared to those of Edvardsson et al (1993), Reddy et al (2002) and Chen et al (2000). The results of comparisons are listed in Table 6. The mean differences are lower than 0.05 and the standard deviations lower than 0.1 dex for all elements proving the good consistency of our study with previous works relying on similar methods. 6. Stellar kinematics, metallicity, elemental abundances Alltheselectedstarsarebrightandnearbyenoughtohave parallaxes and proper motions measured by Hipparcos (ESA 1997). These astrometric quantities have been combined with radial velocties measured on the ELODIE spectra by cross-correlation(with an accuracy better than 100 ms 1 ) to compute the 3 components (U,V,W) of the spatial velocities with respect to the Sun. Combining the measurement errors on parallaxes, proper motions and radial velocities, the resulting errors on velocities are of the order of 1 kms 1. Our first concern was to analyse the content of our sample in terms of stellar populations. As we were interested in the characterization of abundance patterns in the thindiskandthethickdisk,wehaveselectedoursampleto span the metallicity range 1.0 <[Fe/H] < +0.3 in order to define 2 subsamples representative of these two populations. A discrimination of thin disk and thick disk stars is possible using the fact that the two disks are known to be distinct by their spatial distribution and local density, and by their velocity, metallicity and age distributions. However we have performed this classification with a pure kinematical appoach because velocity distributions are the best constrained from observations reported in the literature. If the thin disk velocity distribution is well known thanks to Hipparcos data, authors do not yet agree on the properties of the thick disk. A review of the current knowledgeofthethickdiskisgiveninnorris(1999).there is still a controversy between the adherents of a flat and dense thick disk, with typical scale height of 800pc and

8 8 Mishenina et al.: [Mg/Fe] vs kinematics in the disk Table 7. Kinematical parameters adopted for the thin disk and the thick disk from Soubiran et al (2003) thin disk thick disk V lag (kms 1 ) σ U (kms 1 ) σ V (kms 1 ) σ W (kms 1 ) local relative density of 6-7% (Reylé & Robin 2001) to 15% (Soubiran et al 2003), and the adherents of a thick disk with a higher scale height, typically 1300pc and a lower local relative density, of the order of 2% (Reid & Majewski 1993). Velocity dispersions are generally found to span typical values between 30 and 50 kms 1, with an asymmetric drift between -20 and -80 kms 1. Recently, Soubiran et al (2003) determined the kinematical parameters of the thick disk and the thick disk from an unbiased sample of clump giants. These values, listed in Tab. 7, are used in the present study. The metallicity distributions of the thin and thick disks are also matter of debate. Recent thick disk metallicity determinations quoted mean values from -0.36(Bell 1996) to (Robin et al 1996, Gilmore et al 1995), passing by (Soubiran et al 2003). Morrison et al (1990) brought to the fore low-metallicity stars ( 1.6 <[Fe/H] < 1.0) with disk-like kinematics. Chiba & Beers (2000) estimate that 30% of the stars with 1.6 <[Fe/H] < 1.0belong to the thick disk population. Similarly, Bensby et al (2003) found in their sample stars with super-solar metallicities and thick disk kinematics. It remains unclear whether these extreme populations are separate from the thick disk or their metal-weak and metal-rich tails. Moreover, due to the great overlap of the thin disk and thick disk distributions, it is very difficult to estimate where the transition occurs. Authors involved in the study of metallicity distributions are often concerned with the chemical evolution of the Galaxy and generally do not attempt to separate the thin and thick disks but rather consider the thick disk as an integral part of the disk. It was the case for instance in the work of Edvardsson et al (1993) on the age - metallicity relation. Haywood (2001) performed a revision of the solar neighbourhood metallicity distribution and concluded that it peaks at [Fe/H] =0 and that only 4% of the nearby stars have [Fe/H] < These considerations made us think that the selection of thin disk and thick disk stars on a pure kinematical criterion would be more reliable. Fig. 6 showsthe distribution ofthe samplein the plane metallicity - W velocity. It is clear from this plot that a transition occurs at [Fe/H] : above this value we measure a dispersion of the W velocities of 16 kms 1 while below this value the dispersion is 38 kms 1. These values are typical of the thin disk and the thick disk respectively (see Tab 7). They are related to the different scale heights of the two populations and testify that our sample is indeed a mixture of thin and thick disk stars. Fig.6. Distribution of the 174 stars in the plane metallicity - W velocity. In order to compute the probability of each star to belong to the thin disk or the thick disk on the basis of its spatial velocity, we need to know the kinematical parameters of the two disks and their proportion in our sample. The kinematical parameters are listed in Tab. 7, but the proportionsinoursampleareofcoursedifferent ofthereal proportion in the solar neibourhood because our selection in the range 1 <[Fe/H] < +0.3 is supposed to have increased the number of thick disk stars comparatively to their relative density in the solar neighbourhood (2% to 15%). Moreover the ELODIE library was build from various observing programs which may have biased its content and consequently our sample. To estimate the proportion of thin and thick disk stars in our sample, we have applied an algorithm of deconvolution of multivariate gaussian distributions, on (U,V,W) velocities. This algorithm, SEM (Celeux & Diebolt 1986), solves iteratively the maximum likelihood equations, with a stochastic step, with no a priori information on the mixture of the gaussian distributions. The results are unambiguous, our sample is consistent with a mixture of 2 gaussian populations with parameters very similar to those listed in Tab. 7, but with relative densities of 75% and 25% respectively for the thin disk and the thick disk. Accordingly, we have computed the probability of each star, with a measured velocity (U,V,W), to belong to the thin disk (Pr 1 ) and to the thick disk (Pr 2 ) : F i = F 1 Pr 1 = p 1 F, Pr F 2 2 = p 2 F 1 3 exp 0.5[ U 2π σui σ Vi σ Wi σu 2 + (V V i) i σv 2 + W i σw 2 ] i F = p 1 F 1 +p 2 F 2 where p 1 = 0.75,p 2 = 0.25 are the relative densities of the 2 populations in our sample and V i,σ Ui,σ Vi,σ Wi,i = 1,2 are the kinematical parameters of the thin and thick disks listed in Table 7. From these probabilities we have selected two subsamples kinematically representative of each population in order to highlight their respective properties. Putting the

9 limit at Pr i > 80% ensures a minimal contamination of each subsample by the other population. The thin disk subsample includes 109 stars whereas the thick disk subsample include 30 stars. In the figures, unclassified stars are represented by small dots, thick disk stars by larger black dots, thin disk stars by open dots. The only halo star of the sample, HD ([Fe/H] =-1.21, V= kms 1 ) isnot representedin the plots forabetter clarity. Fig. 7, 8 and 9 show the trend of the abundances of Mg, Si and Ni as a function of [Fe/H]. It can be seen that : Mishenina et al.: [Mg/Fe] vs kinematics in the disk 9 there are a few thin disk stars with [Fe/H] < 0.30 metal poor stars ([Fe/H] < 0.60) are all enriched in Mg and Si ([Mg/Fe] > +0.20, [Si/Fe] > +0.15) on the contrary at solar metallicity, the enrichment of Mg and Si does not exceed at a given metallicity thick disk stars have higher [Mg/Fe] on average than thin disk stars the dispersion of [Si/Fe] is remarkably small for [Fe/H] > 0.30 but quite high at lower metallicity [Mg/Fe] and [Si/Fe] decline with metallicity from about +0.40,+0.35 to 0.0, respectively there are stars with thick disk kinematics at solar metallicities, their abundance trends follow the thin disk a rise and upturn is visible for [Ni/Fe] for [Fe/H] > 0 Fig.8. Same as Fig. 7 for Si. Fig.9. Same as Fig. 7 for Ni. All these features are in very good agreement with those described in Bensby et al (2003). Fig. 10. Correlation between [Mg/Fe] and [Si/Fe] Fig.7. [Mg/Fe] vs [Fe/H] for the whole sample. Black dots indicate thick disk stars, open dots thin disk stars, small dots represent the unclassified stars. Fig. 10 represents the correlation between Mg and Si relative to iron. Mg and Si are α elements which are supposedtobemainlyproducedinsnii. Itcanbeseeninthis plotthatatransitionoccursat[mg/fe] +0.2.Starswith [Mg/Fe]< +0.2haveameanabundanceof[Si/Fe]=+0.08 with a very low dispersion of ±0.03, lower than our error estimates. [Mg/Fe] is more dispersed (±0.06) around the mean of [Mg/Fe] = For stars with [Mg/Fe] > +0.2 the distribution is consistent with a linear correlation : [Si/Fe] =0.7[Mg/Fe] (rms=0.06). 7. Discussion 7.1. Stars with thick disk kinematics and thin disk metallicity A remarkable feature in the Fig. 7 to 10 is the behavior of stars with thick disk kinematics and thin disk metallicity. They follow exactly the chemical trend of the thin disk. This led us to wonder if these stars were real thick disk stars. According to the W vs [Fe/H] distribution shown in Fig. 6, metal-rich stars are expected to have a smaller scale height than metal poor stars. To have confirmation, we have computed the orbital parameters of our stars, by integrating the equations of motion in the galactic model

10 10 Mishenina et al.: [Mg/Fe] vs kinematics in the disk of Allen & Santillan (1991) over an age of 8 Gyr. The adopted velocity of the Sun with respect to the LSR is (9, 5, 6) kms 1, the solar galactocentric distance R = 8.5 kpc and circular velocity V lsr = 220 kms 1. Table 8 lists the spatial velocities and the parameters of the orbits : apogalactic and perigalactic distances (R min,r max ), the maximumdistancefromtheplane(z max )andeccentricity (e). The distribution of Z max vs [Fe/H] in Fig. 11 confirms the small scale height of stars with thick disk kinematics and solar metallicity, ruling them out of the thick disk. They have been classified into the thick disk because of they have a large eccentricity, but they have a scale height which is typical of the thin disk. Thus the origin of these stars cannot be in the thick disk neither in the local thin disk despite the same chemical behaviour. The origin of these stars has to be found elsewhere.the only ambiguous star, HD having [Fe/H] =+0.12 and Z max =1 kpc, appears as an outlier in the distribution. This star is discussed farther. to demonstrate if they belong or not to the thick disk we have calculated the orbits of their stars and classified them with the same criterionas ours (Pr i > 80%). The distribution of Z max versus [Fe/H] is shown on Fig. 12. Globally their distribution is similar to ours despite their sample is much smaller. There is a sudden rise in Z max at [Fe/H] < 0.3. They have also 3 outliers with solar metallicities and Z max 1 kpc, HD being in common with us. The five other stars with thick disk kinematics and [Fe/H] > 0.30 have a low scale height typical of the thin disk in agreement with our findings. The question wether HD190360, HD and HD are really part of the thick disk is an open question because they appear as outliers in the distribution of kinematics vs metallicity. It is only with complete samples with no kinematical bias that the relation of such stars with the thick disk will be clarified. Fig.11. The maximum distance to the plane of orbits, Z max, versus [Fe/H] Fig.12. The maximum distance to the plane of orbits, Z max, versus [Fe/H] for Bensby et al s sample. The classification thin disk / thick disk is the same as ours (Pr i > 80%). Among the 8 stars that verify Pr 2 > 80% and [Fe/H] > 0.30, HD030562, HD139323, HD are K dwarfs with very similar velocities and metallicities (Tab. 8, [Fe/H] HD has been identified to be part of the HR 1614 moving group by Feltzing & Holmberg (2000). It is very likely that we have discovered 2 new candidates of this moving group whose age is estimated to be 2 Gyr. The origin of this group is however unknown. HD010145, HD135204, HD have U > 80 kms 1. An hypothesis is that they have been thrown out from the inner disk by the galactic bar. The over density of metal-rich stars with excentric orbits confined near the plane have have been interpreted as the signature of the galactic bar in the solar neighbourhood by Raboud et al (1998) for instance. Bensby et al (2003) have also found several stars with thick disk kinematics and thin disk metallicity on the basis of a similar probability calculation than ours, but using the ratio of Pr 2 to Pr 1. If in our case these stars have a clearly a flatter distribution than the thick disk, it was not mentioned in Bensby et al s study. As it is important 7.2. The thin disk Several thin disk stars are found at low metallicity, down to [Fe/H] = Several previous studies have also established the existence of such stars, however the exact contribution of metal-poor stars to the thin disk is difficult to evaluate from uncomplete samples. Reddy et al (2002) have found in their samples a significant number of stars with [Fe/H] < 0.40 that they identify as belonging to the thin disk. They have analysed 181 FG dwarfs, a sample size similar to ours. Their sample spans the metallicity interval 0.70 to +0.15, with a peak at [Fe/H] = Only 3 stars were identified as belonging to the thick disk, based onthe criterion[mg/fe] > +0.20,all the otherstars being identified as thin disk stars. The peak at [Fe/H] = 0.40 shows clearly that they have favoured in some way the metal-poor tail of the thin disk. Unfortunatly their selection criteria are not described. On the contrary, the way Chen et al (2000) have favoured metal-poor thin disk stars in their sample is very clear : they have selected

11 Mishenina et al.: [Mg/Fe] vs kinematics in the disk 11 equal number of stars in every metallicity interval of 0.1 dex in the range 1.0 <[Fe/H] < According to the low density of the thick disk in the solar neighbourhood, thin disk stars are still dominating at moderatly low metallicity. Haywood (2001) has very carefully taken into accountthe selectionbiases ofhis sample to study the metallicity distribution of the local disk but his metallicities are photometric and no special attention was paid to a kinematical distinction between the thin disk and the thick disk. Nevertheless he found, contrary to Wyse & Gilmore (1995), a low contribution of the thin disk at [Fe/H] < Our observations suggest that the distribution of [α/fe] in the thin disk is very narrow, specially for Si, at [Fe/H] > On this point we are in perfect agreement with Reddy et al (2002) andbensby et al (2003). For the metal-rich part, we obtain as mean values and dispersions : [Mg/Fe] =+0.05, σ [Mg/Fe] =0.07, [Si/Fe] =+0.07, σ [Si/Fe] =0.03.Suchanarrowchemicaldistributionimplies that the stars formed from an homogeneous gas. At lower metallicity, the Mg enhancement is higher. However at a given metallicity the thin disk shows a lower Mg enhancement than the thick disk. The behaviour of Si is quite surprising, showing high dispersion, whereas Bensby et al (2003) find a similar behaviour than Mg. We notice also that the Sun would be slightly Si-poor as compared to the thin disk. However this offset stays within the error bars and might be related to the adopted zero point of our abundance scales. Many previous studies (Edvardsson et al 1993, Chen et al 2000, Reddy et al 2002, Bensby et al 2003) exhibit also a positive offset in [Si/Fe] of the thin disk. A rise of [Ni/Fe] from 0.1 at [Fe/H] =-0.20 to +0.1 at [Fe/H] =+0.20 is visible in Fig. 9 as well as an upturn for the most metal rich star. A similar feature is also observed in Bensby et al (2003). which was also mentioned in Prochaska et al (2000, has to be compared to the results of Fuhrmann (1998) and Gratton et al (2000) who found in the thick disk a constant overabundance of α-elements with respect to iron, with a value similar to that of the halo. In Fuhrmann (1998), only thin disk stars exhibit a decrease of [Mg/Fe] with increasing[fe/h]. We believe that we havebeen able to detect this trend, as Feltzing et al (2003) and Bensby et al (2003), because we have used a more precise criterion to isolate a sample a pure thick disk stars. We agree with most of the previous works that a discontinuity exists between the chemical distributions of the thin and thick disks. Our data suggest that the star formation in the thick disk stopped when the gas had the composition [Fe/H] = 0.3, [Mg/Fe] =+0.20, [Si/Fe] = The distribution of Z max vs [Fe/H] (Fig. 11) does not exhibit a clear vertical gradient in the thick disk metallicity. We have searched for a vertical gradient in the α- element abundances. Fig. 13 represents Z max vs [α/fe] =0.5([Mg/Fe] +[Si/Fe] ) for the thick disk (Pr and [Fe/H] -0.3). A gradient may exist, but should be confirmed with a larger sample. The existence of a vertical gradient would imply a significant timescale for the formation of the thick disk. This implication is also true if the decline of α enhancement is interpreted by the onset of SNIa The thick disk Having eliminated from the thick disk the 8 stars having high metallicity and a flat distribution (Sect. 7.1), we see from Figure 7 that [Mg/Fe] decreases from the halo value (+0.40)at [Fe/H] =-1.00to+0.20at [Fe/H] =-0.30.Si has the same behaviour, decreasing from to +0.17, with one deviating star HD having a lower enrichment in Si. These findings are nicely consistent with the trends found by Bensby et al (2003) that they interpreted as the signature of the chemical enrichment by SNIa. Feltzing et al (2003) claimed the existence of a knee in the [α/fe] abundance trends tracing the beginning of the contribution of the SNIa. This knee may also exist in our data, especially in [Mg/Fe] but the dispersion is too high to conclude. Definitively a larger sample of thick disk stars is necessary to clarify this point. Contrary to Bensby et al (2003) and Feltzing et al (2003), our data suggest that the metallicity distribution of the thick disk does not continue at solar metallicity. The trend of declining ehancement, Fig. 13. A possible vertical gradient in [α/fe] in the thick disk? 8. Summary and perspectives We have investigated the chemical behaviour of the thin disk and the thick disk with a large sample of FGK dwarfs spanning the metallicity interval [-1.00 ; +0.30]. We have used an homogeneous set of high resolution, high S/N echelle spectra to derive abundances of Fe, Mg, Si and Ni with detailed NLTE calculations for Mg. All sources of errors have been analysed and we have compared our abundances with several studies taken in the literature, showing a good agreement. Orbits and space velocities were computed and used to estimate the probability of each star to belong to the thin disk and the thick disk, taking into account the selection bias that affects our local sam-