ULySS: a full spectrum fitting package
M. Koleva^{1,2}  Ph. Prugniel^{1}  A. Bouchard^{1,3}  Y. Wu^{1,4}
1  Université de Lyon, Lyon, 69000; Université Lyon 1,
Villeurbanne, 69622; Centre de Recherche Astronomique de
Lyon,
Observatoire de Lyon, St. Genis Laval, 69561; CNRS, UMR 5574;
École Normale Supérieure de Lyon, Lyon, France
2 
Department of Astronomy, St. Kliment Ohridski University of Sofia, 5 James
Bourchier Blvd., 1164 Sofia, Bulgaria
3 
Department of Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa
4 
National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, 100012, Beijing, PR China
Received 3 December 2008 / Accepted 15 March 2009
Abstract
Context. We provide an easytouse fullspectrum fitting package and explore its applications to (i) the determination of the stellar atmospheric parameters and (ii) the study of the history of stellar populations.
Aims. We developed ULySS, a package to fit spectroscopic observations against a linear combination of nonlinear model components convolved with a parametric lineofsight velocity distribution. The minimization can be either local or global, and determines all the parameters in a single fit. We use
maps, convergence maps and MonteCarlo simulations to study the degeneracies, local minima and to estimate the errors.
Methods. We show the importance of determining the shape of the continuum simultaneously to the other parameters by including a multiplicative polynomial in the model (without prior pseudocontinuum determination, or rectification of the spectrum). We also stress the usefulness of using an accurate linespread function, depending on the wavelength, so that the lineshape of the models properly matches the observation. For simple models, i.e., to measure the atmospheric parameters or the age/metallicity of a singleage stellar population, there is often a unique minimum, or when local minima exist they can be recognized unambiguously. For more complex models, MonteCarlo simulations are required to assess the validity of the solution. Methods. The ULySS package is public, simple to use and flexible. The full spectrum fitting makes optimal use of the signal.
Key words: techniques: spectroscopic  methods: data analysis  galaxies: stellar content  stars: fundamental parameters
1 Introduction
Spectroscopic data from astronomical sources contain most of the information that we can get from the remote universe, such as physical conditions, composition or motion. Often, the information extraction relies on the identification and analysis of spectral signatures due to more or less well understood physical processes. This is usually an interactive task and requires specialised expertise, although this is becoming less true with time (e.g., Sousa et al. 2007).The recent data avalanche prompted a search for automated, objective (and generally more efficient) methods. One possibility is to directly compare an observed spectrum with a set of models or an empirical library of objects with known characteristics on a pixel by pixel basis. This can be done for the analysis of stellar atmospheres, as in BailerJones et al. (1997), Katz et al. (1998, TGMET), Prugniel & Soubiran (2001), Snider et al. (2001), Willemsen et al. (2005), Shkedy et al. (2007) and RecioBlanco et al. (2006, MATISSE). It has also been used for the determination of the history of stellar populations, as in Heavens et al. (2000, MOPED), Panter et al. (2003), Cid Fernandes et al. (2005, Starlight), Moultaka (2005), Ocvirk et al. (2006a,b) and Tojeiro et al. (2007, VESPA).
The main advantage of these methods is that, rather than picking specific features, they make optimal usage of the entire measured signal. This comes with a price, however, as we forfeit our ability to draw direct relations between the strength of a specific feature and a physical characteristic. This lack of simplicity should ultimately be regarded as inherent to the nature of the spectra, where the information is redundant and distributed (possibly uniformly) over a large wavelength range.
This paper presents ULySS (Université de Lyon Spectroscopic analysis Software)^{}, a new software package enabling full spectral fitting for two astrophysical contexts: The determination of (i) stellar atmospheric parameters and (ii) star formation and metal enrichment history of galaxies. Many similarities between these two cases allowed us to build a single package capable of handling both.
In ULySS, an observed spectrum is fitted against a model expressed as a linear combination of nonlinear components, optionally convolved with a lineofsight velocity distribution (LOSVD) and multiplied by a polynomial function. A component is a nonlinear function of some physical parameters (e.g., , and [Fe/H]). The multiplicative polynomial is meant to absorb errors in the flux calibration, Galactic extinction or any other cause affecting the shape of the spectrum. It replaces the prior rectification or normalization to the pseudocontinuum that other methods require. This model is compared with the data through a nonlinear leastsquare minimization.
ULySS has been used to study stellar populations (Koleva et al. 2008c; Michielsen et al. 2007; Koleva et al. 2009,2008a; Koleva 2009) and determine atmospheric parameters of stars (Prugniel et al. 2009).
The goal of this paper is to describe the package and to explore the potential and caveat of direct comparison between an observation and a composite model. A previous implementation of the same idea, briefly described earlier and named NBURSTS (Chilingarian et al. 2007; Koleva et al. 2007), was a variant of the pPXF program (Cappellari & Emsellem 2004, hereafter CE04) applied to stellar populations. Another variant of pPXF was developed by Sarzi et al. (2006) to fit at the same time the stellar and the gas kinematics. ULySS is widening the scope of the package to the measurement of atmospheric parameters and potentially other applications. In addition, the algorithm has been optimized in its mathematical details to improve its precision, robustness and performance. We also wrote documentation and tutorials and made the package userfriendly to facilitate its public distribution.
The paper is organized as follow: in Sect. 2 we describe the method and the package. In Sect. 3 we illustrate the approach with the case of the determination of the atmospheric parameters of a star. In Sect. 4, we give examples of the analysis of a stellar population. Section 5 gives conclusions.
2 Description of the method and package
The method consists of minimizing the between an observed spectrum and a model. The model is generated at the same resolution and sampling as the observation and the fit is performed in the pixel space. This optimizes the usage of the information and simplifies the rigorous treatment of the error on each spectral bin.
The method proceeds in a single fit to determine all the free parameters in order to handle properly the degeneracies between them (e.g., the temperaturemetallicity for a stellar spectrum). Other methods estimate the parameters in different steps. For example, one may measure the temperature using criteria almost insensitive to the metallicity and in turn, using this temperature, derive the metallicity. If in fact the criteria to obtain the temperature is metallicitydependent, the absolute minimum will not be reached unless the fit is iterated. Using a single minimization does not require this additional complexity.
Ocvirk et al. (2006a,b) presented another implementation of this approach, using a nonparametric regularized fit. It performs very well in reconstructing the evolutionary history of a stellar population. But one of its limitations is the difficulty to grant a degree of confidence for the solution (and it is tailored to stellar populations only). This is the main reason which lead us to work on a ``simpler'' parametric minimization, allowing us to better understand the degeneracies and the structure of the parameter space by constructing maps.
The parametric local minimization presented here has been introduced and compared with STECKMAP in Koleva et al. (2008c) in the case of single stellar populations.
2.1 Parametric model
The observed spectrum,
,
is approximated by a linear
combination of k non linear components, CMP, each with weight W.
This composite model is possibly convolved with a LOSVD and multiplied
by an nth order polynomial,
,
and summed with another
polynomial,
:
The LOSVD is a function of systemic velocity, , velocity dispersion and may include GaussHermit expansion (h3 and h4, van der Marel & Franx 1993). is the logarithm of the wavelength (the logarithmic scale is required to express the effect of the LOSVD as a convolution). The CMP_{i} must be tailored to each problem. For instance, to study a stellar atmosphere, the CMP will be a function of the effective temperature, , surface gravity, g, and metallicity, [Fe/H]. Or, for a stellar population, it will be a function of age, [Fe/H], and [Mg/Fe], returning the spectrum of a single stellar population (SSP). In general, a CMP is a function of an arbitrary number of physical parameters and of .
The importance of the multiplicative polynomial has been discussed in Koleva et al. (2008c), as it absorbs the effects of an imprecise flux calibration (a common issue in smallaperture spectroscopy) and of the Galactic extinction (which could also be explicitly included in the model). The effect of is studied in Sect. 3.2 in the particular case of the determination of atmospheric parameters of a star. A similar study was made in Koleva et al. (2008c) in the case of stellar populations. In all the practical cases that we studied, was not degenerate with the parameters of the CMP, and high orders could be used (though is often sufficient to obtain an unbiased estimate of these parameters). The optimal value of n, which depends mostly on the resolution and wavelength range, can be chosen with ULySS.
The additive polynomial is certainly more subtle to use, and is, in most cases, unnecessary. It is indeed degenerate with the intensity, depth or equivalent width of the lines and may bias determinations of the CMP parameters ( , age or [Fe/H], depending on the context). Such a term may be included to account for dataprocessing errors (under or oversubtraction of a smooth background), or may have a physical origin. When determining the stellar kinematics by fitting one or several constant CMPs (i.e., fixed spectra, like stars or models without free parameters), the additive term may be required to absorb the template mismatch (Koleva et al. 2008b). Omitting the additive term would in that case bias the measurement of the velocity dispersion. For this reason, an additive term is generally present in the packages used to determine the internal kinematics of a stellar population. For instance, in CE04 it is an explicit polynomial, or in Fourier quotient programs (Sargent et al. 1977), the spectrum is continuum or meansubtracted before the fit and the amplitude or the features is a free parameter (this is equivalent to an additive term). Another case needing an additive term is the analysis of a stellar population in the presence of the nebular continuum emission of an AGN (in that case a powerlaw is more appropriate than a Legendre polynomial; see the tutorial distributed with the package). The biases on the stellar population parameters may be determined from simulations (Prugniel et al. 2005; Prugniel & Koleva 2007).
In summary, we can only advise careful usage of the additive terms. They are only needed in rare circumstances and including them without valid reasons (or failing to do so when they are required) may severely bias the results of the analysis. The lowering of the that is likely to result from the inclusion of additional degrees of freedom should not be the only criterion to evaluate the validity of the model. The user should check the effects of these terms on the CMP parameters using simulations.
2.2 Algorithm
As written in Eq. (1), the problem is a fit to a multilinear combination of nonlinear functions. The parameters of each CMP_{i} are in general nonlinear and are evaluated together with those of the LOSVD with a LevenbergMarquardt (Moré et al. 1980; Marquart 1963) routine (hereafter LM). To measure the h3 and h4 coefficients of the LOSVD, the package implements the penalization proposed by CE04 to bias them to 0. This option, not illustrated here, is useful to analyse low S/N data.
The (linear) coefficients of the and polynomials are determined by ordinary leastsquares at each evaluation of the function minimized by the LM routine. The weights of each component, W_{i}, are also determined at each LM iteration using a bounding value leastsquare method (Lawson & Hanson 1995) in order to take into account constraints on the contribution of each component, like forcing positivity.
An alternative would have been to use the variable projection method (Golub & Pereyra 1973) to solve this separable nonlinear leastsquares problem (i.e. where the model is a linear combination of nonlinear functions). We did not use this solution because the bounding of the linear parameters is not possible in the public implementation of this method (Bolstad 1977, the VARPRO program); in addition there was no version of this algorithm in the IDL/GDL language used for this project. Therefore we preferred to adjust the linear parameters in a separate layer at each LM iteration, as in Cappellari & Emsellem (2004). We stress that separating the linear variables is important, not only for performance but also for stability.
The optional bounding on the W_{i} has to be adapted to each particular problem. For example, when decomposing a stellar population in a series of bursts, each component must have a positive or null weight. The LM implementation by Moré et al. (1980, in MINPACK^{}# allows the user to define limits for the values of any of the parameters of the .
For and , we use Legendre polynomial developments of respectively orders n and m, for their orthogonality properties. However, the developments intervening in our problem are not strictly orthogonal because (1) the errors are not uniform along the axis and (2) there may be gaps in the signal corresponding to masked pixels (missing or damaged values). For these reasons, the and are determined by an ordinary least square fit. Defining ad hoc orthogonal polynomials would have been equivalent both from a mathematical and a performance point of view.
Note that as the LM minimization approaches to the solution using the local gradients of the CMPs, there is no need or advantage to apply a linear transformation to these functions. For example, applying a rotation in the age  metallicity plane in order to minimize on two orthogonal parameters does not ease the convergence or avoid local minima: The path to the solution is not affected. But a proper choice of the parameters may lead to a better convergence. For instance, to fit a stellar spectrum, minimizing on log( ) or on is preferable to directly using .
2.3 Line spread function
To compare a model to an observation, both must have the same spectral resolution, or we must first transform either the model or the observation in order to match their respective resolution.
The spectral resolution is characterized by the line spread function, LSF, which is the analog in the spectral direction of the PSF (point spread function) for images. This is the impulse response describing the wavelength distribution of the flux of an unresolved spectral line. The LSF results from the convolution between the intrinsic resolution of the spectrograph and the slit function. In a first approximation, it is represented by , where l is the wavelength (linear scale) and is the FWHM of the LSF.
In practice, the LSF may not be defined by a single number. It is not necessarily a Gaussian, and may change with wavelength and position in the field for integralfield or longslit spectroscopy. Usually, the model has a higher resolution than the observation. Otherwise, the analysis will not make optimal usage of the available information (i.e., the high resolution details in the observed spectrum will not be exploited).
To make the model comparable to the observations, we proceed in two steps. First we have to determine the LSF and then inject it in the model.
2.3.1 Determination of the linespread function
Figure 1: Effect of using a precise LSF, illustrated with a fit of a Vazdekis/Miles spectrum with a Pegase.HR/Elodie.3.1 SSP component. The top panel shows the spectrum (in black) and the best fit in blue (both are almost superimposed and the black line can be seen only when zooming on the figure), the red line is the multiplicative polynomial. The yellow regions were rejected from the fit (rejection of flagged telluric lines and automatic rejection of outliers). The middle and bottom panels are respectively the residuals from the best fit when (i) assuming a constant Gaussian LSF (in ) or (ii) a matched LSF. The continuous green lines mark the 1 deviation, and the dashed line is the zero axis. The right side of the figure expands a small wavelength region, around Ca4227. 

Open with DEXTER 
The function that we are seeking is the relative LSF between the model (which has a finite resolution) and the observation. It should normally be determined using calibration observations.
Three types of calibrations can be considered:
 1.
 Arc lamp spectra. They are routinely produced during the observations and are used to adjust the dispersion relation and to achieve wavelength calibration. The slit of the spectrograph is uniformly illuminated with a discharge lamp (for example HeAr) producing narrow emission lines. The position of chosen unblended lines are used to fit the dispersion, and the width and shape can be used to determine the LSF.
 2.
 Standard star. Normally any star, except some hot stars with featureless spectra used for the spectrophotometric calibration, can be used to determine the LSF. The observed spectrum may be compared with a highresolution spectrum of the same star, or with a model of this star, to determine the broadening due the the spectrograph. ULySS can be used to measure this broadening. (Beware that sometimes spectrophotometric standards are observed with a widened slit, and are not usable for LSF calibration.)
 3.
 Twilight spectrum. Spectra of the twilight sky are often used to determine the variation of the sensitivity over the field of the spectrograph. These spectra result from the uniform illumination of the slit by a Solar spectrum and can therefore by used as any stellar spectrum to measure the broadening.
The two other solutions use stellar spectra. As the physical models presented in this article are based on empirical libraries of stellar spectra, a proper choice of the calibration stars can make the task of determining the LSF straightforward (the solar spectrum is included in most libraries). If the exact star is not available at the model's resolution, a similar star, or an interpolated spectrum, may be used. Using a stellar spectrum bypasses some of the difficulties met with arc spectra and can directly gives the relative LSF.
There is still an important phenomenon to consider when describing the LSF. Often, the intrinsic resolution of the spectrograph is significantly higher than the actual resolution, which is limited by the slit width. As a consequence, the distribution of light within the slit has an effect on the spectrum. In particular, the apparent broadening of a star observed under excellent seeing conditions (seeing smaller than the slit) will be smaller than the broadening observed for an extended object (or a star with poor seeing conditions). The effective resolution results from the product of the light profile of the object and the slit function, convolved by the intrinsic resolution of the spectrograph.
This effective resolution depends on the observing conditions and on the light profile of the source. It may vary between consecutive exposures. This problem may be difficult to correct, and by limiting the knowledge of the instrumental broadening, it hampers the measurement of the physical velocity dispersion. In the types of analysis discussed in this paper, this effect preserves the determination of the other parameters (metallicity, age or ).
Hints to this effect of resolution of an object within the slit may be obtained when comparing the LSF derived with various standard stars and those obtained with twilight spectra.
A practical means to measure the LSF is to determine the broadening function (cz, and possibly h3 and h4) in a succession of small wavelength intervals. For spectra with R = 1000 to 3000, we typically use segments of 200 Å, separated with 100 Å steps (they overlap by half their length).
2.3.2 Injection of the LSF in the model
Because the LSF varies with wavelength, it cannot be injected in the model as a simple convolution. The method we use consists of convolving the models by the series of LSFs determined at some wavelength and then interpolating linearly in wavelength between the convolved models.
To illustrate the importance of matching the LSF in the spectrum fitting process, we show in Fig. 1 the analysis of a spectrum with a Pegase.HR population model based on the Elodie.3.1 library. The analysed spectrum is a stellar population model from the library of SSPs computed by Vazdekis^{} with the Miles library (SánchezBlázquez et al. 2006). The first fit simply assumes a purely Gaussian and constant LSF. The second uses an optimal LSF. The residuals of the first fit are minimal in the center, but become larger at the edges; a zoom in a small region shows that this is due to a misfit (not to noise). Using the proper LSF corrects this defect and the residuals become smoother. In this particular example, despite the considerable effect on the residuals, the incidence on the stellar population parameters is marginal. The precision of the LSF mostly affects the parameters of the LOSVD. By suppressing the LSF mismatch, we can search for other signatures in the residuals which could have been masked otherwise.
The LSF injection also corrects possible inaccuracies in the wavelength calibration. An example of this is shown in Koleva et al. (2008c), where a wavelength calibration systematic distortion affecting the Bruzual & Charlot (2003) models is corrected.
2.3.3 Example: the SDSS LSF
As an example of LSF analysis, we use the velocity dispersion template stars from the SDSS copied from http://www.sdss.org/dr5/algorithms/veldisp.html. These 32 G and K giant stars from M 67 were used to determine the velocity dispersion of the galaxies as an average between estimates obtained by Fourierfitting and directfitting methods.In Fig. 2 we show the LSF relative to the Elodie.3.1 library obtained with ULySS (Elodie.3.1 is restricted to the wavelength range 39006800 Å). It was determined using wavelength intervals of 600 Å spaced by 300 Å. The variation of the instrumental velocity dispersion ( ) with wavelength is significant: from 50 to 75 km s^{1}.
The classical methods for measuring the (physical) velocity dispersion, , (Rix & White 1992; van der Marel & Franx 1993; Bender 1990; Tonry & Davis 1979; Franx et al. 1989; Sargent et al. 1977; Cappellari & Emsellem 2004) compares the observation to stellar templates observed with the same setup. As, in general, the LSF changes are moderate, the redshift of the galaxy does not significantly affect for nearby galaxies. But for a distant galaxy, neglecting the shift of the LSF may give a measurable effect. For a lowmass galaxy with km s^{1} measured with the SDSS, the bias would be about 0.1 km s^{1} at z=0.03but 2 km s^{1} at z=0.4.
2.4 Description of the ULySS package
ULySS, available at http://ulyss.univlyon1.fr/, has been programmed in the IDL/GDL language. This is the language of the widely used proprietary IDL (Interactive Data Language)^{} software. The open source GDL (Gnu Data Language)^{} interpretor can also be used to run ULySS. The choice of this programming language is essentially historical: Many of the required routines were already available as public libraries and development was started from the pPXF package (Cappellari & Emsellem 2004). We may offer a version written in a modern language (most likely Python) in the future.
Figure 2: Line spead function of the SDSS. 

Open with DEXTER 
ULySS contains various programs and subroutines that can be used to:
 define the array of components to fit;
 perform the fit;
 visualise the results;
 make maps, convergence maps and MonteCarlo simulations;
 read data from FITS files;
 test the package.
The package contains routines to define various types of CMPs, notably to analyse a stellar atmosphere (TGM) and to analyse a stellar population (SSP), and the construction of other CMPs by the user was made as easy as possible. To fit a composite model, i.e., a linear combination of components, one can simply concatenate several CMPs in a single array.
The core of the package is the local minimization described in Sect. 2. Such a minimization starts from a point (guess) in the parameter space, whose choice may be critical. To release this constraint, ULySS makes it easy to perform a global minimization by providing vectors instead of scalars as guesses.
The most valuable aspect of the package is the possibility of exploring and visualizing the parameter space. The tools offered for this purpose are (i) MonteCarlo simulations (ii) convergence maps, and (iii) maps.
MonteCarlo simulations are performed to estimate the biases, the errors and the coupling (degeneracies) between the parameters. A simple option in the main fitting program allows us to perform a series of minimizations with random noise added to the data. This noise has the same characteristics as the one estimated in the data. Normally the noise should be carried throughout the datareduction process, starting from the statistical noise of the detector, that can usually be securely estimated. During the data reduction, the signal is likely to be resampled, when the spectra are extracted from the initial 2D frame, and when they are calibrated in wavelength (or logarithm of wavelength). This operation introduces a correlation between the pixels (see de Bruyne et al. 2003), which is represented in ULySS by the ratio between the actual number of pixels and the number of independent pixels. Using it, the MonteCarlo simulations reproduce the correct noise spectrum and gives a robust estimates of the errors.
Convergence maps are tools to evaluate the convergence region, i.e., the domain of the parameter space from which guesses converge to the absolute minimum of the . These maps can be generated using the global minimization approach explained above.
maps are visualizations of the parameter space. They are generated by (i) choosing a 2D projection (e.g., age and metallicity) and a node grid in this plan and (ii) performing an optimization over all the other axes of the parameter space for each node of this grid. These maps allow the identification of degeneracies and local minima.
Typically, when approaching a new problem, like using a new CMP, or a new wavelength range or region of the parameters space, these three tools can be used to understand the reliability of the results before proceeding to the analysis of a massive dataset. Their usefulness is presented in the next sections.
3 Determination of stellar atmospheric parameters
The effective temperature, , surface gravity, , and metallicity, [Fe/H], are fundamental characteristics that can be derived from spectroscopic analysis (Cayrel de Strobel et al. 2001). Full spectrum fitting, as provided by ULySS, can be used for this purpose: the program will identify the best matching , and [Fe/H] by comparing to a model.
ULySS carries out a parametric minimization. So, the core of the problem is to obtain a parametric model, i.e. a function returning a spectrum given a set of atmospheric parameters. The reference spectra are either a grid of theoretical models (e.g. Coelho et al. 2005; Munari et al. 2005) or a set of observed stars whose parameters are known from the analysis of individual high resolution spectra (e.g. Soubiran et al. 1998). ULySS requires an interpolator of this grid. In the present paper, we use the one presented in Prugniel & Soubiran (2001) for the ELODIE library. In brief, it consists of polynomial approximations of the library. Three overlapping ranges of are considered (hot, warm and cold) and linearly interpolated to produce the final function. In each range each pixel of the spectrum is computed as a 21 term polynomial in , and [Fe/H]. The coefficients of these polynomials were fitted over the 2000 spectra of the library. The choice of the terms, of the limits and of weights were finetuned to minimize the residuals between the observations and the interpolated spectra. In this paper, we use the interpolator built on the continuum normalized spectra. An alternative solution to this global polynomial representation of the library would have been to use a local approximation based on a gaussiankernel smoothing, as in Vazdekis et al. (2003, Appendix B).
ULySS defines a model component (CMP in Eq. (1)) for this model. The TGM component, as we named it, allows to perform the minimization on the three atmospheric parameters. With the current version of the ELODIE library (Prugniel et al. 2007b, version 3.1), the temperature range is limited to < < 30 000 K. In future versions (Wu 2009, in preparation), a greater number of hot ( > 10 000 K) and cold ( < 4200 K) stars will be included, extending the current validity range of the interpolator.
3.1 TGM fit example
We analysed the 18 stars from the CFLIB (indoUS, Valdes et al. 2004) library of spectra in common with the study by da Silva et al. (2006) of G & K stars using spectroscopy and LTE models. We performed a fit with a TGM component starting from a grid of guesses in order to be independent of the prior knowledge of the parameters. The LSF was determined by using several stars in common between CFLIB and the ELODIE library. The results, presented in Fig. 3, are consistent with those of da Silva et al. (2006), except for HD 189319 where we find a significant discrepancy in metallicity. The measurements from da Silva et al. (2006) are: = 3978 K = 1.10 g cm^{2}, [Fe/H] = 0.29 and from ULySS: = 3904, = 1.77, [Fe/H] = 0.10. It is the most discrepant star for both metallicity and gravity; it is also the lowest gravity and coolest star of this sample. Another recent spectroscopic analysis gives: = 4150 = 1.70 [Fe/H] = 0.41 (Hekker & Meléndez 2007); interferometric measurements tend to indicate a lower : 36503800 K, and hence probably lower gravity = 0.9 (Wittkowski et al. 2006; Neilson & Lester 2008). This discrepancy is not inconsistent with ULySS, but the ELODIE interpolator surely deserves to be improved in this region of the HR diagram.
Figure 3: Comparison of the atmospheric parameters determined by ULySS with those from high resolution spectra (da Silva et al. 2006). The abscissas are the measurements from da Silva et al. (2006). On the top panels, the ordinates are from ULySS and the green lines are the diagonals. On the bottom panels the ordinates are the differences ULySS literature. 

Open with DEXTER 
For this sample, we found (i.e., K), and . Excluding the discrepant M star we obtain: . The deviations reported here are the standard deviations on individual measurements. The temperatures found by ULySS are systematically cooler by 60 K than those of da Silva, consistent with the offset mentioned by these authors in their own comparison to the literature. They found a systematic difference of 39 to 50 K, their measurements being hotter. The three atmospheric parameter error estimates from our program are similar to those given by MC simulations and are about 20 times smaller than the ``external'' errors, so we did not draw error bars in Fig. 3, nevertheless the deviations (external errors) are identical to those reported by da Silva et al. (2006). We can conclude that the measurements performed with ULySS are precise and reliable.
Table 1: Stellar atmospheric parameters for six CFLIB stars of different spectral types.
In Table 1 we selected six CFLIB stars representative of the various spectral types and luminosity classes. Figure 4 presents the fit for a Solar type star from this list. A detailed discussion of the CFLIB stellar library will be made in a separate work.
Figure 4: Fit of the CFLIB star HD 76151 with a TGM component. The symbols and conventions are the same as in Fig. 1. The order of the multiplicative polynomial is n=200. The right side expands a small wavelength region around Mg_{b}. 

Open with DEXTER 
3.2 Multiplicative polynomial continuum
Most stellar analysis programs first require the observed spectrum to be normalized to a pseudocontinuum, which can be determined either interactively or automatically. By contrast, ULySS determines this normalization in the same fitting process by including a multiplicative polynomial, in Eq. (1), in the model. This singlestep fitting procedure insures that the minimum can be reached and allows one to check the possible dependences between this continuum and the measured physical parameters.
Figure 5 presents the results of a series of fits of the six representative stars from Table 1, varying the order of . The observations consist of 8300 independent pixels in the wavelength range 39006800 Å. As there is no external estimate of the noise, we gave a constant weight to all the pixels (except those rejected as outliers), and computed the noise in order to reach = 1 for n=200. We explored the multiplicative polynomial order range 0 < n < 800; while n is increasing, the value of the decreases as a power law.
Figure 5: The evolution of the stellar atmospheric parameter fit results (, log( ), , and [Fe/H]) with increasing Legendre polynomial degree for 6 example CFLIB stellar spectra. Black, green, red, blue, violet and gold colors are for each star listed in Table 1. 

Open with DEXTER 
The atmospheric parameters converge rapidly to their asymptotic values (defined here as the mean of the solutions for n > 25). For the F, G, K and O stars the plateau is reached between n = 10 to 15 (the stability of the solution is lower for the O star, in particular for metallicity). For the M star, the plateau is reached at n = 35, but the fit is not stable above n=150. The A1Ib spectral type CFLIB star HD 195324 displays a significant dependence between n and the measured metallicity; it did not stabilize to a plateau. This is likely due to the limited quality of the ELODIE V3.1 interpolator in this underpopulated region of the parameter space. The ELODIE library counts only five Atype supergiants and therefore the interpolation is not secure. In Atype stars the ELODIE continuum is taken in the flanks of the Balmer lines and the analysis relies on the multiplicative polynomial to fit them. Using a fluxcalibrated model improves the situation.
Note that the variations of the parameters with nare slightly larger than the error bars. On , and [Fe/H] the errors are typically 0.1%, 0.006 and 0.005 while the dispersion of the values for n>20 are 0.2%, 0.01 and 0.01, i.e. about twice the errors. The extremely small internal errors hide some potential biases of either observational or physical origin.
Further evidence for the nondegeneracy between the atmospheric parameters and n, even for values of nconsiderably larger than what is used in practical cases, is given by the MonteCarlo simulations of the next section. In the case of degeneracy, the error bars computed by the fitting program would be underestimated.
In order to check if the high values of n are overfitting the data (i.e. fit the noise), we carried out similar experiments with noise spectra having the same characteristics as the data. The measured decreases as expected, to reach 0.99 for n = 100, 0.98 for n = 400 and 0.96 for n = 800. It is clear that the trend seen for the observation is not due to the overfitting, as the slope is much larger. The decrease of is probably due to the shape of the continuum being increasingly better fitted when n rises, and to some extent the physical effects ignored by this simple model. The effect of rotation for hot stars may be the main influence, moreover detailed abundances cannot be mimicked by the multiplicative polynomial.
The importance of this multiplicative polynomial for stellar population studies is discussed in Koleva et al. (2008c). Within the same wavelength range, the fits reach a plateau for lower n, probably because the models are flux calibrated.
3.3 Convergence, maps and MonteCarlo simulations
ULySS also includes the tools to assess the significance and validity of the results. Figures 6 and 7 illustrate the usage of convergence maps, MonteCarlo simulations and maps to explore the parameter space.
The convergence map, Fig. 6, shows two basins. In the wide region defined by < 10 000 K, any choice of initial guess will converge toward the correct solution, while hotter guess may fall into a local minimum in an unphysical region.
Figure 6: Convergence maps on different projections of the parameters space for the CFLIB star HD 76151 inverted with the TGM component. Red crosses stand for the global minimum solutions found by ULySS. 

Open with DEXTER 
MonteCarlo simulations are used to estimates the errors when the different parameters are not independent. In Fig. 7 the errors determined by MonteCarlo simulations are compared with those given by the minimization procedure. Though both are in approximate agreement, only the simulations can show the effect of the coupling of the errors between the parameters. The simulations consist of series of analyses of a spectrum plus a random noise corresponding to the estimated noise. The added noise has a Gaussian distribution and takes into account the correlation between the pixels introduced along the processing, as stressed by de Bruyne et al. (2003). This latter effect is modeled by keeping track of the number of independent pixels during the steps of the processing, and then generating a random vector of independent points that is finally resampled to the actual length of the spectrum.
Figure 7: MonteCarlo simulations and maps for the CFLIB star HD 76151 inverted with the ELODIE library as presented in Fig. 4. The three projections of the parameters space are presented. The 1000 MonteCarlo simulations are performed adding a random noise equivalent to the estimated one. The superimposed crosses give the internal errors estimated by ULySS, and the ellipses the standard deviation computed from the simulations. 

Open with DEXTER 
The maps complete the MonteCarlo simulations by revealing the degeneracies and the presence of local minima. We built the map by choosing a grid of nodes in a 2 dimensional projection of the parameter space, and performing the minimization for each node (hence optimizing n2 parameters). Any local minimum can be detected, providing that the grid is fine enough. The topology of these maps also indicates the degeneracies. In Fig. 7, showing the measurement of the three atmospheric parameters for a Solar type star, the maps are regular, with weak dependences between the parameters and a single minimum. When using more complex models, like a composite stellar population, the maps often show local minima.
4 History of stellar populations
By using the SSP component (single stellar population) provided with ULySS, one can evaluate many evolutionary parameters from an integrated spectra. As in the case of the TGM component, the SSP component describes the fitting boundaries and the overall recipe to create SSP spectra and fit them to the observed data. This time, the CMP is characterised by age, [Fe/H] and [Mg/Fe] (this last dimension is currently only available in semiempirical models under development, see Prugniel et al. 2007a). A grid of SSPs given as input is splineinterpolated to provide a continuous function.
The CMPs can be linked to a number of population synthesis models. Koleva et al. (2008c) tested 3 of them: Galaxev (Bruzual & Charlot 2003), Pegase.HR (Le Borgne et al. 2004) and VazdekisMiles (SánchezBlázquez et al. 2006; Vazdekis 1999). They concluded that Pegase.HR and VazdekisMiles are reliable and consistent.
The first step towards reconstructing the star formation history (SFH) of an object is to calculate its SSPequivalent parameters by using a single CMP that interpolates a grid of SSP in age, [Fe/H] and possibly [Mg/Fe]. These SSPequivalent parameters to some approximation correspond to the luminosityweighted average over the distribution in age and metallicity (but see Trager & Somerville 2009, for a discussion of this simplification). The present method has been used by Koleva et al. (2008c) who have shown that reliable information can be retrieved. The metallicity of Galactic globular clusters can be compared to the determinations derived from spectroscopy of individual stars with a precision of 0.1 dex, which is the actual precision of these latter measurements.
If the object has a complex SFH, with several epochs of star formation, SSPequivalent parameters are essentially representative of the star formation burst that dominates the light (often the most recent). The ULySS package can be used to reconstruct a detailed SFH, generally limited to 2 to 4 epochs of star formation because of the degeneracies and the finite quality of the models and observation (see also Tojeiro et al. 2007).
4.1 Complex stellar population: application to NGC 205
The galaxies have in general a complex SFH and retracing the star formation rate along the history is a fundamental piece of information to understand the physics of the galaxies and for the cosmology. In principle, one can access such information by directly fitting a positive linear combination of many SSPs, but such approach would be unstable because of the degeneracies between the components and would require a regularization.To circumvent these degeneracies, we start with simple physical assumptions, like the presence of an old stellar population, then divide the time axis in intervals (by setting limits in two or more intervals). As the number of free parameters increases, local minima appear and a global minimization is required; the maps help to understand the structure of the parameter space. Usually the fit is performed several times, with an increasing number of components and varying limits on the population boxes. Then, doing MonteCarlo simulations and checking the residuals of the fits, it is possible to assess the relevance of the solutions and select the most probable SFH.
As an example, we analyse the star formation in the inner 2 (roughly the size of the nucleus, Butler & MartínezDelgado 2005) of NGC 205. The data, taken from Simien & Prugniel (2002), have in the central region and a spectral resolution of . For this present demonstration, we analyse this spectrum in terms of two epochs of star formation (i.e., 2 CMPs, each one being an SSP): one ``young'' (age <800 Myr) and one ``old'' ( Myr). This hypothesis is not bound to an a priori knowledge of the stellar population; it is essentially a choice of time resolution. Depending on S/N, resolution and wavelength range, the number of components may be increased; e.g., in Koleva (2009) the same data are decomposed in four epochs of star formation.
Figure 8: Fit result for the central 2 of NGC 205 with 2 SSP components. In the top panel the black line represents the observed spectra, the blue line is the fitted model and the red line shows the multiplicative polynomial ( ). The bottom panel shows the residuals ( black), mean ( green dashed) and 1 deviation ( solid green). In both panels, some of the data were not considered for the fit ( yellow). 

Open with DEXTER 
Table 2: Ages and metallicities of the central 2'' of NGC 205.
Figure 9: MonteCarlo simulation ( dots) and map ( contours) for a 2 CMP fit to spectra of the inner 2 of NGC 205. The blue contours represent the levels for the young stellar population (age constrained to be <800 Myr) and the red contours are for the old population (age between 800 and 14 000 Myr). 

Open with DEXTER 
Figure 8 and Table 2 shows the results of our analysis. For the young component we find an age of 130 Myr, consistent with photometric results (J, H, K photometry) from Davidge (2003) and with Cepa & Beckman (1988) who found from orbital considerations that NGC 205 crossed the disk of M 31 about 100 Myr ago. It represents about 25% of the light and only about 7% of the mass.
The MonteCarlo simulations of Fig. 9 show that we can distinguish between the two populations, in the sense that the two clouds corresponding to these populations are well separated. However, the existence of two bursts was one of our hypotheses and to test its validity we will apply the same analysis to the best fit SSP of the first NGC 205 experiment. We find that a young burst would be detected in about 10% of the cases in MonteCarlo simulations, but is easily rejected as the solution does not cluster around a marked minimum.
The mean solutions estimated from the MonteCarlo simulations with two bursts, given in Table 2 (lines noted ``MC'' in the column), are compatible with the direct fit, but the errors are significantly larger (because of the degeneracies). It is also interesting to see that the analysis with a 2burst model of the bestfit SSP produces an unbiased solution (lines noted ``MC1'').
Examining Fig. 9 in more detail, we note that the distribution of the MonteCarlo solutions are not simple Gaussians centered on the direct solution. For the old burst, there is a small and concentrated secondary cloud with 10% of the solutions at an age of about 3 Gyr and [Fe/H] = 1. The solutions belonging to that cloud also form a tail at the old side of the young burst. This feature is associated with a local minimum detected on the map whose depth is almost similar to the solution. Because of the random noise, this local minimum can become the actual solution. This means that the two minima cannot be statistically distinguished. The morphology of the map gives the impression that the oldest cloud is an artifact due to a discontinuity in the model. But there is no objective argument to reject this (less probable) solution.
4.2 Exponentially declining or constant SFR populations
Our last series of experiments is the analysis of simulated populations with an extended star formation history. We consider two scenarios: a constant star formation rate (SFR) stopped 500 Myr ago and an exponentially declining SFR with a characteristic time, , of 5 Gyr. These spectra were computed with Pegase.HR and ELODIE.3.1. In both cases we analysed the simulated galaxies with a model made of 3 CMPs defined by the age limits: [12, 800], [800, 5000], [12 000] Myr. The goal is to see if we can distinguish between different scenarios of star formation.
Figure 10: Star formation histories of two simulated populations with constant ( upper panel) and exponentially decreasing ( lower panel) star formation rates. The vertical axis is the SFR normalized to of total mass of stars formed. The blue dashed lines represents the simulated star formation. The red squares are the direct solution of the fit to a 3 SSP model. The green crosses are the results from MC simulations (200 inversions with ). 

Open with DEXTER 
We fitted the direct solutions (without noise) and we made a series of 200 inversions with different realizations of the added noise, corresponding to . Figure 10 presents the results as SFR versus time. For both simulated populations, we find a significant contribution (more than 10% in light) in the youngest component box, where the star formation had the lowest rate. The SFR reconstructions reproduce the simulated histories well and the MC simulations agree with the direct solution within the uncertainties. The direct solution underestimates the error bars due to the degeneracies between the ages and the relative weights of each component.
The error bars on the age are larger than in the case of NGC 205. The individual solutions of the MC realizations span the whole age range allowed for each of the components. This may be an indication that by contrast, the star formation history of NGC 205 is marked by a recent short burst rather than a smoothly declining SFR.
5 Conclusion
We presented the ULySS packages and its applications to the analysis of stellar atmosphere and stellar population spectra. This package is simple to use and is an efficient tool to determine the atmospheric parameters of stars ( , and [Fe/H]). The convergence region and the degeneracies can be studied in detail and the errors on estimated parameters are robustly determined. ULySS is also used to recover the history of star formation in galaxies and stellar clusters by decomposing the observed spectrum as a series of SSPs. The simultaneous analysis of the kinematics and of the stellar mix of a population allows us to break degeneracies and increase the reliability and precision of both the kinematics and the star formation history.
ULySS is available for download (http://ulyss.univlyon1.fr). Beside the applications described in the present paper, it contains other components (e.g. LINE, used to fit emission lines) and can easily be extended to other applications.
Acknowledgements
M.K. acknowledges a grant from the French embassy in Sofia and Y.W. acknowledges a grant from China Scholarship Council. We are grateful to Craig Markwardt (MPFIT) and Michele Cappellari (pPXF & BVLS) who freely distribute IDL/GDL routines which made this project possible. We thank also David Fanning (graphic library) and the contributors to the IDL Astronomical library. We acknowledge the help of Nicolas Bavouzet, Paul Blondé, Igor Chilingarian, Maela Collobert, and Martin France in the development and tests of the method over the years. We thank the anonymous referee for constructive comments that helped to improve the manuscript.
References
 BailerJones, C. A. L., Irwin, M., Gilmore, G., & von Hippel, T. 1997, MNRAS, 292, 157 [NASA ADS] (In the text)
 Bender, R. 1990, A&A, 229, 441 [NASA ADS]
 Bolstad, J. 1977, VARPROA general nonlinear leastsquares fitting code, Stanford Univ. (In the text)
 Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] (In the text)
 Butler, D. J., & MartínezDelgado, D. 2005, AJ, 129, 2217 [NASA ADS] [CrossRef] (In the text)
 Cappellari, M., & Emsellem, E. 2004, PASP, 116, 138 [NASA ADS] [CrossRef] (In the text)
 Cayrel de Strobel, G., Soubiran, C., & Ralite, N. 2001, A&A, 373, 159 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Cepa, J., & Beckman, J. E. 1988, A&A, 200, 21 [NASA ADS] (In the text)
 Chilingarian, I., Prugniel, P., Sil'Chenko, O., & Koleva, M. 2007, in IAU Symp. 241, ed. A. Vazdekis, & R. F. Peletier, 175
 Cid Fernandes, R., Mateus, A., Sodré, L., Stasinska, G., & Gomes, J. M. 2005, MNRAS, 358, 363 [NASA ADS] [CrossRef] (In the text)
 Coelho, P., Barbuy, B., Meléndez, J., Schiavon, R. P., & Castilho, B. V. 2005, A&A, 443, 735 [NASA ADS] [CrossRef] [EDP Sciences]
 da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Davidge, T. J. 2003, ApJ, 597, 289 [NASA ADS] [CrossRef] (In the text)
 de Bruyne, V., Vauterin, P., de Rijcke, S., & Dejonghe, H. 2003, MNRAS, 339, 215 [NASA ADS] [CrossRef] (In the text)
 Franx, M., Illingworth, G., & Heckman, T. 1989, ApJ, 344, 613 [NASA ADS] [CrossRef]
 Golub, G., & Pereyra, V. 1973, SIAM, 10, 413 (In the text)
 Heavens, A. F., Jimenez, R., & Lahav, O. 2000, MNRAS, 317, 965 [NASA ADS] [CrossRef] (In the text)
 Hekker, S., & Meléndez, J. 2007, A&A, 475, 1003 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Katz, D., Soubiran, C., Cayrel, R., Adda, M., & Cautain, R. 1998, A&A, 338, 151 [NASA ADS] (In the text)
 Koleva, M. 2009, Ph.D. Thesis, University of Lyon
 Koleva, M., Bavouzet, N., Chilingarian, I., & Prugniel, P. 2007, in Science Perspectives for 3D Spectroscopy, ed. M. KisslerPatig, J. R. Walsh, & M. M. Roth, 153
 Koleva, M., Gupta, R., Prugniel, P., & Singh, H. 2008a, in Pathways Through an Eclectic Universe, ed. J. H. Knapen, T. J. Mahoney, & A. Vazdekis, ASP Conf. Ser., 390, 302
 Koleva, M., Prugniel, P., & De Rijcke, S. 2008b, Astron. Nachri., 329, 968 [NASA ADS] [CrossRef] (In the text)
 Koleva, M., Prugniel, P., Ocvirk, P., Le Borgne, D., & Soubiran, C. 2008c, MNRAS, 385, 1998 [NASA ADS] [CrossRef]
 Koleva, M., De Rijcke, S., Prugniel, P., Zeilinger, W. W., & Michielsen, D. 2009, MNRAS, in press [arXiv:0903.4393]
 Lawson, C. L., & Hanson, R. J. 1995, Solving Least Squares Problems, Classics in Applied Mathematics, No. 15 (Philadelphia, Penn.: SIAM) (In the text)
 Le Borgne, D., RoccaVolmerange, B., Prugniel, P., et al. 2004, A&A, 425, 881 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Luck, R. E., & Bond, H. E. 1980, ApJ, 241, 218 [NASA ADS] [CrossRef]
 Marquart, D. W. 1963, SIAM, 11, 431 [CrossRef]
 Michielsen, D., Koleva, M., Prugniel, P., et al. 2007, ApJ, 670, L101 [NASA ADS] [CrossRef]
 Moré, J. J., Garbow, B. S., & Hillstrom, K. E. 1980, User Guide for MINPACK1, Report ANL8074, ANL, ANL
 Moultaka, J. 2005, A&A, 430, 95 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Munari, U., Sordo, R., Castelli, F., & Zwitter, T. 2005, A&A, 442, 1127 [NASA ADS] [CrossRef] [EDP Sciences]
 Neilson, H. R., & Lester, J. B. 2008, A&A, 490, 807 [NASA ADS] [CrossRef] [EDP Sciences]
 Ocvirk, P., Pichon, C., Lançon, A., & Thiébaut, E. 2006a, MNRAS, 365, 74 [NASA ADS] [CrossRef]
 Ocvirk, P., Pichon, C., Lançon, A., & Thiébaut, E. 2006b, MNRAS, 365, 46 [NASA ADS] [CrossRef]
 Panter, B., Heavens, A. F., & Jimenez, R. 2003, MNRAS, 343, 1145 [NASA ADS] [CrossRef] (In the text)
 Prugniel, P., & Koleva, M. 2007, in Spectral Line Shapes in Astrophysics, ed. L. C. Popovic, & M. S. Dimitrijevic, AIP Conf. Ser., 938, 27
 Prugniel, P., & Soubiran, C. 2001, A&A, 369, 1048 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Prugniel, P., Chilingarian, I., & Popovic, L. C. 2005, Mem. Soc. Astron. It. Supp., 7, 42
 Prugniel, P., Koleva, M., Ocvirk, P., Le Borgne, D., & Soubiran, C. 2007a, in IAU Symposium, ed. A. Vazdekis, & R. F. Peletier, IAU Symp., 241, 68 (In the text)
 Prugniel, P., Soubiran, C., Koleva, M., & Le Borgne, D. 2007b[arXiv:astroph/0703658] (In the text)
 Prugniel, P., Singh, H. P., Gupta, R., Wu, Y., & Koleva, M. 2009, in preparation (In the text)
 RecioBlanco, A., Bijaoui, A., & de Laverny, P. 2006, MNRAS, 370, 141 [NASA ADS] (In the text)
 Rix, H.W., & White, S. D. M. 1992, MNRAS, 254, 389 [NASA ADS]
 SánchezBlázquez, P., Peletier, R. F., JiménezVicente, J., et al. 2006, MNRAS, 371, 703 [NASA ADS] [CrossRef] (In the text)
 Sargent, W. L. W., Schechter, P. L., Boksenberg, A., & Shortridge, K. 1977, ApJ, 212, 326 [NASA ADS] [CrossRef] (In the text)
 Sarzi, M., FalcónBarroso, J., Davies, R. L., et al. 2006, MNRAS, 366, 1151 [NASA ADS] (In the text)
 Shkedy, Z., Decin, L., Molenberghs, G., & Aerts, C. 2007, MNRAS, 377, 120 [NASA ADS] [CrossRef] (In the text)
 Simien, F., & Prugniel, P. 2002, A&A, 384, 371 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Snider, S., Allende Prieto, C., von Hippel, T., et al. 2001, ApJ, 562, 528 [NASA ADS] [CrossRef] (In the text)
 Soubiran, C., Katz, D., & Cayrel, R. 1998, A&AS, 133, 221 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Soubiran, C., Bienaymé, O., Mishenina, T. V., & Kovtyukh, V. V. 2008, A&A, 480, 91 [NASA ADS] [CrossRef] [EDP Sciences]
 Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Monteiro, M. J. P. F. G. 2007, A&A, 469, 783 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Takada, M. 1977, PASJ, 29, 439 [NASA ADS]
 Takeda, Y. 2007, PASJ, 59, 335 [NASA ADS]
 Tojeiro, R., Heavens, A. F., Jimenez, R., & Panter, B. 2007, MNRAS, 381, 1252 [NASA ADS] [CrossRef] (In the text)
 Tonry, J., & Davis, M. 1979, AJ, 84, 1511 [NASA ADS] [CrossRef]
 Trager, S. C., & Somerville, R. S. 2009, MNRAS, 395, 608 [NASA ADS] [CrossRef] (In the text)
 Valdes, F., Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 251 [NASA ADS] [CrossRef] (In the text)
 van der Marel, R. P., & Franx, M. 1993, ApJ, 407, 525 [NASA ADS] [CrossRef] (In the text)
 Vazdekis, A. 1999, ApJ, 513, 224 [NASA ADS] [CrossRef]
 Vazdekis, A., Cenarro, A. J., Gorgas, J., Cardiel, N., & Peletier, R. F. 2003, MNRAS, 340, 1317 [NASA ADS] [CrossRef] (In the text)
 Venn, K. A. 1995, ApJS, 99, 659 [NASA ADS] [CrossRef]
 Willemsen, P. G., Hilker, M., Kayser, A., & BailerJones, C. A. L. 2005, A&A, 436, 379 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
 Wittkowski, M., Hummel, C. A., Aufdenberg, J. P., & Roccatagliata, V. 2006, A&A, 460, 843 [NASA ADS] [CrossRef] [EDP Sciences]
 Wu, Y. 2009, in preparation (In the text)
Footnotes
 ... Software)^{}
 ULySS is available at: http://ulyss.univlyon1.fr/
 ... MINPACK^{}
 http://www.netlib.org
 ... Vazdekis^{}
 http://www.iac.es/galeria/vazdekis/vazdekis_models_ssp_seds.html
 ... Language)^{}
 http://www.ittvis.com/idl/
 ... Language)^{}
 http://gnudatalanguage.sourceforge.net/
