Quantitative estimates of the constraints on solarlike models imposed by observables^{⋆}
^{1}
Observatoire de Paris, LESIA, CNRS UMR 8109, 92195 Meudon, France
email: Nesibe.Ozel@obspm.fr
^{2}
Institut d’Astrophysique et de Géophysique de l’Université de Liège, Allée du 6 Août 17, 4000 Liège, Belgium
email: MA.Dupret@ulg.ac.be
Received: 17 June 2009
Accepted: 3 May 2011
Context. Seismic parameters such as the large Δ_{0} and small δ_{02} frequency separations are now being measured in a very large number of stars and begin to be used to test the physics of stellar models.
Aims. We estimate the influence of different observed quantities (oscillation frequencies, interferometry, etc.) and the impact of their accuracy in constraining stellar model parameters.
Methods. To relate the errors in observed quantities to the precision of the theoretical model parameters, we analyse the behaviour of the χ^{2} fitting function around its minimum using the singular value decomposition (SVD) formalism. A new indicator called “weighting” quantifies the relative importance of observational constraints on the determination of each physical parameter individually. These tools are applied to a grid of evolutionary sequences for solarlike stellar models with varying age and mass, and to a real case: HD 49933 – a typical case for which seismic observations are available from space using CoRoT.
Results. The mass ℳ is always the best determined parameter. The new indicator “weighting” allows us to rank the importance of the different constraints: the mean large separation Δ_{0}, the radius R/R_{⊙}, the mean small separation δ_{02}, the luminosity L/L_{⊙}, the effective temperature T_{eff}. If the metallicity and age parameters are known, for example in an open cluster, using either individual or mean frequency separations yields the same uncertainties for masses less than 1.1 M_{⊙}. For HD 49933 the combination of ℳ and Y_{0}: ℳ ^{2}Y_{0} is well determined because of their correlation. However, they are poorly constrained individually. The frequency difference δ_{01}, if known with an error of about 0.3%, can determine the size of the convective core overshooting with about 3% accuracy.
Key words: methods: numerical / stars: fundamental parameters / stars: interiors / stars: oscillations
Appendices A and B are available in electronic form at http://www.aanda.org
© ESO, 2011
1. Introduction
The most promising way to probe the internal structure of a star is asteroseismology. The Sun has shown the way with helioseismology: observations of pmodes in the Sun have greatly improved solar modelling. Presently, pmode oscillations are being detected in solarlike stars. From the ground, the organization of coordinated campaigns has been successful (Bouchy & Carrier 2002; Carrier & Bourban 2003; Bedding et al. 2004; Kjeldsen et al. 2005; Fletcher et al. 2006; Bazot et al. 2007). From space, the results are now pouring in on many types of stars. MOST (Walker et al. 2003) and WIRE (Hacking et al. 1999) have detected oscillations in classical pulsators, while CoRoT (Baglin & The CoRoT Team 1998) and Kepler missions (Borucki et al. 1997) are observing the pmode spectrum in many solarlike stars.
These quantities, independent of the distance of the star, have a very important diagnostic power. However, to be really efficient, seismic data have to be complemented as much as possible by classical measurements, such as surface abundances, effective temperature, and surface gravity, all of which are obtained by spectroscopy. For nearby stars, radii and luminosities are sometimes measurable by interferometry, photometry, and the astrometry. For some binary systems, especially for eclipsing binaries, masses and radii are available with high accuracy.
To be able to derive the greatest amount of information from the data of the observed frequencies, many authors have developed diagnostic tools for probing the stellar interior and constraining the model parameters (e.g., ChristensenDalsgaard 1984; Gough 1987). These tools are the large frequency separation Δ_{n,ℓ}, and the small frequency separation δ_{n,ℓ}, defined as (1)These quantities were inferred from the lowdegree modes asymptotic expression. The physical interpretation of these seismic parameters can be found in Tassoul (1980) and Gough & Novotny (1990). According to the asymptotic theory, the large separation is proportional to the characteristic frequency of a star Δ_{0} (or Δ_{n,ℓ = 0}): Δ_{0} ∝ (M/R^{3})^{1/2}. Thus, in a homogeneous sequence of stars, the different characteristic frequencies of the oscillation spectrum are scaled by Δ_{0}, which is closely related to the star’s mean density (). The small separation depends on the variation in the sound speed in the central part of the star and may provide information about the composition of the star in its nuclear region. Consequently, the large and small separations can be combined to give estimates of both the mass and the age of a field star (Ulrich 1986; ChristensenDalsgaard 1984, 1988), and if other information about the star (for example its luminosity L) is available, then something might be learned about the convective mixing length or the composition (Gough 1987).
Thanks to the space missions CoRoT and Kepler, we have moved from ambiguous detections to firm measurements of seismic data in solarlike stars. It is therefore of interest to estimate the influence of constraints from seismic data and those obtained by photometry, spectroscopy, and interferometry, and to study in detail the effect of their precision on the determination of the fundamental stellar structure parameters like mass, age, and chemical composition. The knowledge of these basic parameters of stellar structure for a suitable sample of stars has a direct bearing on our understanding of the age and chemical evolution of the Galaxy.
To investigate the potential utility of all of the observations in an optimal way for determining the stellar parameters, we use a formalism based on singular value decomposition (SVD). This formalism relates errors in observed quantities to those in model parameters. In this context, we here address the following questions:

What is the importance of classical and seismic constraints on the determination of stellar parameters?

What is the behaviour of the χ^{2} function, defined by differences between observations and model, around its minimum in the parameter space?

Are the parameters correlated?

