arxiv: v1 [astro-ph.sr] 3 Sep 2012

Transcription

1 Astronomy & Astrophysics manuscript no. Posbic2012 c ESO 2014 May 5, 2014 SPADES: a stellar parameters determination software H. Posbic 1, D. Katz 1, E. Caffau 2,1, P. Bonifacio 1, A. Gómez 1, L. Sbordone 2,1, and F. Arenou 1 1 GEPI, Observatoire de Paris, CNRS, Université Paris Diderot, Place Jules Janssen, 92190, Meudon, France 2 Zentrum für Astronomie der Universität Heidelberg, Landessternwarte, Königstuhl 12, 69117, Heidelberg, Germany Received 14 April 2012/Accepted 16 July 2012 ABSTRACT arxiv: v1 [astro-ph.sr] 3 Sep 2012 Context. As increasingly more spectroscopic data are being delivered by medium- and high-resolving power multi-object spectrographs, more automatic stellar parameter determination softwares are being developed. The quality of the spectra collected also allows the determination of elemental abundances. Aims. SPADES is an automated software for determining: the radial velocity (Vr), the effective temperature (T eff ), the surface gravity (log g ), the metallicity ( ), and most importantly, the individual abundances. In this first version it is targeted on the analysis of mid-f-g dwarfs, but is meant to evolve to analyze any type of single stars. Methods. SPADES relies on a line-by-line modeling to determine the stellar parameters. Results. The internal systematic and random errors of SPADES were assessed by Monte Carlo method simulations with synthetic spectra and the external systematic errors by analysing real ground-based observed spectra. For example, by simulating the Giraffe setups HR13 and HR14B with synthetic spectra for a dwarf with T eff = 5800 K, log g=4.5, =0.0 dex and with a signal-tonoise ratio (snr) of 100, the stellar parameters are recovered with no significant bias and with 1-σ precisions of 8 K for T eff, 0.05 for log g, for, for [Ti/Fe] and 0.01 for [Ni/Fe]. Key words. Methods: data analysis - Techniques: spectroscopic - Stars: fundamental parameters - Stars: abundances 1. Introduction The present and future large spectroscopic surveys are going to significantly increase the number of spectroscopic data to be analyzed. A few examples are RAVE with some stars observed so far (Boeche et al. 2011; Siebert et al. 2011), the Gaia-ESO Survey with about 10 5 stars to be observed, Gaia with about 200 million stars (Katz et al. 2004; Katz 2009). To analyze these quantities of data, automatic spectra analysis softwares are needed. The development of automated spectra analysis software can be traced back to the early 1990s. Since then, many programs were developed based on a wide variety of methods: e.g. orthogonal vector projection (Cayrel et al. 1991; Perrin et al. 1995), minimum distance method, e.g. TGMET (Katz et al. 1998; Soubiran et al. 2003), and ETOILE (Katz 2001), artificial neural networks (Bailer-Jones 2000; Allende Prieto et al. 2000), feature fitting softwares such as Abbo (Bonifacio & Caffau 2003), projection vector softwares such as MATISSE (Recio-Blanco et al. 2006; Bijaoui et al. 2008), principal component analysis softwares such as MAχ (Jofré et al. 2010), and forward modeling algorithm ILIUM (Bailer-Jones 2010). One can also cite DAOSPEC (Stetson & Pancino 2008, 2010) and FITLINE (François et al. 2003), which automatically measure equivalent widths. DAOSPEC and FITLINE are not automated per-say parameter determination softwares, but coupled with programs that derive atmospheric parameters from equivalent widths, they are used by many astronomers to automatically parameterize stars. The work presented here describes the development of a new automatic stellar spectra analysis software. This is based on line-by-line analysis. The reference lines are modeled. One of the motivations to develop SPADES (Stellar PArameters DEtermination Software) is the analysis of a set of 200 Giraffe HR13 and HR14 spectra, mainly F-G dwarfs, that were collected in the context of a study of the structure of the Galactic disk. Therefore, the first version of SPADES is illustrated and tested with simulated Giraffe-like spectra of mid-f-g dwarfs. SPADES is coded with Java. SPADES is described in Sect. 2. It relies on auxiliary data: a line list and a grid of reference spectra. They are presented in Sect. 3. The performances of SPADES were evaluated by Monte Carlo method simulations against synthetic spectra and by comparison to observed spectra of the Sun,νAnd,βVir,µHer, and σ Dra (see Sect. 4). 2. SPADES 2.1. General concept SPADES automatically determine: the radial velocity (Vr), the effective temperature (T eff ), the surface gravity (log g ), the metallicity ( ), and individual abundances ([X/Fe] for the element X) for single stars. This first version of SPADES is focused on mid-f-g dwarfs and the micro-turbulence is set to 1 km s 1. Many of existing automatic spectra analysis softwares uses global methods, meaning they process all pixels of the studied spectrum (Katz et al. 1998; Katz 2001; Soubiran et al. 2003; Recio-Blanco et al. 2006). With SPADES we aimed to explore another possibility, which is the line-by-line analysis. Each parameter is determined using a pre-defined set of lines. The merits 1