All Tables
Table 1: Stellar atmospheric parameters for six CFLIB stars of different spectral types.
Table 2: Ages and metallicities of the central 2'' of NGC 205.
All Figures
Figure 1: Effect of using a precise LSF, illustrated with a fit of a Vazdekis/Miles spectrum with a Pegase.HR/Elodie.3.1 SSP component. The top panel shows the spectrum (in black) and the best fit in blue (both are almost superimposed and the black line can be seen only when zooming on the figure), the red line is the multiplicative polynomial. The yellow regions were rejected from the fit (rejection of flagged telluric lines and automatic rejection of outliers). The middle and bottom panels are respectively the residuals from the best fit when (i) assuming a constant Gaussian LSF (in ) or (ii) a matched LSF. The continuous green lines mark the 1 deviation, and the dashed line is the zero axis. The right side of the figure expands a small wavelength region, around Ca4227. 

Open with DEXTER  
In the text 
Figure 2: Line spead function of the SDSS. 

Open with DEXTER  
In the text 
Figure 3: Comparison of the atmospheric parameters determined by ULySS with those from high resolution spectra (da Silva et al. 2006). The abscissas are the measurements from da Silva et al. (2006). On the top panels, the ordinates are from ULySS and the green lines are the diagonals. On the bottom panels the ordinates are the differences ULySS literature. 