What precision can be obtained on each of these parameters for given errors on the observables?
2. Methods
The first step in studying the behaviour of the χ^{2} space consists of computing the best model that satisfies a set of observational constraints, called the reference model (RM).
2.1. The reference model (RM)
Given a set of n measurements y_{obs,i} (e.g. effective temperature, luminosity, seismic constraints) with associated error (σ_{i}), we first determine the bestfit model, or the RM, depending on a set of m parameters x_{j} (e.g. age, mass, etc.). This RM minimizes the χ^{2} fitting function defined as (2)where x = { x_{j} } (j = 1,...,m) is the parameter vector and y_{the,i}(x) (i = 1,...,n) are the model functions. The model functions, which we call observables, are the theoretical predictions for the observed quantities y_{obs,i} such as effective temperature T_{eff}, luminosity L, etc.
There are several local approaches available to find this minimum (Press et al. 1992), two notable ones are the LevenbergMarquardt (LM) algorithm and the singular value decomposition (SVD) method, etc. To automatically minimize the χ^{2} fitting function, we use the LM algorithm. This algorithm is more robust than the SVD method because it is designed for cases where the convergence is difficult. The LM method is a combination of the Newton and the steepest descent methods. This algorithm behaves like a Newton method when convergence is good and like the steepest descent method when convergence is difficult.
There are also algorithms that allow us to find the global minimum, for example the genetic algorithms as applied to white dwarf stars by Metcalfe & Charbonneau (2003). An advantage of our local approach is that the link between the different constraints and free parameters is naturally established through the evaluation of the derivatives, (∂y_{the,i}/∂x_{j}), at each step of the iteration. Problems with this method can arise when there are several local minima in the same region of the HR diagram. Therefore, it is impossible to find the global minimum and convergence problems may occur. In this case, either one can use global methods or construct a stellar grid and analyse the behaviour of χ^{2} function near the minima. Our study is restricted to wellconstrained solarlike stars with useful precision of seismic (oscillation frequencies) and classical information (spectroscopy, etc.). We assume a fixed identification of the oscillation modes and that they are not very far away from the asymptotic regime. In this case, the local approach works very efficiently, particularly when one has a good idea beforehand about the χ^{2} behaviour. It is then possible to restrict the parameter domain to a region where there is only one local minimum.
2.2. Erroruncertainty problem
To analyse the behaviour of the χ^{2} around its minimum, we use the singular value decomposition SVD method. This method was first introduced in the context of asteroseismology by Brown et al. (1994) and was subsequently used by several authors (e.g. Creevey et al. 2007). We recall the basic lines because of its importance in this paper. We first take the first order Taylor expansion of the model function or observable y_{the,i}(x) in the neighbourhood of a reference set of parameters x_{j0} that minimize the χ^{2} fitting function: (3)where y_{the,i0} is the set of observables resulting from the reference set of parameters x_{j0}, δx_{j} = x_{j} − x_{j0} and the derivatives ∂y_{the,i}/∂x_{j} are evaluated at x_{j} = x_{j0}. To determine the precision obtained on the model parameters and to show which combinations of parameters are well or poorly determined, we analyse the behaviour of χ^{2} around its minimum.
After substituting the expression (3) into Eq. (2) and some manipulation, the behaviour of χ^{2} around its minimum can be expressed thus (4)where χ_{min} is a minimum at x_{0} (x_{0} = x_{j0}) and D is called the design matrix (Brown et al. 1994): (5)This matrix relates small changes in the parameters to corresponding changes in the observables. This error analysis for solarlike stars is carried out using the SVD method.
Any n × m dimensional matrix D may be decomposed as (6)Inserting expression (6) into Eq. (4) and carrying out the matrix multiplications, one finds an m dimensional ellipsoidal equation, as shown in Brown et al. (1994) (7)where V^{(1)},...,V^{(m)} denote the columns of the matrix V. In the SVD method, the columns of V are an orthonormal set of m vectors that are the principal axis of the error ellipsoid Δχ^{2} = 1, while the corresponding values of W^{1} are the lengths of these axes.
The major advantage of the method is that it is more convenient to analyse the error ellipsoid in the parameter space, which gives information concerning the origin of the uncertainties on the obtained parameters. A vector V^{j} corresponding to a small singular value describes a direction in which χ^{2} varies little. However, χ^{2} increases very rapidly in the direction V^{j}, corresponding to a high singular value. If W_{j} is very high, the combination V^{1,j}δx_{1} + V^{2,j}δx_{2} + ... is welldetermined. However, if W_{j} is low, the combination V^{1,j}δx_{1} + V^{2,j}δx_{2} + ... is very poorly determined.
The estimation of the variancescovariances matrix of the free parameters due to the measurement errors on the n observables is expressed by (8)The diagonal elements of Cov give the uncertainty in each of the parameters: (9)To rule out the ambiguity between the observable errors and the uncertainties in the parameters, we denote the former by σ and the latter by ϵ.
The variancecovariance matrix quantifies the correlation between the parameters j and k, given by (10)If r_{jk} is close to 1, then the parameters jk are highly correlated. This means χ^{2} is partly degenerate and the parameters j and k are poorly constrained individually.
2.3. Weighting and significance
Many studies (e.g., Brown et al. 1994; Creevey et al. 2007; Ozel et al. 2010) have already asked the question of the role of the constraints on the results of the fitting. The usual approach consists of defining a quantity that estimates the importance of each observable on the global solution, called the “significance”. This has been already done several times, in particular by Brown et al. (1994). Here we try to go a step further and evaluate the importance of the observables on the determination of each parameter of the solution.
Weighting of an observable: the approach used here allows us to quantify the importance of each observable on the determination of parameters. It evaluates the impact of a given observable to constrain each individual parameter on the given confidence limit. To measure the weighting of a given observable y_{obs,i}, we compare two situations, completing the SVD analysis with and without a given observable. If the uncertainty on the jth parameter ϵ(x_{j}) increases abruptly when we remove this constraint, we conclude that this observable has a large weight and is important in constraining this parameter. Otherwise, it has no effect in determining the parameter. Furthermore, this approach gives a quantitative error estimate that exactly shows how the uncertainties of the parameters change when one removes/includes one of the constraints.
Significance: the “significance” approach has been proposed by Brown et al. (1994). The significance points out the importance of each observable for the full parameter solution. It is defined by the following equation: (11)In the remainder of the paper, the term “significance” represents this indicator. It can be understood as follows. Let us assume that for a set of observables y_{obs,0} we have found the minimum of χ^{2} at x_{0}. Then, we continuously modify the value of the ith observable y_{obs,i} while keeping the other observables constant. At each new value of y_{obs,i}, we can associate a set of parameters x minimizing the χ^{2} function. While doing this, we move from x_{0} to a point such that Δχ^{2} = χ^{2}(x_{1}) − χ^{2}(x_{0}) = 1 (recall that the ellipsoid Δχ^{2} < 1 can be interpreted as a confidence region, see Press et al. 1992, Sect. 14.5: Δχ^{2} as a function of confidence level and degrees of freedom). We have (12)and thus, (13)Hence, (14)We can distinguish two cases.