2 of this technique are, first, the possibility to select lines that are particularly sensitive to the stellar parameter to be determined, and second, this allows one to exclude lines that are not correctly modeled. A drawback is that only a fraction of the total available information is used. Often, whether automated or manual, methods using a lineby-line analysis rely on equivalent width measurement and use a curve-of-growth analysis. We chose to rely on synthetic spectra modeling and profile fitting. The motivation is twofold. First, this takes care of small blend problems. Assuming that the blend is correctly modeled, its contribution to the observation and to the modeled profile will (partly) cancel out. Second, some parameters can be derived by line profile fitting, e.g. T eff from the wings of the Balmer lines. A drawback of this method is the need for a very large library of synthetic spectra. Reading and processing many synthetic spectra is computationally demanding and therefore the analysis is slower. In the line-by-line analysis the accuracy and precision of the results is very sensitive to the determination of the continuum or pseudo-continuum. Yet, this is often difficult particularly in cool and/or super metal-rich stars. SPADES determines local continuum from the comparison of the reference (synthetic) spectra and the studied spectrum (see sect ). As opposed to mean metallicities, SPADES determines iron and individual elements abundances. To derive of the individual abundances SPADES calculates on-the-fly atmospheric models and synthetic spectra for the values of T eff, log g, it had determined (see sect. 2.6). SPADES can simultaneously analyze various spectral domains (possibly with different resolving powers) of the same target, e.g. different Giraffe setups, or different orders of an echelle spectrograph, or even spectral domains collected with different spectrographs Processing overview A SPADES analysis proceeds in three steps: it first determines the radial velocity, then the atmospheric parameters, and finally it determines the individual abundances. The radial velocity is determined by cross-correlation with a synthetic template (see sect. 2.4). The three atmospheric parameters (effective temperature, surface gravity, iron over hydrogen ratio) are determined independently (see sect. 2.5). Currently, each parameter is associated to one diagnostics. SPADES was conceived to aggregate more diagnostics along its life span and more will be added in the future Current diagnostics are effective temperature: Balmer line wings fitting (suited for T eff derivation in F-G dwarfs), surface gravity: ionization equilibrium, iron abundance: neutral iron line profile fitting. The individual abundances are determined element per element by line profile fitting. The reader will notice that some of these diagnostics are similar to those used in the so-called detailed analysis. Two differences should be emphasized: (i) SPADES relies on profile fitting and modeling rather than equivalent width measurement, and (ii) the SPADES precision is not limited by the step in the grid of atmospheric models and spectra (see Sect ) Pre-processing Each of the three processing steps includes pre-processing tasks. Pre-processing a spectrum includes shifting to rest frame, excluding cosmics, excluding telluric lines, convolving the synthetic spectrum to the resolving power of the ground-based observed spectrum, normalizing to continuum and excluding possible unexposed edges of the spectrum. Every functionality remains an option: it can be switched on and off in the configuration file by the user. Shifting to rest frame is performed by using either a radial velocity provided by the user or that determined by SPADES (see sect. 2.4). The exclusion of cosmics is performed by detecting the fluxes that lie significantly above the continuum and replacing their value by interpolation with the neighboring pixels. The wavelengths of the telluric lines were pre-tabulated. These lines are not cut. Their pixels are simply flagged as not valid and are excluded from all subsequent analysis. The spectral resolution of the synthetic spectra should be matched to the resolution of the ground-based observed spectrum. To do so, the reference spectra are convolved using a Gaussian profile (valid for slow rotators). Flux Normalised Raw spectrum Wavelength (A) Final spectrum Wavelength (A) Fig. 1. Observed Giraffe/VLT (HR14B setup) Sun spectrum before normalization, with the polynomial fit shown as a thick blue line (top) and after normalization (bottom) The synthetic spectra are already normalized to continuum. The studied spectrum is iteratively fitted with a polynomial function, using a sigma-clipping technique to exclude the lines and cosmic rays from the fit. The degree of the polynomial is defined by the user, depending on the large-scale shape of the spectrum. The studied spectrum is normalized to 1 at the level of the continuum by division by the polynomial. Figure 1 presents an observed Giraffe/VLT (HR14 setup) Sun spectrum before normal- 2

