Asteroseismic model fitting by comparing ϵ_{nℓ} values
Astronomy Unit, Queen Mary University of London,
Mile End Road,
London
E1 4NS,
UK
email:
I.W.Roxburgh@qmul.ac.uk
Received: 25 May 2015
Accepted: 26 October 2015
We present an asteroseismic model fitting algorithm based on comparing model and observed ϵ_{ℓ}(ν) values defined in terms of frequencies by ν_{nℓ} = Δ [ n + ℓ/ 2 + ϵ_{ℓ}(ν_{nℓ}) ] where Δ is an average large separation. We show that if two stellar models have the same interior structure but different outer layers then the difference between their ϵ_{ℓ}(ν) values, interpolated to the same frequencies, collapses to a function only of frequency, independent of angular degree ℓ. The algorithm tests the goodness fit by comparing the difference in model and observed ϵ values after having subtracted off a best fit ℓ independent function of frequency ℱ(ν), and only requires interpolation in model values and not in observed values so the errors on the observed values are uncorrelated; it is independent of the n values assigned to the radial ordering of the frequencies and does not require the calculation of inner phase shifts of the model. We contrast this to a proposed direct frequency matching technique which minimises the difference between observed and model frequencies after having subtracted off an ℓ independent fit to these differences. We show this technique is flawed in principle, that all models with the same dimensionless structure but any mass and radius have the same quality of fit to an observed data set, and that it can give erroneous best fit models. We illustrate the epsilon matching technique by comparing stellar models and then apply it to data on HD 177153 (aka Perky). On comparing observations with a set of main sequence evolutionary models we find that models which satisfy constraints on the luminosity, radius, Δ, and on ϵ matching, have masses in the range 1.155 ± 0.035 M_{⊙} and ages in the range 4.486 ± 0.250 × 10^{9} yr. Since the large separation and the radius are surface layer dependent we examine “pure surface layer independent” model fitting where the only constraints on the model fitting are on the luminosity and epsilon matching, and show that the best fit models have M/M_{⊙} = 1.13 ± 0.06 and age = 4.62 ± 0.39 × 10^{9} yr.
Key words: stars: oscillations / asteroseismology / stars: interiors / methods: analytical / methods: numerical
© ESO, 2015
1. Introduction
The seemingly straight forward way to find a model that has the same internal structure as an observed star is to seek to match the model and observed frequencies. The problem that besets such an approach is that it requires accurate modelling of the outer layers of a star, since these layers make make a significant contribution to the values of the frequencies. However our understanding of the physics of these layers is poor, e.g. modelling convection, the equation of state, diffusion/levitation, mixing ..., which renders it difficult to reliably model these outer layers.
In previous papers we have presented two model fitting techniques which overcome this difficulty by seeking to subtract off the unknown contribution of the outer layers and just fit the interior structure: comparison of the ratio of small to large separations (Roxburgh & Vorontsov 2003a, 2013; Roxburgh 2005) and phase matching (Roxburgh & Vorontsov 2003b; Roxburgh 2015a). The former suffers from having to interpolate in the observed frequencies and from correlation of errors, the latter from having to calculate inner phase shifts of the stellar models.
We here present a technique based on the result that if a star and a model have the same interior structure then the difference between their ϵ_{ℓ}(ν) values, evaluated at the same frequencies, collapses to a function only of frequency and not of angular degree ℓ. The ϵ are defined in terms of frequencies ν_{nℓ} by ν_{nℓ} = Δ(n + ℓ/ 2 + ϵ_{nℓ}) where Δ is an average large separation. This procedure only requires knowledge of observed and model frequencies, and interpolates only in model frequencies so the errors on observational values are uncorrelated.
2. Properties of the frequencies and ϵ_{nℓ}
A set of frequencies ν_{nℓ} of a (slowly rotating) mainsequence star can be represented in terms of “phases” ϵ_{nℓ} as (1)where Δ is some arbitrarily chosen reference large separation usually taken as an average value Δ = ⟨ ν_{nℓ} − ν_{n − 1,ℓ} ⟩. For a given Δ the ϵ_{ℓ}(ν) are therefore known at the discrete set of frequencies ν_{nℓ}, but clearly their values depend the choice of Δ. An example is given in Fig. 1 for Model A of mass 1.15 M_{⊙}, with 14 ℓ = 0,1,2 frequencies in the range 1500−3000 μHz with an average large separation Δ = 104.3 μHz.
Model A, and all other models used in this paper, were calculated using the STAROX code (cf. Roxburgh 2008a), with EOS5 equation of state (Rogers & Nayfonov 2002), NACRE reaction rates (Angulo et al. 1999), OPAL/Whichita opacities (Iglesias & Rogers 1996; Ferguson et al. 2005), GS98 relative abundances (Grevesse & Sauval 1998), the mixing length model of convection, an Eddington atmosphere out to an optical depth of 10^{3}, and do not include diffusion or convective overshooting. Frequencies were calculated using the OSCROXL code (Roxburgh 2008b). Model A has an initial hydrogen abundance X = 0.72, heavy element abundance Z = 0.021, mixing length parameter α = 1.8, no core overshooting (β = 0), and is evolved to a central hydrogen abundance X_{c} = 0.096 at an age of 4.6 × 10^{9} yr.
Fig. 1 ϵ_{nℓ} values for the frequencies of a mainsequence model of 1.15 M_{⊙} (Model A). The continuous curves ϵ_{ℓ}(ν) are obtained by interpolation in the discrete values. 