1.
S_{i} ≪ 1. This implies δy_{obs,i} ≫ σ_{i} (Eq. (14)). It is very unlikely to encounter a measurement several standard deviations away from the mean value. A modification of the constraint δy_{obs,i} ≃ σ_{i} would imply a negligible displacement of the minimum compared to the size of the Δχ^{2} < 1 domain, which means that this observable does not significantly affect the location of the minimum;

2.
S_{i} ~ 1. Here, a variation δy_{obs,i} ~ σ_{i} moves the minimum near the error ellipsoid edge (Δχ^{2} = 1), such an observable significantly affects the location of the minimum.
3. Stellar models and observables
3.1. Computation of stellar models and their seismic properties
A model of a solarlike star can be described with five free parameters: mass ℳ, age τ, mixing length parameter α to describe the outer convective zone, initial helium abundance Y_{0}, and initial ratio between heavyelement abundance and hydrogen Z/X_{0}. The reference values α and Y_{0} are set to the values obtained by calibration of the Sun (see Table 1) and Z/X_{0} corresponds to the solar abundances of Grevesse & Noels (1993).
Models are computed using CESAM2k (Morel & Lebreton 2008) for each parameter set (ℳ ,τ,α,Y_{0},Z/X_{0}). The code uses the OPAL equation of state and OPAL96 opacities (Iglesias & Rogers 1996), completed at low temperatures with the opacities of Alexander & Ferguson (1994). The physical description of the convective transport is the standard mixinglength theory (MLT, BöhmVitense 1958). Diffusion is not included. The oscillation frequencies are calculated with LOSC (Scuflaire et al. 2008) for modes of n = 15−25 and ℓ = 0−3, which could correspond to the detectable ones in the space asteroseismology era. CESAM2k calculates the following observables: effective temperature T_{eff}, luminosity L/L_{⊙}, metallicity Z/X, radius R/R_{⊙}, and seismic quantities: mean and individual large and small frequency separations Δ_{0}, δ_{02}, Δ_{0,i}, and δ_{02,i}.
Parameters of the models.
Observables y_{i,obs} and expected errors σ_{i}.
The model parameters are given in Table 1. The observables and their expected standard errors that define the reference models are given in Table 2. Depending on the evolutionary state of the model, the values of the observables y_{i,obs} vary between the intervals shown in the second column. The expected errors of the reference models σ_{i} are given in percent.
Each derivative of the matrix D is computed from differences of δx centred on the RM values x_{j0} given in Table 1, i.e., x_{j0} ± δx. The interval δx has to be sufficiently small so that the linear approximation is still valid, but also large enough to guarantee sufficient numerical accuracy. The increments δx_{j} for the grid of RM are reported in Table 3.
Increments δx_{j} for the RMs.
3.2. Observables and their errors
Spectroscopy provides effective temperature T_{eff}, surface gravity log g, and surface abundances. The intrinsic luminosity L_{⋆} is derived from the apparent magnitude, the bolometric correction, and the distance. If the distance is not known, one has to rely on photometric calibrations. If the star is an eclipsing binary, mass M and radius R can be determined at a high accuracy. If it is sufficiently close, the limbdarkened angular diameter θ_{ld} can be estimated, and combined with the parallax, can give a direct estimate of the radius. The oscillation frequencies are obtained using photometric and/or spectroscopic observations over a long time series.
The relative errors of the observables in Table 2 are chosen as the limit of the present instrumental techniques. The oscillation frequencies can be determined with an accuracy of about 0.1 μ Hz. This is representative of current CoRoT observations lasting for about 120 days. The number of large individual separations used as constraints is 23, and that of individual small separations is 6.
To investigate how seismic data contribute to determining physical stellar parameters, we consider three cases:

1.
using the mean large Δ_{0} and small δ_{02} frequency separations as seismic constraints;

2.
including the individual large and small separations (Δ_{0,i},δ_{i,02});

3.
reducing the errors on some of the observables by considering a future situation.
We decided not to use the individual oscillation frequencies for several reasons.