3 ization, with the polynomial shown as a thick blue line (top) and after normalization (bottom). Edge exclusion is necessary if the ground-based observed spectrum suffers from underexposed or unexposed pixels. SPADES can automatically detect and flag as not valid the underexposed/unexposed pixels First step: Radial velocity The radial velocity determination is based on a cross-correlation of the studied spectrum and a synthetic template spectrum. The atmospheric parameters of the template spectrum are chosen by the user. At this stage in the processing, SPADES has not yet provided estimates for the atmospheric parameters of the studied source. If no external information is available to guide the choice of the template, it is likely that there will be a mismatch between the ground-based observed spectrum and the selected template. Monte Carlo method simulations were performed to assess the impact of selecting a template 400 K and K warmer than the studied spectrum. The tests were conducted with 200 snr=50 synthetic studied spectra of solar parameters. In both cases (400 and K offsets), the mismatch leads to a bias of about 0.2 km s 1 and no significant change in the precision of the estimated radial velocities (see table 1). A bias of 0.2 km s 1 represents about 1/60 th of the full width at half maximum (FWHM) and has a negligible impact on the following steps of the analysis. If, a posteriori, the parameters estimated by SPADES and those of the selected template are significantly different, the radial velocity can be re-estimated with a better matching template. The studied and the template spectrum are pre-processed. The reference spectrum is then shifted step by step in radial velocity over an interval defined by the user. In the tests, we used an interval of -500 to 500 km/s. At each step, the synthetic spectrum is re-sampled to the same sampling as the studied spectrum. At each step, the correlation coefficient of the two spectra is calculated: CC(v)= 1 N N i=0 (S (i) <S>) (R(i) <R>), (1) stdv(s ) stdv(r) where v is the velocity shift applied to the template spectrum, N is the number of valid pixels in the studied spectrum, S (resp. R) the studied (resp. reference) spectrum counts, < S > (resp. < R>) the mean studied spectrum (resp. reference) counts, and stdv the standard deviation: stdv(x)= ( 1 N N (X(i) < X>) 2) 1/2, (2) i=0 The correlation coefficients as a function of the radial velocity shifts form the cross-correlation function. Figure 2 shows the cross-correlation function obtained for a synthetic studied spectrum of solar parameters (T eff = 5777 K, log g= 4.44 and = 0.0 dex) and S/N=100 (Poisson noise) and a noiseless synthetic template with the same atmospheric parameters. The radial velocity of the observed source corresponds to the maximum of the correlation function. The maximum is derived in two steps. First, the highest correlation coefficient is found and then a second-degree polynomial is adjusted to the correlation function in the vicinity of the highest coefficient. Correlation Coefficient Correlation function Radial velocity (km/s) Fig. 2. Correlation function for a synthetic studied spectrum of solar parameters (T eff = 5777 K, log g = 4.44 and = 0.0 dex) and S/N=100 (Poisson noise) and a noiseless synthetic template with the same atmospheric parameters This allows one to estimate the radial velocity with a precision better than the step used to shift the template. The error on the radial velocity determination,σ Vr is estimated with the Zucker (2003) formula: σ 2 Vr = N C (m) C(m) C(m) 2 1 C(m) 2 1, (3) where C(m) is the value of the correlation function at its maximum, C (s) the value of the second derivative of the correlation function at the maximum of the correlation function and N the number of valid 1 pixels in the ground-based observed spectrum. The performances of the radial velocity estimation of the radial velocity were assessed with a Monte Carlo method. For three S/N: 30, 50 and 100, Giraffe HR13 and HR14B spectra were simulated (200 per SNR) by injecting a synthetic spectrum with Poisson noise. The synthetic spectrum was a solar-like spectrum in the rest frame. The same spectrum, but noiseless, was used as template. Table 1 presents the median and dispersion of the distributions of radial velocities derived by SPADES for each SNR. Because the true radial velocity of the studied spectra is 0 km s 1, the median of a distribution is an estimate of the systematic error and the dispersion an estimate of the random error. 1 Some non-valid pixels affected by, e.g. telluric lines, could be excluded from the derivation of the correlation function (see Sect. 2.3). 3

4 Table 1. Median and dispersion (in km s 1 ) of the radial velocity distribution derived by SPADES for 200 studied spectra for SNR 30, 50 and 100. SNR Vr median σ Vr (km s 1 ) (km s 1 ) Second step: Stellar parameter determination Overview SPADES determines the effective temperature (T eff ), the surface gravity (log g ), and the iron abundance ( ). In this first version targeted on mid-f-g dwarfs, we assumed a microturbulence of 1 km s 1 for all our stars, including the tests with the observed spectra. The determination of the atmospheric parameters is iterative. Initial values for the three parameters are provided either by the user or by SPADES (see Sect ). Each parameter is determined independently of the two others, leaving these two at their initial values. At the beginning of a new iteration, the values of the three parameters are updated, therefore providing new initial values. This is repeated until the convergence condition is met: i.e. absolute differences between the parameters obtained at iteration n and n 1 all simultaneously lower than 10 K for T eff, 0.1 for log g, and for. The iterative approach was chosen over a direct approach (e.g. scanning the whole grid of 2000 reference spectra over all parameters at once), because it was computationally faster. Yet, the iterative approach could potentially be more complex in terms of convergence. The convergence performances of SPADES are assessed in Sect At any iteration, the estimation of each parameter follows the same logic: SPADES first calculates a 1D grid of synthetic spectra, varying only along the parameter that is estimated. This 1D-subgrid is calculated from the large grid of synthetic spectra described in Sect The grid is linearly interpolated to the fixed values of the other parameters. As an example, at the start of iteration n, the starting values for the atmospheric parameters derived from iteration n-1 are T eff = 5743K, log g = 4.42, and = The first parameter to be determined in iteration n is T eff. A 1D-subgrid is linearly interpolated for T eff ranging from 5200 to 6200 K by step of 200 K and log g = 4.42 and = 0.02 dex. Each synthetic spectrum is then compared (not globally, but the selected lines) to the ground-based observed spectrum according to the diagnostic appropriate for the parameters to be determined. As an example for T eff, the Hα profile of the groundbased observed spectrum will be compared to that of each synthetic spectrum. The result of the comparison is a similarity indicator for each node in the 1D-subgrid. This series of similarity indicators as a function of the parameter to be determined form a function that is analyzed (minimized or nulled, depending on the function) for the best similarity value, leading to the best estimate of the parameter to be determined. Interpolation or polynomial fit allows one to determine the parameter with a precision that is not limited by the step of the 1D-subgrid. SPADES is iteratively searching for the best matching synthetic spectrum in a 3D space (T eff, log g and ). In this first version of SPADES, the search is split in three 1D minimizations, rather than using a multi-dimensional search method. The motivation is that this facilitates visualizing and interpreting the residuals and similarity functions (See e.g. Fig. 5) that depend on a single parameter. In future versions, it would be interesting to experiment with multi-dimensional minimization algorithms such as Gauss-Newton, gradient, or Levenberg-Marquardt methods Starting parameters SPADES proceeds iteratively to derive the atmospheric parameters. It therefore needs starting values for the effective temperature, surface gravity, and iron over hydrogen ratio prior to the first iteration. Those can be deduced from photometric or bibliographic information if available. If no external information is available, SPADES offers the possibility to perform a first parameterization of the studied source, using the TGMET/ETOILE method (Katz et al. 1998; Katz 2001; Soubiran et al. 2003). The results of this first analysis are used as starting values for the iterative parameterization with SPADES. The TGMET/ETOILE method is direct (not iterative). It compares a studied spectrum to a library of reference spectra, looking for those most similar in the least-squares-sense. The parameters deduced for the studied spectrum could either be the parameters of the most similar spectrum or a combination of the parameters of the most similar spectra. Because SPADES is only looking for starting values, it simply adopts the parameters of the most similar spectrum. The performances of the SPADES implementation of the TGMET/ETOILE method (hereafter referred to as SPADES- TE)) were assessed with a Monte Carlo method method. Starting from a synthetic spectrum with T eff = 5800 K, log g = 4.5, = 1 dex (which is a node of the grid of synthetic spectra), three series of 200 studied spectra each were simulated for S/N=100, 50 and 30, respectively, and the resolving power and wavelength range of the Giraffe HR13 and HR14B setups. At S/N=50 and 100, the correct node in the grid is recovered in 100% of the cases. At S/N=30, SPADES-TE returns most of the time a solution biased by one grid step in T eff (i.e.+200 K) and two grid steps in (i.e.+0.4 dex). SPADES-TE was also tested on a Sun spectrum observed with Giraffe HR13 and HR14B setups in twilight conditions. SPADES-TE found T eff = 5800 K, log g = 3.5 and = 0.2 dex, to be compared to the canonical parameters for the Sun T eff = 5777 K, log g= 4.44, and = 0 dex. We recall that the version of SPADES presented here is limited by the GRID of synthetic spectra and by the list of reference lines to the mid-f-g dwarf stars. Future versions of the program will rely on larger grid(s) and on several fine-tuned line lists, which will extend the scope of SPADES Line and continuum identification and vertical adjustment For each spectral line, the pixels constituting the line are dynamically identified using the synthetic spectra. Using the rest wavelength of the line, SPADES starts from its center and agregates pixels in each wing as long as the value of the signal 4