Open with DEXTER 
The ϵ_{nℓ} measure the departure of the structure of a model from that of a uniform sphere and are determined by integrals over the structure of the star, the dominant contributions coming from the outer layers and the deep interior. As shown in previous work (cf. Roxburgh & Vorontsov 1994, 2000, 2013; Roxburgh 2015a, and references therein) these two contribution can be separated by solving the equations governing the oscillations of a star in terms of inner phase shifts δ_{ℓ}(ν,t) which satisfy the central boundary conditions and are determined by the structure of the inner layers, and outer phase shifts α_{ℓ}(ν,t) which satisfy the outer boundary conditions and are determined by the structure of the outer layers, where t is the acoustic radius ^{∫}dr/c. Both α_{ℓ}(ν,t) and δ_{ℓ}(ν,t) are continuous functions of both t and ν. The eigenfrequencies ν_{nℓ} are determined by continuity of the solutions at any arbitrarily chosen t and satisfy the eigenfrequency equation (Roxburgh & Vorontsov 2000) (2)where T is the total acoustic radius of the star.
The important properties of this phase shift representation, on which all surface layer independent model fitting techniques rest, is that in the outer layers of a star, the outer phase shifts α_{ℓ}(ν,t) are almost independent of ℓ and beneath the near surface layers are almost independent of t, and can therefore be replaced by a single ℓ independent function α(ν). This is illustrated in Fig. 2 where we plot the α_{ℓ}(ν,t) for Model A against acoustic radius t for 3 values of the frequency ν; the curves for ℓ = 0,1,2 at a given ν lie on top of each other until deep inside the star. The inner phase shifts are dependent on ℓ but are almost constant independent of acoustic radius in the outer layers beneath the near surface layers. (For a detailed analysis see Roxburgh 2015a.)
Suppressing the matching radius t, which can be anywhere in the outer layers, Eq. (2) can be written in the form (3)The surface phase shift, and hence the contribution of the outer layers, can be eliminated by interpolating the ϵ_{ℓ}(ν) values to the same frequency and subtracting so that (4)which only depend on the inner phase shift differences and hence on the inner structure of the star. The separation ratios (interpolated to the same frequency) are an approximation to these differences (cf. Roxburgh & Vorontsov 2003a,b, 2013).
Fig. 2 Superposed curves for inner phase shift differences δ_{ℓ0} from the model, ϵ differences ϵ_{0ℓ}(ν) from interpolation in the ϵ_{nℓ} values from the frequencies, and the separation ratios r_{010},r_{02} for Model A. 