First, for a given absolute error on the frequenciesδν, the relative error δν/ν (used in our analysis) decreases as the frequency ν increases, which means that more weight would be given to higher frequencies. By high frequency we mean a frequency typically higher than the frequency of maximum power ν_{max}.

There are systematic errors in the estimation of theoretical frequencies. These errors can be caused by the adiabatic approximation or by not taking into account the dynamical effects of convection in the treatment of the physics of the outermost layers (ChristensenDalsgaard & Thompson 1997). These systematic errors are more dominant at high frequency. Therefore, high frequencies should not be given a large weight.

These systematic errors can be partly removed by using frequency separations.
It is accepted at present that the effective temperature T_{eff} is determined with a typical error of 100 K by photometry, and of about 50 K by spectroscopy (Bruntt 2009). σ [Fe/H], determined by a detailed spectroscopic analysis, may be of the order of 0.05 dex if the star is bright enough, does not rotate very rapidly (vsini < 25 km s^{1}), and has a spectrum with sufficient quality (Bruntt 2009). We choose σ [Fe/H] = 0.1 dex. The metallicity (Z/X) of the star is derived as . Consequently, the relative error on (Z/X)_{ ⋆ } is 23%. The accuracy of the surface gravity determinations is of the order of 0.2 dex (Kupka & Bruntt 2001).
The error on luminosity depends on the errors on the apparent magnitude (m_{V}), bolometric correction (BC), and parallax π. We can consider the relative error value on the parallax of 1% for bright stars, which translates to a relative error of 2% in luminosity.
The radius of a few single bright stars can be measured with an accuracy of about 1% (North et al. 2007) with interferometric methods or for eclipsing binaries. The error on radius depends on the errors on the parallax π and the angular size θ: D(D_{⊙}) = 107.47θ(mas)/π(mas). Stellar masses are measured directly for the eclipsing binary systems, with a relative error of 1%. However, they should be in general considered as unknown.
Table 2 lists the observables, including their relative errors in percent for the reference models. Scales and units used for errors on the observables and uncertainties on the parameters are expressed in logarithmic scale (relative error) so that errors are independent of the value of the measured observables and to facilitate comparison between errors in the different parameters. Using the logarithmic scale is not mandatory. However, it is the more appropriate choice if one does not want the results to depend on the units, the aim instead is to consider a broad part of the HR diagram, and the relative errors, not the absolute errors, are not expected to change much from one model to the other.
4. Results
We present two approaches for studying the parameter uncertainties. The first is a global analysis using all of the five free parameters and the second is restricted to wellchosen subparameter spaces.
Predicted uncertainties on the parameters ϵ(x).
The global analysis corresponds to real cases where all global parameters are unknown. However, it often appears in practical cases that the number of free parameters is larger than the number of observables. The SVD method is not designed for these underdetermined (n < m) and illposed problems. On the contrary, it is particularly powerful for overdetermined (n > m) systems. In our global analysis, there are five free parameters and six constraints. However, these constraints are not independent because the luminosity L, effective temperature T_{eff}, and radius R are related to each other by the StefanBoltzmann law. Thus, there are as many independent observables as there are free parameters. This means that the system is not overdetermined in the global analysis and it is not possible to use the “weighting” indicator.
Thus, for a detailed analysis of the relation between observables and parameters, we later restrict the parameter space to three parameters (for an analysis with two free parameters, see Ozel et al. 2010). This approach allows us to determine the “weighting” of each observable. We also note that some parameters are well known in some particular situations. For example, if the star is a member of a cluster, the age or the chemical composition could be well determined. Along the same lines as discussed in Brown et al. (1994), where the concept of pseudoobservable is introduced, the number of free parameters could be reduced in this situation. Furthermore, if one of the parameters is uncorrelated with all other parameters and has a fairly direct connection with one observable, then this parameter and the corresponding constraint can be safely ignored. A typical example is the metallicity, as we show below.
The difficulty related to the small number of observables compared to the number of parameters was partly overcome thanks to the new seismic observables. However, the situation here is quite critical. It can also be overcome by reducing the number of parameters and trying to reach a less complex situation.
Finally, reducing the number of parameters facilitates the interpretation of the results. To understand the connection between different parameters and observables, we take the three best constrained parameters, shown in black in Table 4, assuming that the other parameters are precisely known, as in open clusters. This is the approach adopted in this paper.
4.1. Global analysis
The relative uncertainties on the parameters ϵ(x) resulting from the SVD analysis for the global analysis are shown in Table 4. In this table we give the results for the models with the central hydrogen abundance X_{c} = 0.6, 0.4, and 02, corresponding to the beginning, the middle, and the end of the main sequence. The three bestdetermined parameters are shown in black. The precisions on the metallicity parameter are shown in grey; they are the same for all reference models. The least determined parameter is shown in red. These results are interpreted as follows.
4.1.1. Mass (ℳ)
In all analyses, ℳ is the bestdetermined parameter. For a relative error of 1% on the radius, we can achieve a precision of less than 3% on ℳ for models less than 1.3 M_{⊙}. For models grater than 1.3 M_{⊙}, the precision decreases towards more massive solarlike stars. This increase of ϵ(ℳ) comes from the correlation between ℳ and other parameters, particularly mixing length parameter α and age τ. For example, the correlations between ℳ and α (r_{ℳα} = 0.18), and ℳ and τ (r_{ℳ τ} = 0.33) for a 1.25 M_{⊙} star increase to r_{ℳ α} = 0.61, and r_{ℳ τ} = 0.69 for a 1.31 M_{⊙} star.
4.1.2. Mixing length parameter (α)
The relative uncertainty on this parameter ϵ(α) varies between 14 and 6% for all masses. It also decreases during the evolution for a 1.1 M_{⊙} star owing to the well constrained radius. For stars larger than 1.1 M_{⊙}, ϵ(α) increases. For a 1.55 M_{⊙} star, the uncertainty on this parameter is very large. This is because for more massive stars, M/M_{⊙} > 1.55, the change of α does not affect the general characteristics of the star (i.e., the radius R, the effective temperature T_{eff}), because the stellar structure becomes insensitive to the value of α. Thus, the behaviour of χ^{2} in the parameter space does not change in the direction corresponding to this parameter. It is therefore poorly determined.
4.1.3. Initial helium abundance (Y_{0})
The uncertainty on Y_{0}, ϵ(Y_{0}), lies between 10 and 13% for stars less than 1.4 M_{⊙} and is more or less constant for all considered evolutionary states. ϵ(Y_{0}) increases for a 1.55 M_{⊙} star as the star evolves.
In fact, the uncertainty of this parameter as a function of the mass does not show a very clear trend because the relationship between helium abundance and internal structure is indeed complicated. If we increase the abundance of hydrogen (and thus decrease Y_{0}),