Open with DEXTER  
In the text 
Figure 4: Fit of the CFLIB star HD 76151 with a TGM component. The symbols and conventions are the same as in Fig. 1. The order of the multiplicative polynomial is n=200. The right side expands a small wavelength region around Mg_{b}. 

Open with DEXTER  
In the text 
Figure 5: The evolution of the stellar atmospheric parameter fit results (, log( ), , and [Fe/H]) with increasing Legendre polynomial degree for 6 example CFLIB stellar spectra. Black, green, red, blue, violet and gold colors are for each star listed in Table 1. 

Open with DEXTER  
In the text 
Figure 6: Convergence maps on different projections of the parameters space for the CFLIB star HD 76151 inverted with the TGM component. Red crosses stand for the global minimum solutions found by ULySS. 

Open with DEXTER  
In the text 
Figure 7: MonteCarlo simulations and maps for the CFLIB star HD 76151 inverted with the ELODIE library as presented in Fig. 4. The three projections of the parameters space are presented. The 1000 MonteCarlo simulations are performed adding a random noise equivalent to the estimated one. The superimposed crosses give the internal errors estimated by ULySS, and the ellipses the standard deviation computed from the simulations. 

Open with DEXTER  
In the text 
Figure 8: Fit result for the central 2 of NGC 205 with 2 SSP components. In the top panel the black line represents the observed spectra, the blue line is the fitted model and the red line shows the multiplicative polynomial ( ). The bottom panel shows the residuals ( black), mean ( green dashed) and 1 deviation ( solid green). In both panels, some of the data were not considered for the fit ( yellow). 

Open with DEXTER  
In the text 
Figure 9: MonteCarlo simulation ( dots) and map ( contours) for a 2 CMP fit to spectra of the inner 2 of NGC 205. The blue contours represent the levels for the young stellar population (age constrained to be <800 Myr) and the red contours are for the old population (age between 800 and 14 000 Myr). 

Open with DEXTER  
In the text 
Figure 10: Star formation histories of two simulated populations with constant ( upper panel) and exponentially decreasing ( lower panel) star formation rates. The vertical axis is the SFR normalized to of total mass of stars formed. The blue dashed lines represents the simulated star formation. The red squares are the direct solution of the fit to a 3 SSP model. The green crosses are the results from MC simulations (200 inversions with ). 

Open with DEXTER  
In the text 
Copyright ESO 2009