Open with DEXTER 
Equation (3) relates the ϵ_{nℓ} to the value of the acoustic Δ_{T} but it can readily be transformed to any other chosen reference large separation Δ; suppressing the fitting radius t and taking α_{ℓ}(ν) = α(ν) independent of ℓ in Eq. (2), we have (5)and since ν_{nℓ} [ 1 − Δ / Δ_{T} ] is a function only of ν it can be absorbed into α(ν); setting (6)we have (on dropping the asterisk), (7)The inner phase shifts δ_{ℓ}(ν) as a function of ν, which are determined by the inner structure of the star, are unchanged by this redefinition of α.
3. Models with the same interior structure
Consider 2 models with frequencies , with the same interior structure but different outer layers; their ϵ values at their frequencies in terms of inner and outer phase shifts are If the two models have the same interior structure then their inner phase shifts δ_{ℓ}(ν) as a function of ν are identical but are evaluated at different frequencies , . If we interpolate in the values to determine the values at the frequencies then Since the two models have the same interior structure the inner phase shifts , subtracting Eq. (9b) from Eq. (9a) gives (10)where ℱ(ν) is an ℓ independent function of frequency.
So if the two models have the same interior structure then the difference in their ϵ values evaluated at the same frequency ℰ(ℓ,ν) collapses to an ℓ independent function of frequency, whereas if they do not have the same interior structure then their differences are a function of both ℓ and ν. This is the basis of the model fitting algorithm presented here.
4. Model fitting algorithm
Given a set of observed frequencies with error estimates σ_{nℓ} and the frequencies of a model, calculate their ϵ_{ℓ}(ν_{nℓ}) and interpolate in to find the model values at the observed frequencies (e.g. by cubic spline interpolation). Then form their differences (11)make an M parameter ℓ independent fit ℱ(ν) to ℰ(ℓ,ν) (e.g. by a series in Chebychev polynomials ∑ C_{k}T_{k}(ν), or in terms of basis Bsplines) and determine the goodness of fit by the reduced χ^{2}(12)The errors s_{nℓ} are the errors in where the σ_{nℓ} are the errors on the observed frequencies and Δ^{o} the large separation of the observed frequencies. N is the total number of observed frequencies and M is the number of parameters in ℱ(ν) which minimises the χ^{2}, subject to being less than the number of ℓ = 0 frequencies.
This algorithm does not require one to determine the inner and outer phase shifts, only to interpolate in the model’s ϵ_{ℓ} values which are known over a wide range of frequencies (minimising end effects in interpolation), and the errors s_{nℓ} are uncorrelated.
Fig. 3 Superposed curves for the outer phase shifts α_{ℓ}(ν,t) of Model A at 3 frequencies. The values of α_{ℓ} are independent of ℓ until deep inside the star. 

Open with DEXTER 
Fig. 4 Top panel: for Model G, continuous curves ϵ_{ℓ}(ν) and interpolated values for at the “observed” frequencies of Model A (crosses). Bottom panel: residuals of fit of to function of frequency ℱ(ν). 

Open with DEXTER 
5. Model examples
The first example is to compare Model A, (designated as the “observed” star) with frequencies and error estimates σ_{nℓ} = 0.15 μHz, with Model G, with frequencies , which has the same interior structure as Model A but is modified in the outer layers by setting Γ_{1} constant above r = 0.95R; their frequency differences are shown in Fig. 3. In the top panel of Fig. 4 we show the for Model G over a wider range of frequencies than the “observed” set, and their continuous curves ϵ_{ℓ}(ν) obtained by cubic spline interpolation. The interpolated values at the observed frequencies are marked by the crosses.
The residuals of the fit to a function only of ν are shown in the bottom panel; with error estimates on the “observed” frequencies of 0.15 μHz the reduced χ^{2} = 0.016, the departure from zero being primarily due to errors in interpolation, and the fact that the outer layer contribution is not exactly independent of ℓ.
Fig. 5 Top panel: for Model B, continuous curves ϵ_{ℓ}(ν) and interpolated values for at the “observed” frequencies of Model A (crosses). Bottom panel: residuals of fit of to function of frequency ℱ(ν). 