the stellar structure changes in the deeper layers because thenuclear reaction rate increases;

the number of free electrons increases, which then increases the opacity (κ);

the mean molecular weight (μ) decreases, therefore, the equation of state changes.
This parameter is strongly correlated with the other parameters. The uncertainties on the other parameters (e.g., ϵ(α) increases with the mass) will therefore affect the uncertainty of Y_{0}. Thus, it is difficult to interpret the behaviour of ϵ(Y_{0}).
4.1.4. Initial metallicity (Z/X_{0})
Table 5 shows the coefficients of the correlation between Z/X_{0} and other parameters for a 0.9 M_{⊙} star. Z/X_{0} is not correlated with the best determined parameter, the mass ℳ. Thus, any improvement on the accuracy of the observables other than the metallicity (Z/X_{ ⋆ }) allows us to better constrain the mass parameter without affecting the precision on the metallicity parameter Z/X_{0}.
Coefficients of the parameters correlation.
It is normal that for a relative error on the constraint Z/X_{ ⋆ }, σ(Z/X_{ ⋆ })/(Z/X_{ ⋆ }) = 23%, we recover the same precision ϵ(Z/X_{0}) = 23% for all reference models, which are shown in grey in Table 4. As just mentioned, we deal in reality with as many independent constraints as free parameters. The constraint of the metallicity Z/X_{ ⋆ } determines the metallicity parameter Z/X_{0}, and independently (nocorrelation) the other constraints determine the other parameters. This is the case for all considered values of σ(Z/X_{ ⋆ }).
To have an overdetermined system, we remove one parameter in the analysis. For example, Table 6 shows the parameter uncertainties when we remove the initial helium abundance Y_{0} parameter. Obviously, here the uncertainty on the metallicity parameter is smaller than the error on the constraint. Thus, one needs more independent observables than parameters to improve the precision on this parameter.
Predicted uncertainties on the parameters ϵ(x).
4.1.5. Age (τ)
The uncertainty on the age is expressed as the absolute error normalized by the lifetime of the star on the main sequence (δτ/τ_{MS}). The precision on this parameter increases during the evolution particularly from a 1.1 M_{⊙} star.
It is well known that the small separation δ_{02} is sensitive to the evolutionary state of the star and consequently its age (ChristensenDalsgaard 1984). For a more detailed analysis, we reduce the relative error on this constraint by a factor of 2 (σ(δ_{02})/δ_{02} = 12.5 → 6.25%) for a reference model of 0.9 M_{⊙}. The results of the SVD analysis are shown in Table 7. A comparison with the results from Table 4 shows that improving the accuracy of δ_{02} mainly gives a better determination of the age parameter.
Predicted uncertainties on the parameters ϵ(x).
4.2. Subspace of the three bestconstrained parameters and relative importance of each observable
To quantify the “weighting” of each observable, we give the uncertainties of the subspaces of the parameters for a reference model of 0.9 M_{⊙} in Table 8 by including or excluding a given constraint. In Table 8, for example, the combinations of the varying parameters during the evolution of this model are ℳ, Y_{0}, and α, assuming that Z/X_{0} and τ are fixed (i.e., precisely known). The results for each reference model are shown in Tables A.1 − A.3 in the online material. In section ⓢ_{A} in these tables we present the uncertainties when we include all observables Y_{obs} for each reference model in the SVD analysis. In section ⓢ_{B} we include the individual small and large separations (Δ_{0,i},δ_{02,i}) instead of their mean values (Δ_{0},δ_{02}). In section ⓢ_{C} we remove each constraint one by one in the SVD analysis to determine its weighting.
The major results are the following:

1.
Mean large separation Δ_{0}: it has the largest weight on the determination of the mass parameter ℳ. Indeed, removing this constraint strongly increases the uncertainty of the mass ϵ(ℳ). For example, for a 0.9 M_{⊙} star in Table 8, ϵ(ℳ) increases from 2.70% (general case: using all observables Y_{obs} in the SVD analysis) to 759.77% (removed case: removing each constraint one by one in the SVD analysis).
Table 8Predicted uncertainties on the parameters ϵ(x).
In addition, there is a strong correlation between ℳ and other parameters (r_{ℳ,Y0} ≃ r_{ℳ,α} ≃ r_{ℳ,τ} ~ 1). Therefore, removing this constraint not only reduces the precision of ℳ but also the precision of the other parameters because of the correlation.