5 This parameter is adjustable by the user. Figure 3 shows two examples of sets of pixels associated to a line, i.e. green crosses on the figures. Flux Flux Fig. 3. Two examples: (i) pixel identification in the reference line (green crosses), (ii) pseudo-continuum selection (blue circles), and (iii) vertical adjustment of the synthetic reference spectrum (red line) to the studied spectrum (black line). In the top plot, the pseudo-continuum area (blue circles) were selected automatically by SPADES and were used for the vertical adjustment. In the bottom plot, no valid pseudo-continuum area was found. The pixels in the line (green crosses) were used for the vertical adjustment. SPADES relies on profiles fitting. This requires adjusting the observed and synthetic spectra on the same vertical scale. It is possible, in the pre-processing, to normalize the overall continuum to 1. This is a global adjustment on the whole spectrum. This is not enough for the comparison of two lines. Therefore, SPADES systematically adjusts locally, line by line, the vertical level of the synthetic line to the one of the observed line. SPADES first identifies flat areas (pseudo-continuum) in the synthetic spectrum in the vicinity of the studied line and selects one on each side of the line. It then checks that these two areas also correspond to pseudo-continuum in the studied spectrum. If the check is positive, it uses these areas to vertically scale the synthetic spectrum on the studied one, by linear regression, i.e. there are two degrees of freedom, vertical scaling, and slope. If there are no flat areas in the synthetic spectrum or the studied spectrum contains (a) line(s) in the selected pseudocontinuum, another method is used for the vertical adjustement. The method depends on the diagnostic. For all diagnostics relying on the quadratic sum of the residuals (e.g. the iron abundance, Sect ), the pixels in the lines are directly used for the vertical adjustment (allowing here for a single degree of freedom, the vertical scale). For the diagnostic that relies on the the sum of residuals (i.e. the surface gravity, Sect ), SPADES uses pre-tabulated pseudo-continuum. An exception to this scheme is the H α line used to derive of the effective temperature, for which the pseudo-continuum and line area are pre-tabulated. Figure 3 shows the local vertical adjustement of a synthetic spectrum (red) to a studied spectrum (black). In the top plot, the pseudo-continuum area (blue circles) were selected automatically by SPADES and were used for the vertical adjustment of the synthetic reference spectrum (red line) to the studied spectrum (black line). In the bottom plot, no valid pseudo-continuum area was found. The pixels in the line (green crosses) were used for the vertical adjustment Effective temperature In this article, the stars considered are mid-f-g dwarfs, the wings of the Balmer lines are very sensitive to the effective temperature (van t Veer-Menneret & Megessier 1996; Fuhrmann 1998, 2004, 2008; Barklem et al. 2002; Cayrel et al. 2011). In this study, SPADES used the wings of H α present in Giraffe setup HR14B to derive the effective temperature. The core of H α was excluded because it is less sensitive to T eff variation and poorly modeled, specially using 1D LTE models, which is the case here. As presented in section 2.5.1, a 1D grid of synthetic spectra was calculated, differing only by the effective temperature. After pre-processing and vertical adjustment, each synthetic spectrum was compared to the studied spectrum. For each synthetic spectrum, the degree of similarity of the synthetic H α wings with the studied one, is calculated as the quadratic sum of the residuals between the pixels in the two lines: increases. If less than five pixels are associated to the line, it is considered too weak and rejected from the subsequent analysis. S im Hα = N (S (i) R(i)) 2, (4) i=0 5