Open with DEXTER 
Our second example is to take a model similar to Model A and test how well it fits the epsilon values of model A. This model (Model B) has a mass of 1.13 M_{⊙}, XH = 0.70,Z = 0.022,α = 1.8,X_{c} = 0.047 and an age of 4.52 × 10^{9} yr. The top panel in Fig. 5 shows the ϵ values for the model and the values in interpolated at the frequencies of Model A, and the bottom panel the residuals to the fit to of the ϵ differences to a function only of frequency. With error estimates on the frequencies of model A of 0.15 μHz and N − M = 35 degrees of freedom, the χ^{2} of the fit is 0.412 which corresponds to a probability of matching of 0.999; this model fits the ϵ values of Model A.
The results in Fig. 5 illustrate the fact that ϵ matching and separation ratios are not exactly equivalent. The points in Fig. 5 are the ϵ differences ( where , whereas the difference between model and observed separation ratios is ( at . The resulting χ^{2} for the two procedures are not in general the same, which is the smaller depends on the particular model and observed frequencies. In this example whereas from ϵ matching χ^{2} = 0.412.
The third model to be tested (Model C) has a mass of 1.15 M_{⊙} with XH = 0.70,Z = 0.019,α = 2.4,X_{c} = 0.200 and an age of 3.6 × 10^{9} yr. The top panel in Fig. 6 shows the ϵ values for this model and the values interpolated at the frequencies of Model A, and the bottom panel the residuals of the fit of the ϵ differences to a function only of frequency. With error estimates on the frequencies of model A of 0.15 μHz; the χ^{2} of the fit is 3.16 which corresponds to a probability of matching of 3 × 10^{9}, this model does not fit the ϵ values of Model A.
Fig. 6 Top panel: ϵ_{nℓ} values of the model C, the continuous curves ϵ_{ℓ}(ν) obtained by interpolation, and the values at the frequencies ν^{o} of Model A (the crosses). Bottom panel: residuals of fit of to a function of frequency ℱ(ν). The χ^{2} of the fit is 3.16, demonstrating that this model does not have the same interior structure as model A. 

Open with DEXTER 
6. Model fitting by comparing frequencies rather than ϵ values
A seemingly simple alternative model fitting algorithm is to compare observed and model frequencies after having subtracted off an ℓ independent best fit ℱ(ν^{o}) to the frequency differences , defining the goodness of fit as (13)where, as in ϵ matching, N is the number of frequencies and M the number of parameters in the function of frequency ℱ(ν). This was proposed as an alternative to the use of the surface layer offset approach advocated by Kjeldsen et al. (2008), with the offset determined as part of the fitting process rather than by scaling the solar value (cf. Roxburgh 2009). This was misguided as I show below!
Consider a (homologous) family of stellar models, that is a family with different masses M, and radii R, but the same dimensionless structure q(x),s(x),p(x),Γ_{1}(x), where x = r/R and(14)This family of models all have the same dimensionless frequencies ω_{nℓ} and the same inner phase shifts δ_{ℓ}(ω_{nℓ}) as a function of dimensionless frequency, but the physical values of the frequencies differ since they are given by (15)the frequencies of any one model in the set being simply a linear scaling of those of any other model in the set.
Now if we have a model whose frequencies satisfy the matching condition in Eq. (13), then so do all other members of the homology set since the difference in frequencies between any such models is simply a linear function of frequency and is therefore absorbed into ℱ(ν) in Eq. (13), giving the same χ^{2} for all models in the set. So the frequency fitting algorithm (Eq. (13)) is not able to distinguish between stars with widely different M, R, which have the same dimensionless interior structure. The same problem applies to comparing separation ratios at the same radial order n, since the ratios – being ratios – are independent of a constant scaling on the frequencies.
This is not the case for the epsilon matching since we compare values at the same physical frequency. nor for comparison of separation ratios at the same frequency. There is one well known exception; all scaled models that have the same M/R^{3} have the same physical frequency and asteroseismology alone cannot differentiate between this set of models.
Moreover the frequency matching algorithm Eq. (13) cannot, in principle, give a χ^{2} = 0 even if the model and observed star have the same interior structure, unless the model and observed frequencies, , are a perfect match. From Eq. (5) we have (16)where is a function only of frequency given by (17)If the observed star and model have the same interior structure then the inner phase shifts have the same functional form so δ^{m}(ν) = δ^{o}(ν) at the same frequencies, but in general at different frequencies and the frequency differences are ℓ dependent, unless the observed and model frequencies are the same.
Fig. 7 Top panel: residuals of the fit of the frequencies of model C to those of Model A using the frequency matching algorithm (Eq. (13)) gives a . Bottom panel: internal structure of the models. 

Open with DEXTER 
As an example of a spurious good fit from frequency matching we apply the algorithm (Eq. (13)) to model C whose ϵ matching residuals with a χ^{2} = 3.16 are shown in Fig. 6; the residuals from frequency fitting are shown in Fig. 7 and have a χ^{2} = 0.77, naively suggesting that model C is a good fit to Model A. This is not the case as can be seen from Fig. 7b where we plot the sound speed squared and density as a fuction of radial distance from the centre (in units 10^{11} cm).
Fig. 8 ϵ_{nℓ} values of the observed frequencies of HD 177153 (aka Perky). The error bars shown are 5σ to enable them to be seen. 

Open with DEXTER 
However if the model is a good fit to the observations as determined by epsilon matching (), it may, or may not, also have value of , depending on whether or not the n values of the model and observed star are in agreement, or the change in the inner phase shifts δ_{ℓ} over a large separation is sufficiently small.
7. ϵ model fitting to HD 177153 (aka Perky)
The star HD 177153 (aka Perky) was observed by Kepler (KIC 6106415) and Silva Aguirre et al. (2013) list a total of 33 frequencies, 11 for each ℓ = 0,1,2, a ν_{max} ~ 2210 μHz and an average large separation Δ = 104.0 ± 0.5 μHz.
Model fitting was undertaken by SilvaAguirre et al. (2013) using two procedures: comparing frequencies with a “surface offset” (cf. Kjeldsen et al. 2008), and comparing separation ratios (cf. Roxburgh & Vorontsov 2003a, 2013); and more recently by Roxburgh (2015a) using the phase matching technique. Here we apply the ϵ matching technique described above to the same star. The ϵ_{nℓ} values are displayed in Fig. 8 – note that the error bars have been multiplied by 5 in order that they can be seen.
HD 177153: 1σ best fit models: Y ≥ 0.255.
As in Roxburgh (2015a) we take constraints on luminosity and radius derived from angular diameter and bolometric flux measurements by Huber et al. (2012) and the revised Hipparcos parallax (van Leeuwen 2007) giving We define an average large separation in terms of the frequencies, setting Δ = [ ν_{n + 3,0} − ν_{n − 2,0} ] / 5 where ν_{n,0} ≤ ν_{max}<ν_{n + 1,0}, since this can be applied to both observations and models, and take an error estimate of 1 μHz, larger than that given by Silva Aguirre et al. to allow for the uncertainty due to the relatively large contribution of the outer layers (Roxburgh 2014) and the fact that even in the solar case Δ from solar models differs from the observed value by ~1 μHz; this gives an observational constraint of Δ^{o} = 104 ± 1 μHz. Since we are using a surface layer independent model fitting algorithm it would be inconsistent to impose a constraint on the surface [Fe/H] (which is anyway poorly constrained) and is unlikely to reflect the interior composition due to diffusion/levitation/mixing in the outer layers.
We searched for best fit models in just one model set with Grevesse & Sauval (1998) abundances computed using my STAROX code (Roxburgh 2008a) as described in Sect. 2 above. The model set has masses in the range 0.85−1.40 M_{⊙}, initial hydrogen abundance XH in the range 0.68−0.74, metal abundance Z in the range 0.014−0.026, mixing length parameter α in the range 1.1−2.4, and evolutionary models from the initial to terminal main sequence. These models do not include convective overshooting or diffusion/levitation.
We define best fit models as those that satisfy the luminosity, radius and large separation constraints to within 1σ, and which have χ^{2} ≤ 1.114 for the epsilon fitting. The limit on χ^{2} corresponds to a probability of 0.317 of being due to a random realisation of errors for N − M = 24 degrees of freedom (i.e. N = 33 frequencies and M = 9 parameters in the representation of ℱ(ν)), the same as the 1σ limit for 1 degree of freedom. Likewise we set a 2σ limit on χ^{2} of 1.535 corresponding to a probability of 0.0455, the same as the 2σ limit for 1 degree of freedom. One additional constraint can be added; namely that the initial Helium abundance (Y) should be greater than the initial helium fraction from the Big Bang; we set a lower limit on Y of 0.255
In Table 1 we give models that all fell with the 1σ limits; the models have masses in the range 1.155 ± 0.035 M_{⊙} and ages in the range 4.486 ± 0.250 × 10^{9} yr (age9 is age in Gyr and Δ_{μ} the large separation in μHz); well within the target accuracy of 10% required for the PLATO mission (cf. Rauer et al. 2014). The first model below the line has a slightly small radius whilst the following two models are within 2σ for Δ – allowing for the fact that the value of Δ, and to a lesser extent the radius, are influenced by the structure of the outer layers of a star.
We have excluded models with low Helium abundance 0.234 ≤ Y< 0.255; however these models have similar masses and ages to those in Table 1 with 1.13 ≤ M/M_{⊙} ≤ 1.18 and ages in the range 4.64−4.89 × 10^{9} yr.
Fig. 9 Top panel: ϵ_{nℓ} values of a model of 1.15 M_{⊙}, the crosses mark the interpolated values at the observed frequencies of HD 177153. Bottom panel: residuals of fit of to ℱ(ν), a function only of frequency. 

Open with DEXTER 
Figures 9 give the results for one of the best fit models which has M = 1.15 M_{⊙}, XH = 0.72,Z = 0.021,α = 1.8, X_{c} = 0.096, an age of 4.6 × 10^{9} yr. and Δ = 104.2 μHz. The upper panel shows the ϵ_{nℓ} values for this model, the points are the model values, the curves are the result of fitting a cubic spline through these model values, and the crosses the interpolated values at the observed frequencies. The bottom panel gives the residuals of the fit of the differences to a function ℱ(ν) only of frequency, which has a reduced χ^{2} = 0.86. This model has an interior structure compatible with that of the observed star HD 177153.
It is worth pointing out that the frequency offsets – that is the difference between the observed frequencies and those of the best fit models , do not, in general, follow a Kjeldsenlike power law aν^{b}, and can be positive and/or negative. This is discussed in detail in Roxburgh (2015b). The one exception is the 1.150 M_{⊙} model below the line in Table 1 which does have a Kjeldsentype offset (cf. Roxburgh 2015b).
8. Relaxing the constraints – pure surface layer independent model fitting
In the above we have followed the procedure generally used in model fitting, requiring models to satisfy observational constraints on luminosity, radius (or effective temperature) and large separation, as well as fitting to the frequencies – here in the form of epsilon matching. But this is not actually surface layer independent model fitting. Whilst the luminosity L (and mass M) are essentially independent of the structure of the outer layers, being determined by the run of density, temperature and composition profiles in the fully ionised inner layers, the large separation Δ and radius R are strongly influenced by the structure of the outer layers which depend on the model of convection, diffusion and the equation of state in regions of ionisation. We therefore explore the consequences of relaxing constraints on Δ and R.
In Table 2 we give the results of such an analysis. An entry of 1 signifies that the constraint was imposed with a 1σ limit, ∞ that the constraint was not imposed, and 2 that the constraint was imposed with a 2σ limit. As mentioned above we take the 1σ limit for χ^{2} as 1.114, and the 2σ limit as 1.535. The first row corresponds to the models in Table 1 where all variables L, R,Δ,χ^{2} are within their 1σ limits; and rows 2–7 have different combinations of constraints on L,R,Δ and rows 8 and 9 have no constraint from ϵ matching. The last 3 rows below the line are for 2σ limits.
From the data in row 4 we see that “pure surface layer independent model fitting”, that is just imposing a 1σ constraint on the luminosity and on epsilon matching, is effective in narrowing down the range of possible models both in terms of age and mass giving M/M_{⊙} = 1.13 ± 0.06 and age = 4.62 ± 0.39 × 10^{9} yr.
Fig. 10 Models that satisfy different combinations of constraints – models that satisfy all constraints (the large red dots) have M/M_{⊙} = 1.155 ± 0.035 and ages 4.486 ± 0.250 × 10^{9} yr. 

Open with DEXTER 
Mass and age for vs fitting criteria Y ≥ 0.255.
It is clear from rows 8 and 9 in this table that ϵ matching is needed to narrow down the range of masses and ages. This is clearly illustrated in Fig. 10 where we plot all models that satisfy, to within 1σ, just the constraints on the luminosity and radius (the black points), those that satisfy to within 1σ the constraints on luminosity, radius and large separation (the blue points), and in red the best fit models from Table 1 that also satisfy the epsilon matching constraints.
9. Conclusions
We have shown that the if a model and observed star have the same interior structure then their ϵ_{ℓ}(ν) values, compared at the same frequency reduce to a function only of frequency and not of angular degree ℓ. We then give an algorithm for determining best fit models to an observed frequency set which only requires interpolation in model ϵ_{ℓ}(ν_{nℓ}) values to determine their values at the observed frequencies, and does not require interpolation in the observed values; as a result the observational errors in the comparison are uncorrelated. We also show that what was once thought to be a possible alternative algorithm, comparing frequencies after subtracting off a best fit function only of frequency, is flawed in principle, gives the same quality of fit for all models with any radii and masses in a family of models which have the same dimensionless structure, and can give erroneous results.
Applying the epsilon matching algorithm to the observed frequency set of HD 177153 (aka Perky) with observational constraints on L,R,Δ we obtain best fit models with masses in the range 1.155 ± 0.035 M_{⊙} and ages in the range 4.486 ± 0.250 × 10^{9} yr, which are consistent with the results from phase matching given in Roxburgh (2015a), the 3 best fit models being identical in the two analyses.
Further, on imposing only “pure surface layer independent constraints”, that is only constraints on L and ϵ, we find best fit models with M/M_{⊙} = 1.13 ± 0.06 and age = 4.62 ± 0.39 × 10^{9} yr, consistent with the values obtained on imposing a full set of constraints.
Acknowledgments
The author gratefully acknowledges support from the Leverhulme Foundation under grant EM2012035/4 and from the UK Science and Technology Facilities Council (STFC) under grant ST/M000621/1.
References
 Angulo, C., & the NACRE consortium 1999, Nucl. Phys. A, 656, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Ferguson, J. W., Alexander, D. R., Allard, D. T., et al. 2005, ApJ, 623, 585 [NASA ADS] [CrossRef] [Google Scholar]
 Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Huber, D., Ireland, M., Bedding, T., et al. 2012, ApJ, 760, 32 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., & ChristensenDalsgaard, J. 2008, ApJ, 683, L175 [NASA ADS] [CrossRef] [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W. 2005, A&A, 434, 655 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W. 2008a, Ap&SS, 316, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W. 2008b, Ap&SS, 316, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W. 2009, Presentation at meetings, Synergies between solar and stellar modelling Rome, June 22–26, New Insights into the Sun, the potential of a new generation of solarstellar diagnostics, Ponte de Lima, Portugal, September 16–18 [Google Scholar]
 Roxburgh, I. W. 2014, A&A, 571, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W. 2015a, A&A, 574, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W. 2015b, A&A, 581, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 1994, MNRAS, 268, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2000, MNRAS, 317, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003a, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003b, Ap&SS, 284, 187 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2013, A&A, 560, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Silva Aguirre, V., Basu, S., Brandao, I. M., et al. 2013, ApJ, 769, 17 [NASA ADS] [CrossRef] [Google Scholar]
 van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
All Figures
Fig. 1 ϵ_{nℓ} values for the frequencies of a mainsequence model of 1.15 M_{⊙} (Model A). The continuous curves ϵ_{ℓ}(ν) are obtained by interpolation in the discrete values. 

Open with DEXTER  
In the text 
Fig. 2 Superposed curves for inner phase shift differences δ_{ℓ0} from the model, ϵ differences ϵ_{0ℓ}(ν) from interpolation in the ϵ_{nℓ} values from the frequencies, and the separation ratios r_{010},r_{02} for Model A. 

Open with DEXTER  
In the text 
Fig. 3 Superposed curves for the outer phase shifts α_{ℓ}(ν,t) of Model A at 3 frequencies. The values of α_{ℓ} are independent of ℓ until deep inside the star. 

Open with DEXTER  
In the text 
Fig. 4 Top panel: for Model G, continuous curves ϵ_{ℓ}(ν) and interpolated values for at the “observed” frequencies of Model A (crosses). Bottom panel: residuals of fit of to function of frequency ℱ(ν). 

Open with DEXTER  
In the text 
Fig. 5 Top panel: for Model B, continuous curves ϵ_{ℓ}(ν) and interpolated values for at the “observed” frequencies of Model A (crosses). Bottom panel: residuals of fit of to function of frequency ℱ(ν). 

Open with DEXTER  
In the text 
Fig. 6 Top panel: ϵ_{nℓ} values of the model C, the continuous curves ϵ_{ℓ}(ν) obtained by interpolation, and the values at the frequencies ν^{o} of Model A (the crosses). Bottom panel: residuals of fit of to a function of frequency ℱ(ν). The χ^{2} of the fit is 3.16, demonstrating that this model does not have the same interior structure as model A. 

Open with DEXTER  
In the text 
Fig. 7 Top panel: residuals of the fit of the frequencies of model C to those of Model A using the frequency matching algorithm (Eq. (13)) gives a . Bottom panel: internal structure of the models. 

Open with DEXTER  
In the text 
Fig. 8 ϵ_{nℓ} values of the observed frequencies of HD 177153 (aka Perky). The error bars shown are 5σ to enable them to be seen. 

Open with DEXTER  
In the text 
Fig. 9 Top panel: ϵ_{nℓ} values of a model of 1.15 M_{⊙}, the crosses mark the interpolated values at the observed frequencies of HD 177153. Bottom panel: residuals of fit of to ℱ(ν), a function only of frequency. 

Open with DEXTER  
In the text 
Fig. 10 Models that satisfy different combinations of constraints – models that satisfy all constraints (the large red dots) have M/M_{⊙} = 1.155 ± 0.035 and ages 4.486 ± 0.250 × 10^{9} yr. 

Open with DEXTER  
In the text 