2.
Radius R/R_{⊙}: with a relative error σ(R/R_{⊙})/(R/R_{⊙}) = 1%, it is one of the most important constraints after the mean large separation Δ_{0}. The error on the radius has an important weight on the determination of ℳ, α, and Y_{0} as the uncertainties of these parameters increase by a factor of about 2 when this constraint is removed. Its impact on the determination of ℳ is larger when it is combined with Δ_{0} because .

3.
Mean small separation δ_{02}: with a relative error of about 13%, the exclusion of this constraint has no influence on the determination of the parameters. However, with a good precision of about 7% on δ_{02} , as illustrated in Table 7, it has an important weight on the age τ.

4.
Luminosity L/L_{⊙} and effective temperature T_{eff}: with a relative error σ(L/L_{⊙})/L/L_{⊙} ~ 2% and σ(T_{eff})/T_{eff} ~ 0.9%, L/L_{⊙} and T_{eff} have much less weight than Δ_{0} and R/R_{⊙} on the determination of the parameters.

Individual (Δ_{0,i},δ_{02,i}) instead of mean (Δ_{0},δ_{02}): both give approximately the same results on the parameter uncertainties for models less than 1.1 M_{⊙}. For example, for a 1.1 M_{⊙} star and X_{c} = 0.6 in Table A.2, ϵ(ℳ) = 2.65% and 2.25%; ϵ(Y_{0}) = 6.26% and 5.30%; ϵ(α) = 3.28% and 2.96% in section ⓢ_{A} and in section ⓢ_{B}, respectively. This is because we are close to the asymptotic regime where the individual large separations are close to the mean values and the errors on the small separations are too large to give a significant additional constraint. For models greater than 1.1 M_{⊙}, the precision on the parameters is better with individual frequency separations because we are farther from the asymptotic regime. For example, for a 1.25 M_{⊙} star and X_{c} = 0.6 in Table A.2, ϵ(ℳ) = 2.53% and 1.29%; ϵ(Y_{0}) = 6.11% and 2.95%; ϵ(α) = 3.14% and 2.43% in section ⓢ_{A} and in section ⓢ_{B}, respectively.
5. Application to a real case: HD 49933
HD 49933 is a typical example of a single field star, in which strong classical constraints do not exist, but for which the seismic data are now available from the CoRoT satellite.
5.1. Observables and their errors
Mosser et al. (2005) have detected solarlike oscillations in this star during a tennight observational run with the HARPS spectrometer. The first spacebased highquality photometric results obtained by CoRoT were presented by Appourchaux et al. (2008). The data consist of 42 individual frequencies of degree ℓ = 0,1 and 2, lying between 1.2 and 2.4 μHz. The average large separation Δ_{0} was found to be 85.9 ± 0.15 μHz derived from the degree l = 0 and l = 2 modes. We estimated the arithmetic mean frequency difference δ_{01} = ν_{n,0} − (ν_{n,1} + ν_{n − 1,1})/2 = −0.48 ± 0.32, by using the set of frequencies reported by Appourchaux et al. (2008). We do not consider here the small separation δ_{02} because the associated error bars are too large.
HD 49933 has a measured iron abundance of [Fe/H] = −0.37 dex (Solano et al. 2005), lower than the Sun and Procyon. An effective temperature of 6780 ± 130 K was determined by Bruntt et al. (2008), who also found [Fe/H] = −0.46 ± 0.08 dex. More recently, Gillon & Magain (2006) found a quite similar value, which we adopted here, [Fe/H] = −0.37 ± 0.03 dex. With the solar value Z/X_{⊙} = 0.0245 given by Grevesse & Noels (1993), it gives the metallicity Z/X_{ ⋆ } and the error σ_{Z/X ⋆ } given in Table 10. Bruntt (2009) reanalysed the star with highquality spectra and determined more accurate atmospheric constraints (T_{eff} = 6570 ± 60 K, log g = 4.28 ± 0.06 dex, [Fe/H] = −0.44 ± 0.03 dex).
The absolute visual magnitude derived by Appourchaux et al. (2008) is based on the revised Hipparcos catalogue by van Leeuwen (2007). Using the bolometric correction from Bessell et al. (1998), they estimated the luminosity to be L/L_{⊙} = 0.53 ± 0.01. Mosser et al. (2005) proposed a mass around 1.15 M_{⊙} by matching theoretical evolution tracks to the observed L − T_{eff} error box in the colour–magnitude diagram.
The reference model of HD 49933 was calculated by Goupil et al. (2011) using the fullspectrumofturbulence treatment (FST Canuto & Mazzitelli 1992). The model indicates that core overshooting is necessary. The overshooting parameter (d_{ov} = 4.02) has been included in the study.
The other input physics are identical to those adopted for the reference models of solarlike stars. The adiabatic oscillation frequencies are calculated for the modes of degrees ℓ = 0 − 2, n = 13 − 27. In Table 9 we show the derivative (or design) matrix of the computed reference models of HD 49933. The complete design matrix is given in Table B.1 in the online material.
Derivative matrix.
5.2. Results
Table 10 summarizes the observational constraints and the characteristics of the reference model of HD 49933; the first column gives the values of the constraints and the corresponding measurement errors used as input for the determination of the reference model by the LevenbergMarquardt algorithm; the second column gives the theoretical values of these constraints and the six model parameters (ℳ , Y_{0}, α, Z/X_{0}, τ, d_{ov}) for the reference model.
Given the uncertainties with the identification of the modes, we discuss three cases:

Case 1: including the individual large separations(Δ_{0,i});

Case 2: using all of the individual large seperations (Δ_{0,i}) and the small frequency differences (δ_{01,i});

Case 3: reducing the errors on the seismic data considered in Case 2, by factors of 5 and 10. This accuracy can be achievable by the missions such as Kepler (Borucki et al. 1997), SONG (Grundahl et al. 2008).
Table 11 presents the parameter matrix V of the SVD solution for case 1 (left panel) and case 3 (right panel) where the errors on the seismic data are reduced by a factor of 10 in case 3. The columns are ordered from the highest to the lowest singular value in Table 11. Table 12 shows the uncertainties on the parameters for HD 49933 for a global analysis in cases 1, 2, and 3.
Observations of HD 49933 and the properties of the RM.
Parameter matrix V for case 1 and case 3 of HD 49933.
Predicted uncertainties on the parameters ϵ(x).

The bestdetermined parameters ℳ and Y_{0}. Thecombination of parameters corresponding to the first column of V(Table 11) is best determined. This vector mainlylies in the ℳY_{0} plane both in case 1 and 3. This shows that the combination 0.88lnℳ + 0.42lnY_{0}, or equivalently ℳ ^{2.1}Y_{0} is determined with highest precision. The relative uncertainty on this combination is given by δ(ℳ^{2.1}Y_{0})/(ℳ^{2.1}Y_{0}) = 1/(0.42W_{1}) = 0.038% for case 1 and 0.0036% for case 3. However, the individual uncertainties on these parameters are much larger, particularly in case 1 (9.48 and 22.48%, Table 12), because these two parameters also have significant components in the poorly determined vectors (e.g., V^{5} for case 1, see Table 11). In addition, these two parameters are very highly correlated in case 1 : r_{ℳY0} = −0.99. These two effects are smaller in case 3, as shown by the significant decrease of their uncertainties from case 1 to case 3 in Table 12.

The α and τ parameters. Note that the α parameter for HD 49933 is described by the fullspectrumofturbulence treatment (FST Canuto & Mazzitelli 1992). In case 1 in Table 11 these parameters (α, τ) appear as dominant components in V^{2} and V^{3}, which explains their relatively small uncertainties. However they do not appear as a dominant component in any vector V^{i} in case 3. This means that the error ellipsoid is oblique with respect to this direction. As a consequence, these parameters are not well constrained individually.

The initial metallicity Z/X_{0}. The results are similar to those of the grid of reference models in Sect. 4.1.4. With the current precision in case 1, the Z/X_{ ⋆ } observable determines Z/X_{0} and the other observables constrain independently the other parameters. The comparison of case 1 and case 3 shows a slight decrease of the uncertainty on Z/X_{0}, indicating that if seismic data become very precise, they could help improving the precision on the metallicity parameter.

The overshooting parameter d_{ov}. The uncertainty on d_{ov} significantly decreases from case 1 to case 2 and then to case 3 (ϵ(d_{ov}) = 37.8, 2.96, and 0.77%, respectively, see Table 12). The huge improvement from case 1 to case 2 reflects the importance of the frequency difference information δ_{01,i} as an indicator of the extension of the convective core. This tendency is well known, e.g., Goupil et al. (2011) for HD 49933, Miglio & Montalbán (2005) for α Cen system, but we quantify it here.
6. Conclusions
Using classical methods (χ^{2} minimisation and SVD analysis), we have studied the characteristics of the χ^{2} surface close to its minimum. The analysis of the axes of the error ellipsoid Δχ^{2} = 1 helps to understand the relation between observable errors and uncertainty of the parameters. In addition to a classical global evaluation using the “significance”, we also use a new indicator, the “weighting”, which is able to quantify the relative importance of the observational constraints on the determination of each physical parameters individually.
In the mass and age range we examined, the “weighting” analysis ranks the relative importance of the different constraint: mean large separation Δ_{0}, radius R/R_{⊙}, mean small separation δ_{02}, luminosity L/L_{⊙}, and finally effective temperature T_{eff}. However, the degree of their importance is different and varies depending on the mass, the evolutionary stage, and the considered combination of parameters. For models greater than 1.1 M_{⊙}, with an accuracy of 0.1 μHz, the importance of Δ_{0} decreases because of the deviation of the frequencies from the asymptotic regime. With a relative error of 1%, the radius R has an important weight on the determination of ℳ, α, and Y_{0} for models less than 1.3 M_{⊙}.

The mean large separation Δ_{0} plays an important role in constraining the mass parameter ℳ. By combining the two constraints Δ_{0} and R, the relative precision on ℳ is significantly improved because .

With a relative error of about 13% on the mean small separation (δ_{02}), this constraint has no significant influence on the determination of the parameters. However, for a relative error of about 7%, the precision obtained on the parameters is significantly improved, particularly for the age.

With a relative error on the effective temperature σ(T_{eff})/T_{eff} ~ 0.9% and the luminosity σ(L/L_{⊙})/L/L_{⊙} ~ 2% , these constraints have no significant weight in determining the parameters when R, Δ_{0}, and δ_{02} are available. However, we must not conclude from this that the nonseismic constraints T_{eff} and L/L_{⊙} are no longer required in the asteroseismology era. When the mode identification remains uncertain and/or few frequencies are detected with large error bars, the nonseismic constraints such as T_{eff}, L/L_{⊙}, and log g remain crucial for the analysis, particularly for distinguishing between multiple solutions.