6 Teff function determination Function Flux Squared residuals Teff (K) determination Function Residuals Ni abundance determination Function Flux Squared Residuals [Ni/Fe] Fig. 5. Examples, for a synthetic studied spectrum, of similarity functions as a function of effective temperature (top left), surface gravity (top right), iron over hydrogen ratio (bottom left) and nickel over iron ratio (bottom right). The polynomial (or linear interpolation for the surface gravity) used to minimize the functions are represented as blue lines. Squared Residuals line. Fig. 4. Comparison of the H α wings in the studied spectrum (black line) with the 1D grid of synthetic spectra (red lines) ranging here from 5000 (narrowest wings) to 6600 K (broadest wing). The pseudo-continuum area is labeled with blue dots and thehave been valid pixels in the H α wings are identified with green crosses. The bottom plot is a zoom on the blue wing. where N is the number of valid pixels in the H α line, S (i) (resp. R(i) the counts of the studied spectrum (resp. reference spectrum) for pixel i. Figure 4 shows the comparison of the H α wings in the studied spectrum (black line) with the 1D grid of synthetic spectra (red lines) ranging here from 5000 (narrowest wings) to 6600 K (broadest wing). The pseudo-continuum area is labeled with blue dots and the valid pixels in the H α wings are identified with green crosses. The bottom plot is a zoom on the blue wing. The degrees of similarity versus the effective temperature define a function that is minimum for the effective temperature of the studied source. The minimum is found by fitting a second-degree polynomial around the lowest similarity value. The top left plots in Figure 5 and 6 give examples of similarity functions as a function of effective temperature, for a synthetic spectrum and for the ground-based observed spectrum of ν And respectively. Each synthetic spectrum in the 1D subgrid corresponds to a red dot. The second-degree polynomial fitted to the similarity function to minimize it is represented as a blue The similarity functions for the synthetic spectra and ν And show similar degrees of regularity and smoothness. The first reason is that the spectrum ofνand has a good S/N ratio. The second reason is that in each point (red dots) the similarity functions derives from the comparison of the same noisy studied spectrum with different noiseless synthetic spectra. The noises on the points of a similarity function are correlated. The noise in the studied spectrum and the systematic differences between the studied and synthetic spectra propagate to the similarity functions, not so much as fluctuation and irregularities in the functions, but as a vertical and horizontal shift and as a flattening of the curvature of the functions Surface gravity by ionization equilibrium The ionization equilibrium method is based on the fact that the abundances determined from lines of the same element, but with different ionization stages, should be the same. Building on this idea, SPADES simultaneously models Fe lines in two ionization stages: FeI and FeII. Similarly as for T eff, a 1D grid of synthetic spectra, differing only by the surface gravity (this time) was calculated. For each a degree of ionization equilibrium (hereafter refered to as gravity similarity ) is calculated as S im grav = 1 N FeI Res l 1 N FeII Res l, (5) N FeI N l=0 FeII l=0 6

7 Teff function determination Function FeI line Squared residuals Residuals Flux Teff (K) determination Function Ni abundance determination Function FeII line Squared Residuals Squared Residuals [Ni/Fe] Fig. 6. ν And similarity functions as a function effective temperature (top left), surface gravity (top right), iron over hydrogen ratio (bottom left), and nickel over iron ratio (bottom right). The polynomial (or linear interpolation for the surface gravity) used to minimize the functions are represented as blue lines. where N FeI (resp. N FeII ) are the numbers of FeI (resp. FeII) lines used and R i the count residuals for the line l: Res l = N (S (i) R(i)), (6) i=0 where N is the number of valid pixels in the line, S (i) (resp. R(i) the count of the studied spectrum (resp. reference spectrum) for pixel i. Figure 7 shows an FeI line (top, green crosses) and an FeII line (bottom, green crosses), used by SPADES to determine of log g. The studied spectrum (log g=4.5) is denoted in black and the red lines are the synthetic reference spectra of the 1D subgrid ranging from log g=2.0 to log g=5.0 by steps of 0.5. The neutral iron line is alsmost insensitive to gravity, while the singly ionized iron line becomes stronger for decreasing gravities. The gravity similarity, as defined above, should be null for the surface gravity of the studied source. The zero of the gravity similarity function is found by interpolation. The top right plots in Figure 5 and figure 6 give an example of similarity functions as a function of surface gravity. Each synthetic spectrum in the 1D subgrid corresponds to a red dot. The linear interpolation used to find the gravity corresponding to the zero gravity similarity is represented as a blue line iron over hydrogen ratio The determination is based on a global FeI line profile fitting. For the other parameters, a 1D grid of synthetic spectra, differing only by the iron abundance, was calculated. For each Flux Fig. 7. FeI line (top, green crosses) and an FeII line (bottom, green crosses), used by SPADES to determine of log g. The studied spectrum (log g=4.5) is given in black and the red lines are the synthetic reference spectra of the 1D subgrid ranging from log g=2.0 to 5.0 by step of 0.5. synthetic spectrum, the global similarity between all synthetic neutral iron lines and all observed neutral iron lines is derived as N FeI S im Fe = Res 2 l, (7) l=0 where N FeI is the number of neutral iron lines in the list of reference lines, Res 2 l is the square of the count residuals for line l, given by formula 6. The iron content of the source corresponds to the minimum similarity degree. It is found by adjusting a third-degree polynomial through the iron similarity function. A third-degree polynomial is preferred over a second-degree polynomial, because the iron similarity function is usually asymmetric: steeper on the metal-rich side than on the metal-poor side. The bottom left plots in Figure 5 and figure 6 give an example of similarity functions as a function of the iron over hydrogen ratio. Each synthetic spectrum in the 1D subgrid corresponds to a red dot. The thirddegree polynomial fitted to the similarity function to minimize it is represented as a blue line Third step: Individual abundances Once the iterative determination of the three atmospheric parameters has converged, the individual abundances are determined. The general principle is very similar to the determination of the iron content, except that the determination is direct and not within an iterative loop with several other parameters. Each element was determined individually, using a different pre-defined list of lines. The difference with the determination method is that for each element, the 1D grid was calculated on-the-fly for the T eff, log g, and determined by SPADES and varying only by the abundance of the element to 7

8 be determined. SPADES interpolates on-the-fly the input files for the Kurucz programs (e.g. the opacity distribution functions - ODF (Castelli & Kurucz 2004)) and calls Kurucz ATLAS9 (Kurucz 2005) to calculate the atmospheric models and Kurucz SYNTHE (Kurucz 2005) to compute the synthetic spectra. The GNU-Linux ported versions of ATLAS9 and SYNTHE were used (Sbordone et al. 2004; Sbordone 2005). Element per element, for each synthetic spectrum the global similarity between all synthetic lines and all observed lines of the element X is derived as N X S im X = Res 2 l, (8) l=0 where N X is the number of lines of the element X in the list of reference lines and Res 2 l is the square of the count residuals for line l, given by formula 6. The abundance of the element X corresponds to the minimal similarity degree. It is found by adjusting a third-degree polynomial through the element X similarity function. The bottom right plots in Figure 5 and figure 6 give an example of similarity functions as a function of nickel over iron ratio. Each synthetic spectrum in the 1D subgrid corresponds to a red dot. The thirddegree polynomial fitted to the similarity function to minimize it is represented as a blue line. 3. Auxiliary data 3.1. Reference grid of models and synthetic spectra The unidimensional reference grids used to derive the stellar parameters are calculated by SPADES by interpolation using a precalculated grid. Models and synthetic spectra were calculated using the Kurucz programs ATLAS9 and SYNTHE (Kurucz 2005) ported to GNU-Linux (Sbordone et al. 2004; Sbordone 2005). About 2000 models and spectra were generated over the following range of parameters: Effective temperature: T eff from 4800 to 6800 K with a step of 200 K. Surface gravity: log g from 2.0 to 5.5 with a step of 0.5. iron abundance: from -3.0 to 1.0 dex with a step of 0.20 dex. Microturbulence:ξ=1km s 1 The models were calculated with the overshooting option switched off and l/h p = The convergence of each model was checked. The spectra were calculated for the spectral domain ranging from 600 to 680 nm, corresponding to the Giraffe/VLT HR13 and HR14 setups and a resolving power R = (SPADES convolves and re-samples the synthetic spectra to the resolving power(s) of the studied spectrum). 4. Performance assessment To assess the performances of SPADES, three series of tests were performed. The aim of the first series was to assess the internal systematic and random errors obtained as a function of the S/N ratio, but independently of the question of the initialization parameters. We also investigated the correlations between the errors on the different parameters. In the second tests, high S/N ground-based observed spectra were analyzed with SPADES to assess the external errors. SPADES proceeds iteratively to derive the atmospheric parameters. The third series of tests aims to assess the convergence properties of SPADES when it is initialized with parameters that are significantly off the true values. All tests were performed for the wavelength range and the resolving power of the HR13 (R=22 500) and HR14B (R=28 000) Giraffe setups. The three series of tests are presented below Result Teff distribution Teff (K) Result distribution Result distribution Result [Ni/Fe] distribution [Ni/Fe] Fig. 8. Distributions of the T eff (top left), log g (top right), (bottom left) and [Ni/Fe] (bottom right), estimated by SPADES for 200 S/N=50 simulated HR13 and HR14B Giraffe spectra. The true parameters are T eff = 5800 K, log g = 4.5, = 0.0 dex, and [Ni/Fe]= 0.0 dex 3.2. List of reference lines The reference lines used to derive the surface gravity, iron abundance, and individual abundances were extracted from the line list of Bensby et al. (2003). We note that the astrophysical oscillator strengths provided by Bensby et al. (2003) were not used in calculating of the synthetic spectra (see Sect. 3.1). For consistency, we used the Kurucz oscillator strengths for all lines (reference or not) Internal errors The internal systematic and random errors were assessed by Monte Carlo method. Starting from a synthetic spectrum from the reference grid with T eff = 5800 K, log g = 4.5, = 0.0 dex,ξ=1km s 1 and solar abundance scale, three series of 200 noisy spectra were generated for SNRs 30, 50 and 100. The noise was injected following a Poisson 8

9 distribution. The 600 spectra were analyzed with SPADES to estimate their T eff, log g, as well as [Ti/Fe] and [Ni/Fe] ratios as examples of abundance ratios estimates. In this first test, SPADES was initialized with the true parameters of the synthetic spectra. Table 2 presents the medians and dispersions of the distributions of residuals, i.e. estimated minus true, of the parameters derived by SPADES for each S/N ratio. Figure 8 shows the distributions of T eff, log g,, and [Ni/Fe] estimated by SPADES for the 200 S/N=50 spectra. Teff Teff The distributions of temperature and gravity residuals present no significant biases. The systematic errors on the iron over hydrogen and titanium and nickel over iron ratios are of about a few milli-dex and are most of the time one order of magnitude smaller than the random errors. The random errors are small, e.g. at S/N=30: less than 30 K for T eff, 0.15 for log g and a few hundredths of dex for the abundance ratios. we emphasize that these are internal errors accounting only for photon noise. External errors are assessed in Sect Table 2. Medians and dispersions of the distributions of residuals, i.e. estimated minus true, on T eff, log g,, [Ti/Fe] and [Ni/Fe] derived by SPADES for the S/N: 100, 50 and 30. Teff Teff Teff Teff S/N=100 S/N=50 S/N=30 med(t eff res ) (K) σ T effres (K) med(log g res ) σ log g res med( res ) (dex) σ res (dex) med([ti/fe] res ) (dex) σ [Ti/Fe] res (dex) med([ni/fe] res ) (dex) σ [Ni/Fe] res (dex) Figure 9 presents the estimate by SPADES of T eff, log g, and versus each other for the 200 stars with S/N=100 (top), S/N=50 (middle), and S/N=30 (bottom). The clearer error correlation is between the effective temperature and the iron abundance: an error of 50 K corresponds to an error of about dex. This correlation, as the others, is intrinsic to the physics of the stellar atmopsheres. While H α is weakly sensitive to metallicity (van t Veer-Menneret & Megessier 1996), the equivalent width of the iron lines is sensitive to the effective temperature External errors: ground-based observed spectra The first series of tests was devoted to internal errors. The source of noise was Poissonian photon noise. The aim of the second series was to assess the external systematic errors. There are several possible sources of external systematics. In particular, the physics of the stellar atmospheres do not reproduce the reality perfectly (e.g. discrepant oscillator strengths, local thermodynamical Equilibrium assumption, 1D modeling of the atmosphere), there are usually instrumental calibration residuals (e.g. from bias, flat-field, possibly fringing, etc.) and residuals from the observing conditions (e.g. from sky-subtraction). Fig. 9. SPADES estimate of T eff, log g, and versus each others for the 200 stars with S/N=100 (top), S/N=50 (middle), and S/N=30 (bottom) The Giraffe solar spectrum To investigate these problems, we retrieved from the Giraffe archive (Royer et al. 2005, 2008) the high S/N ratio HR13 (R=22 500) and HR14B (R=28 000) observed solar spectrum. The Giraffe solar spectrum was acquired in twilight conditions, which could alter the equivalent widths (Molaro et al. 2008). Moreover, the library of reference spectra used by SPADES was calulated with the Kurucz oscillator strengths (log gf), not with astrophysical log gf, even for the reference lines used by the different diagnostics. We therefore proceeded in two steps to analyze the Giraffe solar spectrum. First the effective temperature and surface gravity were fixed to their bibliographic values, i.e. T eff = 5777 K and log g=4.44, and only the iron over hydrogen and titanium and nickel over iron ratios were determined by SPADES. Then, in a second step, the solar spectrum was reanalyzed, this time requesting SPADES to determine the five parameters T eff, log g,, [Ti/Fe], and [Ni/Fe]. Table 3 lists the iron over hydrogen and titanium and nickel over iron ratios derived by SPADES when the effective temperature and surface gravity were fixed to their bibliographic values T eff = 5777 K and log g=4.44. The iron over hydrogen and titanium over iron ratios are underestimated by 0.08 and 0.09 dex respectively, while the nickel over iron ratio is close to the true value. The solar spectrum was then re-analyzed, requesting SPADES to also determine T eff and log g. SPADES was initialized with the bibliographic parameters Sun and converged in five iterations. Table 4 presents the parameters obtained by SPADES for the Giraffe solar spectrum. The effective temperature derived by SPADES, 5689 K, is 88 K lower than the bibliographic temperature T eff = 5777 K. Figure 10 shows a zoom on the red wing of the H α line from the Giraffe HR14B solar spectrum (black line). Two syn- 9

10 Table 3. iron over hydrogen and titanium and nickel over iron ratios derived by SPADES for a solar HR13 and HR14B solar spectrum, when the effective temperature and surface gravity were fixed to their bibliographic values T eff = 5777 K and log g=4.44. [Ti/Fe] [Ni/Fe] (dex) (dex) (dex) Table 4. Effective temperature, surface gravity, iron over hydrogen and titanium and nickel over iron ratios derived by SPADES for the Giraffe HR13 and HR14B solar spectra. T eff log g [Ti/Fe] [Ni/Fe] (K) (dex) (dex) (dex) thetic spectra are overplotted: in red a synthetic spectrum with T eff = 5689 K and in blue T eff = 5777 K (the values of the other parameters are those of the Sun: log g = 4.44 and = 0.0 dex). The red spectrum (T eff = 5689 K) agrees much better agreement with the Giraffe solar spectrum than the blue spectrum (T eff = 5777 K). Normalised flux Fig. 10. Zoom on the red wing of the H α line from the Giraffe HR14B solar spectrum (black line). Two synthetic spectra are overplotted: in red a synthetic spectrum with T eff = 5689 K and in blue T eff = 5777 K (the values of the other parameters are those of the Sun: log g = 4.44 and = 0.0 dex). The surface gravity derived by SPADES, log g = 4.46 is consistent with the Sun s gravity log g= The iron to hydrogen ratio derived by SPADES is underestimated by 0.18 dex with respect to the Sun. There are two main origins for this offset. On the one hand, there is the offset of 0.08 dex reported in table 3 that is likely due to the modeling of the spectra, in particular slightly discrepant oscillator strengths, and also to possible equivalent width alterations due to twilight observing, as reported by Molaro et al. (2008). On the other hand, as discussed in Sect. 4.1, the error on the effective temperature propagates into an error on the iron abundance. The offset of 88 K here propagates into an offset of 0.1 dex. The titanium and nickel over iron ratios derived when T eff and log g are determined by SPADES are similar to those obtained when T eff and log g were fixed to their bibliographic values. The reason is that the iron, titanium and nickel abundances are all similarly affected by the 88 K offset and therefore the effect partly canceled out in the logarithmic ratios of titanium over iron and nickel over iron S4N reference spectra The parameters of the Sun are very accurately known, which makes it a natural target to assess the performances of SPADES. On the other hand, Molaro et al. (2008) reported an alteration of the equivalent widths of the lines in the Giraffe solar spectrum, which was acquired in twilight conditions. To further test the performances of SPADES with a groundbased observed spectrum, we selected 5 stars, the Sun,νAnd (HD9826),β Vir (HD102870),µ Her (HD161797) andσdra (HD185144) in Cayrel et al. (2011). The authors derived the effective temperatures of their stars from angular diameters measured with interferometry. All their stars have diameters known to 2% or better, leading to an accuracy on the effective temperatures of about 25 K. All stars are bright and were extensively studied. Cayrel et al. (2011) adopted for log g and recent determinations from the PASTEL catalog (Soubiran et al. 2010). We proceeded similarly, adopting for log g and the median over the five most recent determinations contained in PASTEL. The titanium and nickel over iron ratios were extracted from Valenti & Fischer (2005). The spectra for the five stars were retrieved from the S4N library (Allende Prieto et al. 2004). They were degraded in resolving power and in sampling and restricted in wavelength to the Giraffe HR13 and HR14B setups. SPADES converged in two (ν And) to four (µ Her andσdra) iterations. Table 5 lists the bibliographic parameters and the parameters derived by SPADES for the five tests stars. Fig. 11 displays the temperatures, gravities, iron over hydrogen and nickel over iron ratios estimated by SPADES versus their bibliographic values. Table 6 presents the mean and dispersion of the residuals (SPADES minus bibliographic) on the estimation of the atmospheric parameters as well as the titanium and nickel over iron ratios. On average, the effective temperature is recovered with a bias of+46 K. For the same stars, Cayrel et al. (2011) found a bias of similar amplitude, but with the opposite sign, of about 90 K. In addition to the small statistics and the details of the H α profile fiting, we note that the physics used to compute the grids of models and synthetic spectra differ in several aspects. SPADES relies on Kurucz ATLAS9 models computed with the mixing length parameter l/h p = 1.25 and SYNTHE spectra (see Sect. 3.1). Cayrel et al. (2011) used the Kurucz ATLAS9 models with l/h p = 0.5 and a modified version of the Kurucz BALMER9 with the Stark broadening treatment of 10

11 Table 5. Bibliographic parameters and parameters derived by SPADES for the five tests stars. Teff SPADES (K) SPADES (dex) name T eff log g [Ti/Fe] [Ni/Fe] (K) (dex) (dex) (dex) Sun bib Sun S PADES ν And bib ν And S PADES β Vir bib β Vir S PADES µ Her bib µ Her S PADES σ Dra bib σ Dra S PADES Teff bib (K) bib (dex) SPADES [Ni/Fe] SPADES (dex) bib [Ni/Fe] bib (dex) Fig. 11. Temperatures (top left), gravities (top right), iron over hydrogen (bottom left), and nickel over iron ratios (bottom right) estimated by SPADES versus their bibliographic values. Table 6. Means and dispersions of the residuals (SPADES minus bibliographic) on the estimation of the atmospheric parameters as well as the titanium and nickel over iron ratios. T eff log g [Ti/Fe] [Ni/Fe] (K) (dex) (dex) (dex) mean σ Convergence properties First test SPADES proceeds iteratively to derive the atmospheric parameters T eff, log g, and. In the first series of tests, SPADES was initialized with the true parameters. Of course, in reality, the true parameters are unknown. The aim of this third series of tests is to assess the behavior of SPADES when it is initialized with parameters that are significantly different from the true values. For this third test, 34 single stars were generated. Their parameters were randomly drawn following uniform distributions over the following intervals: T eff in [5200, 6400] K, log g in [3.0, 4.5], and in [-2.0, 0.6] dex. No star was simulated for below -2.0, because the four Fe II lines present in the Giraffe HR13 and HR14B become too weak to allow derivating the surface gravity. This is illustrated in the top right plot of Figure 12, which shows three synthetic spectra of metallicity = 0.0 dex (blue), = 1.0 dex (red), = 2.0 dex (black). The line at Å is an Fe II line. At = -2.0 dex, the line is barely visible. Normalized flux Normalized flux FeI line Wavelength (A) Ni line Wavelength (A) Normalized flux FeII line Wavelength (A) Ti line Wavelength (A) Fig. 12. Extract of three synthetic spectra of metallicity =0.0 dex (blue), =-1.0 dex (red), =-2.0 dex (black), showing the influence of the metallicity on an Fe I line (top left), an Fe II line at Å (top right), a Ni I line (bottom left) and a Ti I line at Å (bottom right). Normalized flux Stehlé & Hutcheon (1999) and impact broadening by neutral hydrogen collisions of Allard et al. (2008). The iron over hydrogen and nickel over iron ratios are recovered without significant offsets. The titanium over iron ratio, determined from four Ti lines (without astrophysical log gf adjustement), presents an offset of 0.05 dex. To the limit of the small statistics, the error on T eff shows no trend with T eff. For each of the 34 single stars, three noisy HR13 and HR14B Giraffe spectra were simulated for S/N of 30, 50 and 100. The noise was Poissonian photon noise. The three series of 34 spectra were analyzed by SPADES. In all 102 cases SPADES was initialized with the same parameters: T eff = 5800 K, log g= 3.0 and = 1.0 dex. The left side of Figure 13 presents the differences between the true parameters of the 34 synthetic stars and the initialization parameters: T eff (top), log g (middle) and (bottom). 11