As for most stars of CoRoT and Kepler, radius and mass constraints are not available. Such a “real case” represented by HD 49933 shows that the new seismic (Δ_{0,i},δ_{01,i}) constraints, associated with the classical ones (T_{eff},L/L_{⊙},Z/X_{ ⋆ }), allow us to reach a precision of 6% on the mass, 5% on the age, 5% on the mixing length parameter α, 3% on the overshooting parameter, and about 13% on the initial helium abundance Y_{0}.

The combination of ℳ and Y_{0}: ℳ ^{2}Y_{0} is well determined because of their correlationfor all the considered cases of HD 49933. However, they are poorly constrained individually.

As already known, in stars with a convective core the major indicator of the size of the overshooting region is the seismic quantity δ_{01}. This allows the uncertainty of the overshooting parameter d_{ov} to be determined better than 3% if we have a precision of factors of 5 or 10 better than the actual precision on δ_{01}.

The seismic constraints with current precision cannot constrain the metallicity better than the spectroscopic value.
The type of the analysis described here can also be used to select the best targets for future seismic missions and to define the most relevant follow up observations from the ground.
References
 Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baglin, A., & The CoRoT Team 1998, in New Eyes to See Inside the Sun and Stars, ed. F.L. Deubner, J. ChristensenDalsgaard, & D. Kurtz, IAU Symp., 185, 301 [Google Scholar]
 Bazot, M., Bouchy, F., Kjeldsen, H., et al. 2007, A&A, 470, 295 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bedding, T. R., Kjeldsen, H., Butler, R. P., et al. 2004, ApJ, 614, 380 [NASA ADS] [CrossRef] [Google Scholar]
 Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
 BöhmVitense, E. 1958, Zeitschrift fur Astrophysik, 46, 108 [NASA ADS] [Google Scholar]
 Borucki, W. J., Koch, D. G., Dunham, E. W., & Jenkins, J. M. 1997, in Planets Beyond the Solar System and the Next Generation of Space Missions, ed. D. Soderblom, ASP Conf. Ser., 119, 153 [Google Scholar]
 Bouchy, F., & Carrier, F. 2002, A&A, 390, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brown, T. M., ChristensenDalsgaard, J., WeibelMihalas, B., & Gilliland, R. L. 1994, ApJ, 427, 1013 [NASA ADS] [CrossRef] [Google Scholar]
 Bruntt, H. 2009, ArXiv eprints [Google Scholar]
 Bruntt, H., De Cat, P., & Aerts, C. 2008, A&A, 478, 487 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Canuto, V. M., & Mazzitelli, I. 1992, ApJ, 389, 724 [NASA ADS] [CrossRef] [Google Scholar]
 Carrier, F., & Bourban, G. 2003, A&A, 406, L23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 ChristensenDalsgaard, J. 1984, in Liège International Astrophysical Colloquia, 25, Liège International Astrophysical Colloquia, 155 [Google Scholar]
 ChristensenDalsgaard, J. 1988, in Advances in Helio and Asteroseismology, ed. J. ChristensenDalsgaard, & S. Frandsen, IAU Symp., 123, 295 [Google Scholar]
 ChristensenDalsgaard, J., & Thompson, M. J. 1997, MNRAS, 284, 527 [NASA ADS] [Google Scholar]
 Creevey, O. L., Monteiro, M. J. P. F. G., Metcalfe, T. S., et al. 2007, ApJ, 659, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Fletcher, S. T., Chaplin, W. J., Elsworth, Y., Schou, J., & Buzasi, D. 2006, MNRAS, 371, 935 [NASA ADS] [CrossRef] [Google Scholar]
 Gillon, M., & Magain, P. 2006, A&A, 448, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gough, D. 1987, Nature, 326, 257 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O., & Novotny, E. 1990, Sol. Phys., 128, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Goupil, et al. 2011, submitted [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Perfectionnement de l’Association Vaudoise des Chercheurs en Physique, 205 [Google Scholar]
 Grundahl, F., Arentoft, T., ChristensenDalsgaard, J., et al. 2008, J. Phys. Conf. Ser., 118, 012041 [NASA ADS] [CrossRef] [Google Scholar]
 Hacking, P., Lonsdale, C., Gautier, T., et al. 1999, in Astrophysics with Infrared Surveys: A Prelude to SIRTF, ed. M. D. Bicay, R. M. Cutri, & B. F. Madore, ASP Conf. Ser., 177, 409 [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., Butler, R. P., et al. 2005, ApJ, 635, 1281 [NASA ADS] [CrossRef] [Google Scholar]
 Kupka, F., & Bruntt, H. 2001, J. Astron. Data, 7, 8 [Google Scholar]
 Miglio, A., & Montalbán, J. 2005, A&A, 441, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Mosser, B., Bouchy, F., Catala, C., et al. 2005, A&A, 431, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 North, J. R., Davis, J., Bedding, T. R., et al. 2007, MNRAS, 380, L80 [NASA ADS] [CrossRef] [Google Scholar]
 Ozel, N., Dupret, M., & Baglin, A. 2010, Ap&SS, 143 [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing (Cambridge University Press), 2nd ed. [Google Scholar]
 Scuflaire, R., Montalbán, J., Théado, S., et al. 2008, Ap&SS, 316, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Solano, E., Catala, C., Garrido, R., et al. 2005, AJ, 129, 547 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Ulrich, R. K. 1986, ApJ, 306, L37 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Results of ϵ(x) to subspaces of RMs
Predicted uncertainties on the parameters ϵ(x).
Predicted uncertainties on the parameters ϵ(x).
Predicted uncertainties on the parameters ϵ(x).
Appendix B: Derivative matrix for HD 49933
Derivative matrix