Issue 
A&A
Volume 663, July 2022



Article Number  A92  
Number of page(s)  12  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/202243640  
Published online  14 July 2022 
Probing stellar cores from inversions of frequency separation ratios
Observatoire de Genève, Université de Genève, Chemin Pegasi 51, 1290 Versoix, Switzerland
email: Jerome.Betrisey@unige.ch
Received:
25
March
2022
Accepted:
18
May
2022
Context. With the rapid development of asteroseismology thanks to spacebased photometry missions such as CoRoT, Kepler, TESS, and in the future, PLATO, and the use of inversion techniques, quasimodelindependent constraints on the stellar properties can be extracted from a given stellar oscillation spectrum. In this context, inversions based on frequency separation ratios, which are less sensitive to surface effects, appear as a promising technique to constrain the properties of stellar convective cores.
Aims. We developed an inversion based on frequency separation ratios, with the goal of damping the surface effects of the oscillation frequencies. Using this new inversion, we defined a new indicator to constrain the boundary mixing properties of convective cores in solarlike oscillators.
Methods. We verified our inversion technique by conducting tests in a controlled environment, where the stellar mass and radius were known exactly, and conducted an extensive hare and hounds exercise.
Results. The inversion is not affected by surface effects. With the construction of an extensive set of models, favoured and forbidden regions can be highlighted in the parameter space. If the ratios are well fitted, unsurprisingly the inversion does not provide any additional information.
Conclusions. The indicator, coupled with the inversion based on frequency separation ratios, seems promising at probing the properties of convective cores, especially for Ftype stars exhibiting solarlike oscillations.
Key words: stars: interiors / stars: fundamental parameters / asteroseismology
© J. Bétrisey and G. Buldgen 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the SubscribetoOpen model. Subscribe to A&A to support open access publication.
1. Introduction
In the last two decades, the spacebased photometry missions CoRoT (Baglin et al. 2009), Kepler (Borucki et al. 2010), and TESS (Ricker et al. 2015) were launched, enabling a rapid development of asteroseismology. Further progress is also expected from the future PLATO mission (Rauer et al. 2014). The quality of the data from these missions made possible the use of socalled seismic inversion techniques, which had previously been applied with tremendous success in helioseismology (see Buldgen et al. 2019a; ChristensenDalsgaard 2021, for recent reviews). The situation for asteroseismology is, however, very different from that of the Sun, as geometric cancellation forbids the detection of modes with angular degrees higher than three for solarlike oscillators. This limits the available data sets to about 50 frequencies at best, whereas thousands of modes are observed for the Sun. As a consequence, the philosophy of inversions was reoriented. Indeed, the classical inversions for the Sun are local inversions that can scan through its entire structure, but they require a large set of frequencies in order to be efficient. This can only be achieved for the best asteroseismic targets (e.g. 16CygA and 16CygB, Bellinger et al. 2017, 2019, 2021; Buldgen et al. 2022), and with the limitation that the information can only be extracted for the deep layers. Therefore, inversions based on socalled indicators were developed (Reese et al. 2012; Buldgen et al. 2015a,b, 2018), and applied to various cases (e.g. Buldgen et al. 2016a,b, 2017a, 2019b,c; Bétrisey et al. 2022). The idea was to concentrate all the information of the frequency spectrum to conduct a global inversion that constrains a global quantity, the indicator (e.g. the mean density).
In this work, we are interested in the entropy profile in the central adiabatic convective regions to put some constraints on the properties of the convective core. With the current inversions based on individual frequencies, a large number of them are required, especially at various harmonic degrees. In addition, a very accurate characterisation of the mass and radius is also needed (Buldgen et al. 2017b). Entropy inversions were conducted for the Sun (Buldgen et al. 2017c), and for the 16Cyg binary system (Buldgen et al. 2022) using the indicators defined in Buldgen et al. (2018). However, even though the 16Cyg system is one of the best Kepler asteroseismic targets, Buldgen et al. (2022) found that no additional information about the central entropy could be extracted for 16CygB. In addition, inversions based on individual frequencies are sensitive to the surface regions, which is a limiting factor given that the nearsurface effects are poorly modelled. In the case of 16CygA, the surface effects and the quality of the reference models prevented Buldgen et al. (2022) from extracting new constraints on the central entropy. Therefore, we chose to develop a new indicator based on frequency separation ratios instead of individual frequencies, with the aim of carrying out inversions of solarlike stars presenting convective cores. Indeed, frequency separation ratios are defined to damp surface effects (Roxburgh & Vorontsov 2003) and are more efficient at probing central stellar regions (Otí Floranes et al. 2005). Furthermore, Deheuvels et al. (2016) and Farnir et al. (2019) have shown that frequency separation ratios could be used to constrain the extent of convective cores in solarlike oscillators.
In Sect. 2, we recall the stateoftheart of indicatorbased inversions using individual frequencies, and introduce a new inversion based on frequency separation ratios. In Sect. 3, we define a new indicator that probes central stellar regions and test it in a controlled environment. In Sect. 4, we conduct a hare and hounds exercise to simulate how the indicator would behave with a real observed target. Finally, in Sect. 5, we draw the conclusions of our study.
2. Theoretical considerations
In the inversion community, the following terminology is generally applied. The inversion aims to retrieve the properties of an ‘observed’ or ‘target’ model, which can either be a real observation or a synthetic target. To this end, the inversion takes a ‘reference’ model as an input, which is the result of a modelling procedure, and gives a small correction that is added to a quantity of interest computed from the reference model as an output. For simplicity in this manuscript, an ‘inverted’ quantity is defined as the quantity of interest that includes the correction from the inversion.
2.1. Stateoftheart of individual frequencies
The structure inversion equation is based on a perturbative analysis of the stellar oscillations at linear order. This approach is motivated by the work of LyndenBell & Ostriker (1967) and their predecessors (see e.g. Chandrasekhar 1964; Chandrasekhar & Lebovitz 1964; Clement 1964), who showed that the equation of motion fulfils a variational principle. In the case of individual frequencies, the frequency perturbation is directly related to the structural perturbation (Dziembowski et al. 1990):
with ν the oscillation frequency, n the radial order, l the harmonic degree, a and b two structural variables, and the structural kernels, and using the definition
where the index ref stands for reference, and obs stands for observed. Historically, Eq. (1) was derived for the sound speed and density structural pair (c^{2}, ρ), but the equation can be generalised to any combination of physical quantities appearing in the adiabatic oscillation equations (Elliott 1996; Kosovichev 2011; Buldgen et al. 2017b).
Given the frequency differences, the idea is then to solve Eq. (1) to have access to the structure differences and compute a global integrated quantity, the socalled indicator t (e.g. the mean density), which concentrates all the information of the frequency spectrum. The general form of the indicator is
where f is a weight function that depends on the radius, and g and h are two functions of the structural variables. Except in rarer cases (see e.g. Buldgen et al. 2015a), h(b) is set to one, and the indicator is only a function of one of the structural variables. In practice, the options for g are limited, because the target function 𝒯_{t} in Eq. (4) must be reproducible with the number of modes available in asteroseismology. Therefore, simple functions such as g(a) = a, g(a) = 1/a, etc. are good candidates. A more detailed explanation of the determination of f and g is provided in Sect. 3, where we define the indicator used in this work.
This approach has the advantage of being quasimodel independent because it does not rely on the physics assumed to generate the reference model that serves as starting point for the inversion. In addition, the inversion is also quasiindependent from the starting point in the sense that, if we start from another point in the parameter space, the inversion will still correct towards the exact value.
The inversion can be performed with the subtractive optimally localised averages (SOLA) method (Pijpers & Thompson 1994). This minimisation technique is an adaptation of the OLA approach of Backus & Gilbert (1968, 1970), and allows the inversion to be conducted in a numerically more efficient way. In practice, for a given indicator t, the following cost function is minimised:
where the averaging and crossterm kernels are related to the structural kernels,
The normalisation k depends on the properties of the indicator. The balance between the amplitudes of the different terms during the fitting is adjusted with tradeoff parameters, β and θ. Here, the idea is to provide a good fit of the target function 𝒯_{t} while reducing the contribution of the crossterm and of the observational errors to the individual frequencies. We introduce a short notation for the identification pair of a frequency, i ≡ (n, l). The 1σ uncertainty of the relative frequency differences is denoted by σ_{i}, and is defined, with N the number of observed frequencies. The variable λ is a Lagrange multiplier, and c_{i} are the inversion coefficients.
The function ℱ_{Surf}(ν) is an empirical term to describe the surface effects. Besides introducing additional free parameters in the minimisation at the expense of the fit of the target function, the treatment of the surface effects constitutes a limitation of the inversion. Indeed, structural kernels based on individual frequencies have a large amplitude in the surface regions, as shown in the upper panel of Fig. 1, and are therefore sensitive to the treatment of these layers. In practice, this limitation can be translated into a systematic uncertainty (see e.g. Bétrisey et al. 2022), but it can become significant enough to prevent the inversion of from correcting accurately.
Fig. 1. Frequency kernels versus ratio kernels of the model 1.3MoAGSS from Sect. 3.1. For both panels, we show the kernel amplitude as a function of the position, and for the structural variable pair (c^{2}, ρ). Upper panel: loworder (l = 0, n = 11) and highorder (l = 0, n = 26) kernels for the individual frequencies, respectively in orange and blue. Lower panel: loworder (n = 11) and highorder (n = 26) kernels for the frequency separation ratios, respectively in orange and blue. 
2.2. Adaptation for frequency separation ratios
Roxburgh & Vorontsov (2003) introduced frequency separation ratios to damp surface effects. They are defined as the ratio of the small frequency separations divided by the large frequency separations. Because the modes with a higher harmonic degree than l = 2 are difficult to observe in asteroseismology, we only consider the following ratios:
where Δ_{l}(n) are the large separations, and d_{xy}(n) are the small separations:
The first version of a kernel of a ratio was presented in Otí Floranes et al. (2005), with the introduction of the r_{02} kernels, and modulo a different scaling than in this work. We present here the derivation of the r_{01} kernels. The derivations for the r_{10} and r_{02} kernels can be found in Appendix A. As for the individual frequencies, the relative differences of the ratios can be related to the structure differences
Because a frequency ratio is a linear combination of individual frequencies, its kernel is a linear combination of the structural kernels of the frequencies used in the definition of the ratio. For example, the r_{01} kernels are the result of the following combination:
The kernels are computed by switching a and b in the previous expression.
Figure 1 shows a comparison of the kernels for the individual frequencies and for the frequency separation ratios, for a 1.3 M_{⊙} star (model 1.3MoAGSS from Sect. 3.2). The frequency kernels have a large amplitude at the surface of the star, and are therefore well suited to probe this region, but sensitive to surface effects. On the contrary, ratio kernels have an amplitude that is damped at the surface, suppressing the surface effects. In addition, they have a larger amplitude in the central stellar regions, allowing these regions to be probed more efficiently. We chose to display the kernels of the structural pair (c^{2}, ρ) because it nicely illustrates this behaviour, but in practice, another structural pair should be used. Indeed, the form of these kernels is not adapted with our indicator, and the contribution of the crossterm is important for this structural pair. For inversions based on frequency separation ratios, structural pairs such as (S_{5/3}, Y) or (S_{5/3}, Γ_{1}) are better suited. The variable S_{5/3} = P/ρ^{5/3} is a proxy of the entropy, P is the pressure, ρ is the density, Y is the helium mass fraction, and is the first adiabatic exponent.
There are two main differences for the SOLA cost function. We cannot assume that the ratio uncertainties are uncorrelated because a ratio is, by definition, a combination of frequencies, and the surface effects can be neglected by construction (as we show later in Sect. 3.1).
As for the individual frequencies, β and θ are tradeoff parameters to balance the amplitudes of the different terms during the fitting. Again, the idea is to provide a good fit of the target function 𝒯_{t} while reducing the contribution of the crossterm and the observational errors. The variable λ is a Lagrange multiplier, and c_{p} are the inversion coefficients. The normalisation , where is the square of the Frobenius norm of the covariance matrix of the ratios, and N_{r} is the number of observed ratios. The entries of the covariance matrix of the ratios are given by Eq. (B.14) from Appendix B: .
3. A new indicator probing stellar cores
In this section, we aim to define an indicator that is sensitive to the properties of convective cores. For the first variable, the entropy proxy S_{5/3} appears to be a natural choice. Indeed, this proxy is motivated by the SackurTetrode equation for the entropy in the case of an ideal gas,
where k_{B} is the Boltzmann constant, μ is the mean molecular weight, m_{u} is the atomic mass unit, and g(μ) is a function that depends only on the mean molecular weight. They both form a plateau in adiabatic convective regions, and the height and extension of this plateau are related to the boundary mixing properties, such as the temperature and mean molecular weight gradients. Because these boundary regions are generally poorly modelled, an indicator sensitive to their properties can provide new constraints on the physical processes acting in these regions.
For the second structural variable, we chose the helium mass fraction. In comparison to Γ_{1}, the amplitude of the crossterm is significantly smaller, which naturally reduces the crossterm contribution in the SOLA cost function. Even if the structural pair including Γ_{1} is less well suited, an inversion based on these variables is feasible, although with a lessened stability. Indeed, we can define error measures to quantify the sources of error in an inversion (Buldgen et al. 2015b). The first sources of error are the errors on the averaging and crossterm kernels, which quantify how well the associated target functions are reproduced. The last source of error is the residual error that quantifies the unknown remaining error after taking into account the two previous error measures. In an ideal case, the averaging kernel would perfectly reproduce the target function, the crossterm would be zero, and the inversion could therefore retrieve the exact value. In practice, the crossterm error is negligible, and compensations occur when the averaging and residuals errors are of the same magnitude but with opposite sign. In such cases, the inversion is unstable because a significant unknown source of error affects the inversion result.
In the following section, we define a new indicator that probes central stellar regions, and test it in a controlled environment where the stellar mass and radius between the reference and target models are identical. In these conditions, the mean density is well controlled, which strengthens the stability of the inversion.
3.1. Definition and calibration
We define our indicator as follows:
where f(r) is a weight function. The minus sign is purely cosmetic to ensure that the indicator is positive. Because only a small number of modes are available in asteroseismology, the target function must be constructed so that it can be fitted by the structural kernels. Therefore, the weight function is chosen such that the resulting target function is close to the natural form of the kernels.
We compute the eulerian linear perturbation of the indicator:
Hence, the target function is:
The determination of the weight function is the most difficult part when developing a new indicator. Indeed, there is no straightforward procedure to determine it, except trial and error, and one has to find a good balance between the finetuning of the weight function and the tradeoff parameters. In practice, to find the weight function, we start with a guess for f(r) and iterate until we find a target function that is reproducible with the available kernels. There is no recipe for the initial guess, but we started with an r^{2} term (because the kernels follow approximately this behaviour in the very central regions), an exponential decay (to suppress the signal at the surface), and two Gaussian functions. After a few iterations, we found that including a third Gaussian function to match the number of ‘lobes’ of the target function (illustrated in Fig. 2), as well as the hyperbolic tangent (to more efficiently suppress the signal beyond the central region), produces better results. Because we use Y as structural variable, the crossterm is naturally damped and the value of β has very little impact on the inversion. To find an adequate balance, we calibrated the weight function and the tradeoff parameters for one model, the model 1.3MoAGSS of Sect. 3.2. Then, we verified the indicator in a controlled environment (Sect. 3.2) and with an hare and hounds exercise (Sect. 4). We found that adopting β = 10^{−4}, θ = 10^{−4}, and
Fig. 2. Target function and entropy profiles for the calibration model. The observed and modelled entropy profiles are in blue and orange, respectively. The target function is in green. 
where
is a good compromise in most of the cases.
In Fig. 2, we show the target function and the entropy profiles for the calibration model to illustrate where the information is localised and to which stellar regions the target function is the most sensitive. The structural differences are mainly located in the central region, up to onetenth of the stellar radius. Therefore, most of the contribution in the inversion will originate from this region. The first inner peak and the inner part of the second peak of the target function are the most important because they can ‘see’ where the structural differences are localised. Also, it is not necessary to efficiently reproduce the last peak of the target function because it is sensitive to a region where the structural differences are almost nonexistent. In addition, we recall that the kernels are basis vectors that are combined to reproduce the target function. From a linear algebra point of view, some of them will generate the global shape of the target function, while some others will reproduce the details. In our inversion, the global shape is generated by the loworder kernels and they have the highest coefficients. The highorder kernels retrieve some of the details and have smaller coefficients. In this sense, loworder kernels are thus an asset for the inversion. We note that extreme orders (n < 5 or n > 30) showed a numerically unstable behaviour. This might be a numerical effect that could be fixed or the sign of an intrinsic nonlinear behaviour. From a theoretical point of view, further investigations beyond the scope of this article are required. However, it is currently unrealistic to measure modes with n < 5, and highorder modes (n > 30), as well as being very difficult to observe, will be very sensitive to surface effects.
For the calibration model, we tested that the inversion is indeed not affected by surface effects. To this end, we added surface effects on the frequencies of the synthetic observed model using the Sonoi et al. (2015) prescription, which generates strong surface corrections, especially at high order, as shown in Fig. 3. Even though the frequencies are significantly affected, the impact on the ratios is very small and the relative ratio differences are modified on average by about 0.08%, and much less for loworder modes that carry the inversion. Because this effect is much smaller than any optimistic observational uncertainty, especially at low order, the impact of the surface effects on the inversion is completely negligible.
Fig. 3. Impact of the surface effects on the frequencies. The frequencies before (orange) and after (purple) the addition of the surface effects are shown for different harmonic degrees l. 
3.2. Behaviour in a controlled environment
For the preliminary tests in a controlled environment, we generated synthetic ‘observed’ targets in the main sequence with masses of 1.1 M_{⊙}, 1.2 M_{⊙}, 1.3 M_{⊙}, 1.4 M_{⊙}, and 1.5 M_{⊙}, and modified one physical ingredient at a time. For the observed targets, we used the OPAL opacities (Iglesias & Rogers 1996), supplemented by the Ferguson et al. (2005) opacities at low temperature and the electron conductivity by Potekhin et al. (1999). We used the GN93 abundances (Grevesse & Noels 1993), the FreeEOS equation of state (Irwin 2012), and microscopic diffusion was described using the formalism of Thoul et al. (1994), but with the screening coefficients of Paquette et al. (1986), and taking partial ionisation into account. The nuclear reaction rates were from Adelberger et al. (2011). We fixed the mixinglength parameter α_{MLT} at a solar calibrated value of 2.05, and following the implementation of Cox & Giuli (1968). For the atmosphere modelling, we used the T(τ) relation from Eddington (1959). For all the models in this work, we used the Liège Evolution Code (CLES, Scuflaire et al. 2008a), and the frequencies were computed with the adiabatic Liège Oscillation Code (LOSC, Scuflaire et al. 2008b).
For the reference models, we changed the abundances for the AGSS09 abundances (Asplund et al. 2009), the opacities for the OPLIB opacities (Colgan et al. 2016), and the atmosphere by using the T(τ) relation described by ModelC in Vernazza et al. (1981) (hereafter VALC). We also introduced overshooting and undershooting, respectively with α_{ov} = 0.1 and α_{under} = 0.3, and assuming an adiabatic stratification in both cases. We selected the reference models in the evolutionary sequence by matching the radius with the observation. Concerning the uncertainty of the relative ratios differences, we adopted an error of 20% (hereafter equalweighted uncertainties) that corresponds to the expectation for highquality data.
For the inversions, we considered several sets of modes, summarised in Table 1. As shown in Fig. 4, with the addition of loworder modes, the averaging kernel better reproduces the target function, especially in the central regions. Moreover, the contribution of highorder modes is small. Indeed, the target function of sets 1 and 5 is almost identical. For sets 3 and 4, the target function is poorly reproduced by the averaging kernel and it is meaningless to carry out an inversion in such unstable conditions. In practice, we would recommend using mode sets similar to sets 1, 2, or 5, and stress that such mode sets can be observed (and were, for example, for some Kepler targets).
Fig. 4. Averaging kernel considering several sets of modes. Upper panel: averaging kernels over the full radius domain. Lower panel: xoom on the region that contributes to the inversion. 
Mode sets for the inversions in a controlled environment.
The inversion results are shown in Fig. 5 for lower mass stars (1.1 M_{⊙} to 1.2 M_{⊙}) and higher mass stars (1.3 M_{⊙} to 1.5 M_{⊙}). The atmosphere model and undershooting affect surface regions, and have no impact on the central entropy profile. The indicator is therefore insensitive to these quantities. It shows promising results for the abundances, and especially for the opacities and the overshoot. The latter indeed induce more significant structural changes, especially for the higher masses, which can be captured by the inversion. For the models with the AGSS09 abundances and masses of 1.2 M_{⊙}, 1.4 M_{⊙}, and 1.5 M_{⊙}, the inversion was less successful or unsuccessful because the relative ratio differences were too small. An inversion has a resolution limit linked with the averaging error and how well the target function is fitted. If the reference model is too close to the observed model, the error on the averaging kernel dominates the inversion and prevents a meaningful correction. Hence, if the ratios are well fitted, it is not appropriate to carry out an inversion.
Fig. 5. Behaviour of the indicator in a controlled environment. Left panel: results for the lower masses (1.1 M_{⊙} to 1.2 M_{⊙}). Right panel: results for the higher masses (1.3 M_{⊙} to 1.5 M_{⊙}). 
4. Hare and hounds
4.1. Presentation of the hares
In this hare and hounds exercise, we aimed at computing models with different physical ingredients, unknown to the hounds, to reproduce as realistically as possible the presence of additional physical processes in observed targets unaccounted for in the seismic modelling procedure. To avoid biases, the hares were generated independently and only the frequencies, as well as the spectroscopic constraints (effective temperature, surface metallicity, and luminosity), were provided to the modeller of the hounds. After all the inversions were carried out, the modeller was given access to the stellar properties and internal profiles of the hares to draw the conclusions of this hare and hounds exercise.
The properties of the hares are presented in Table 2. We varied the opacity tables using either the OPAL or the OPLIB opacities, the reference solar abundances using either the AGSS09 or the GN93 abundances, and the formalism for convection using either the classical mixinglength theory or the full spectrum of turbulence theory (Canuto et al. 1996). The effects of microscopic and turbulent diffusion were considered, and the effects of radiative accelerations on metals were simulated for one hare by solely taking the diffusion of helium and hydrogen into account. Other ingredients were kept unchanged, such as the equation of state (for which we used the FreeEOS equation of state), and the atmosphere model (for which we used a VALC atmosphere). To generate realistic frequency uncertainties, we searched in the LEGACY sample for stars with similar frequency ranges and adopted their observed uncertainties, namely KIC8394589 (Target 1 and 6), KIC12258514 (Target 2), KIC8228742 (Target 3), KIC6508366 (Target 4), and KIC8938364 (Target 5).
Characteristics of the hares.
4.2. Results of the hounds
To generate the hounds as realistically as possible, we used an advanced modelling strategy that involved first conducting a Markov chain MonteCarlo (MCMC) in a grid, fitting the classical constraints (effective temperature, metallicity, and luminosity) and the individual frequencies. We used the AIMS software (Rendle et al. 2019) to perform the MCMC and two grids that differed in their overshooting value. The global properties of the grids are summarised in Table 3 and the physical ingredients of the grid models are the same as for the observed targets in Sect. 3.2, but with one difference that has no influence; partial ionisation was not taken into account. With our indicator, we study central stellar regions of mainsequence stars, where accounting for partial ionisation (or not) only has a very small impact (Turcotte et al. 1998a,b). If radiative accelerations are included or if the target is close to the end of the main sequence, partial ionisation is then not negligible (Schlattl 2002; Deal et al. 2018). However, inversions are quasimodel independent, especially they do not rely on the choice of the transport processes. It is, therefore, not problematic if the reference model does not include the most detailed physics possible for the hare and hounds exercise. Next, we carried out a mean density inversion on the optimal MCMC model to add the inverted mean density to the set of constraints, with a conservative precision of 0.6%, and we conducted new MCMC fitting the classical constraints, the inverted mean density, and frequency separation ratios (r_{01}, r_{02}, or both). Finally, we performed local minimisations with a LevenbergMarquardt algorithm (see e.g. Roweis 1996) by leaving the overshooting parameter free. We summarised the constraints in Table 4 and included the r_{01} in the constraints of some of the hounds to verify that the correction predicted by the inversion was significantly smaller or negligible, as we expected. This also allows us to estimate the impact of the resolution limit. This advanced modelling strategy provided the optimal model of Kepler93 (which was confirmed by a more elaborated modelling procedure, Bétrisey et al. 2022). For stars with physical ingredients that are too different from the physics of the grids, a series of local minimisations are necessary. This was the case for target 5, for which the MCMC was not successful, and for which we conducted local minimisations by starting with guesses of different masses (models LM, M126, M133, M135, and M140), and with different physical ingredients (models Ov020/GN93 and OPLIB/GN93). In general, if we fit frequency separation ratios, only the r_{02} ratios can be fitted most of the time. Indeed, the fit of both r_{01} and r_{02} ratios tends to be unstable. These ratios are sensitive to different stellar properties and a minimisation including both of them has to deal with too many constraints, preventing a successful convergence into an optimal solution. An inversion based on the r_{01} ratios is therefore relevant at this stage, because it can provide quasimodelindependent constraints that can highlight favoured and forbidden regions in the parameter space, and provide new insight for the modelling and understanding of the star.
Global properties of the grids for the MCMC.
Constraints for the targets 1, 2, 3, 4, and 6.
Two main conditions are required for the inversion to provide a meaningful correction (here, we already assume that the linear regime assumption is fulfilled, as it is typically the case with the usual modelling strategies of solarlike stars). First, the indicator must be sensitive to some structural differences; in our case, the differences in the entropy profile in the central stellar regions. Due to our modelling strategy, some of the hounds presented a quasiidentical entropy profile and were thus identical for our indicator. We therefore discarded the redundant models to avoid biase. In the left panels of Fig. 6, we illustrate the nonredundant structural differences of the hounds of targets 2 and 5. The hare profiles are in orange. The entropy profiles of the hounds differ significantly from each other, the height of the entropy plateau is different, the extent of the convective core is different, and the mixing boundary properties are also different. The corresponding indicator values are displayed in the right panels of Fig. 6 (red squares). We observe that the indicator is indeed sensitive to the structural differences. For example, if the overshooting parameter is changed (e.g. models NuOv000 and NuOv010), the indicator is significantly affected by the resulting structural differences. Second, if the ratios are well fitted, the relative ratios differences are too small. The inversion combines these relative differences to compute the correction and if they are too small, the correction is therefore damped. In an ideal case, this behaviour means that the seismic information is well modelled and that additional information cannot be extracted with an inversion. In other words, if the inversion predicts a tiny correction, it confirms that the seismic information is well described by the reference model. However, inversions are numerical methods, and we must ensure that the correction is not driven by numerical noise. In particular, if the ratios are well fitted, the target function may not be well reproduced because of the limited number of modes available. In that case, a nonnegligible correction is driven by the bad fit of the target function and should not be mistaken as a sign that the physics of the reference model is lacking. In Fig. 7, we show examples of relative ratio differences that can and cannot be used for an inversion. We recall that we are conducting linear inversions, and we therefore discarded relative ratios differences that were larger that 0.5 in absolute values, as these were not compatible with the linear formalism. Such large relative differences occur for highorder modes and are purely a numerical artefact. For these modes, the frequency separation ratios are close to zero and automatically produce a large relative difference. This numerical issue is not limiting because the contribution of highorder modes to the correction predicted by inversion is very small. These modes can therefore be discarded without losing any information.
Fig. 6. Entropy profiles and hare and hounds results of targets 2 and 5. The entropy proxy S_{5/3} is in cgs units . (a) Entropy profiles of target 2. (b) Hares and hounds results of target 2. (c) Entropy profiles of target 5. (d) Hares and hounds results of target 5. 
Fig. 7. Comparison of the relative r_{01} differences for models that can (right) and cannot (left) be used for a meaningful inversion. The relative differences above 0.5 are discarded because they are not appropriate for the linear formalism of the inversion. (a) Model LM1 of target 2, well fitted r_{01} ratios that are not suited for an inversion. (b) Model LM of target 2, relative r_{01} differences suited for an inversion. 
The hares fall into two categories of results. For targets 1, 3, 4, and 6, all the ratios (r_{01} and r_{02}) were well fitted, resulting in relative ratio differences that were too small for the inversion. In this situation, the reference value stays within the 1σ interval of the inverted value, and no new knowledge can be extracted. The most interesting results were obtained for targets 2 and 5, for which it was not possible to simultaneously fit the r_{01} and r_{02} ratios. In this case, the inversion achieved significant improvements.
The entropy profiles and results for the first category of targets are shown in Fig. 8. Target 3 shows the reason for which we stressed that an inversion should not be carried out on well fitted ratios. Theoretically, we would expect a negligible correction in such conditions, as seen for targets 1 and 4. However, some of the hounds of target 3 exhibit nonnegligible corrections, which are a combination of physical structural differences and the resolution limit of the inversion. We observe that the inversion corrects in the right direction, but its magnitude is difficult to trust because it is probably largely dominated by the quality of the fit of the target function. The behaviour of the inversions for target 6 is similar to the results for targets 1 and 4, with one exception; the ratios were not well fitted by the LM model and the entropy profiles were significantly different (Fig. 8g). The conditions for a meaningful inversion were therefore met and the inversion behaves accordingly with our expectations, as illustrated in Fig. 8h.
Fig. 8. Entropy profiles and hare and hounds results of targets 1, 3, 4, and 6. The entropy proxy S_{5/3} is in cgs units . (a) Entropy profiles of target 1. (b) Hare and hounds results for target 1. (c) Entropy profiles of target 3. (d) Hare and hounds results for target 3. (e) Entropy profiles of target 4. (f) Hare and hounds results for target 4. (g) Entropy profiles of target 6. (h) Hare and hounds results for target 6. 
The results for targets 2 and 5 are displayed in Figs. 6b and 6d, respectively. Because the relative ratio differences were significant enough to carry out an inversion, we conducted more intensive tests. For target 2, we tested the impact of the tradeoff parameter θ, and our choice of θ = 10^{−4} appears to be a good compromise. Indeed, the inversion result is slightly affected only if a much higher value is considered. In such cases, the target function is less well reproduced because it is equivalent to assuming high observational uncertainties. We also tested two uncertainty sets, where we considered the standard observational uncertainties and checked that equalweighted uncertainties produce similar results. In this fashion, we could discard the possibility of the presence of highly nonlinear modes. For target 5, we tested two mode sets by removing the lowest order mode. For both sets, the inversion significantly improved the reference values, which shows that it is not necessary to use excessively low order modes that could be difficult to observe. Moreover, model Ov020/GN93 strengthens the consistency of the inversion, because it lies at the opposite side of the exact value with respect to the other models, and therefore corrects in the opposite direction. We point out that, even if we do not have the same convection theory in the hounds of target 5, which strongly impacts the stellar structure, the inversion is still able to correct significantly and in the right direction.
In Table 5, we tried to estimate the improvement achieved by the inversion. We computed the average value^{1} of the indicator and defined an accuracy measure as the relative difference of the reference or inverted value with respect to the exact value. For target 2, the accuracy is improved from 5.3% to 2.2%, and for target 5, from 8.6% to 5.4%. This result should be put in perspective with the precision achieved by the inversion. From a pure statistical point of view, the precision is on the order of ∼0.1% to 0.2%, which is negligible in comparison with the systematic uncertainty due to the choice of the physical ingredients, and the resolution limit of the inversion. In practice, a modelling strategy similar to that of Bétrisey et al. (2022) should be considered. An extensive set of models with various physical ingredients should be constructed and the precision measured with the standard deviation of the inversion results. With this approach, the precision is significantly lower than the statistical precision, but it is generally smaller or equal to the precision of the reference models. Hence, the inversion can shift and in some cases shrink the reference range towards the exact value.
Results for targets 2 and 5.
5. Conclusions
In this work, we adapted the SOLA method for frequency separation ratios in Sect. 2 and introduced a new indicator probing the properties of stellar convective cores in Sect. 3.1. Thanks to this approach, the inversion is not affected by surface effects. In Sect. 3.2, we verified our technique in a controlled environment, where the stellar mass and radius were identical in both the reference and observed models. Then, we conducted an intensive hare and hounds exercise in Sect. 4.
Our inversion technique showed promising results for probing the convective core of mainsequence stars, and is especially suited for Ftype stars. Indeed, they have an adequate mass span for the indicator, and the inversion is not sensitive to the surface effects that are difficult to treat for this type of star. Because the entropy differences are localised in the very central stellar regions that are difficult to probe, even with the ratios, we found that a modelling strategy similar to that of Bétrisey et al. (2022) should be considered in a practical application of the method. An extensive set of models should be constructed, and the properties of the core can be inferred from the global behaviour of the inversion on the models of the set. The ratios are influenced by the overshooting, the radiative opacities, and the chemical composition. It is, therefore, not possible to directly constrain one of these quantities, but favoured and forbidden regions can be highlighted in the parameter space. We stress that if the ratios are well reproduced, an inversion is meaningless due to the resolution limit of the method.
In a broader modelling context, the inversion provides quasimodelindependent constraints on core properties that can help to point out the limitations of the reference models and, for example, be the sign of an incorrect or missing physical process. The inversion can also help to exclude some of the reference models, thus improving the accuracy and precision of stellar fundamental parameters such as mass, radius, and age, a crucial objective for the PLATO mission.
Acknowledgments
J.B and G.B. acknowledge funding from the SNF AMBIZIONE grant No 185805 (Seismic inversions and modelling of transport processes in stars).
References
 Adelberger, E. G., García, A., Robertson, R. G. H., et al. 2011, Rev. Mod. Phys., 83, 195 [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Backus, G., & Gilbert, F. 1968, Geophys. J., 16, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Backus, G., & Gilbert, F. 1970, Philos. Trans. R. Soc. London Ser. A, 266, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2009, in Transiting Planets, eds. F. Pont, D. Sasselov, & M. J. Holman, IAU Symp., 253, 71 [NASA ADS] [Google Scholar]
 Bellinger, E. P., Basu, S., Hekker, S., & Ball, W. H. 2017, ApJ, 851, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Bellinger, E. P., Basu, S., Hekker, S., & ChristensenDalsgaard, J. 2019, ApJ, 885, 143 [Google Scholar]
 Bellinger, E. P., Basu, S., Hekker, S., ChristensenDalsgaard, J., & Ball, W. H. 2021, ApJ, 915, 100 [NASA ADS] [CrossRef] [Google Scholar]
 Bétrisey, J., Pezzotti, C., Buldgen, G., et al. 2022, A&A, 659, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [Google Scholar]
 Buldgen, G., Reese, D. R., Dupret, M. A., & Samadi, R. 2015a, A&A, 574, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Reese, D. R., & Dupret, M. A. 2015b, A&A, 583, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Reese, D. R., & Dupret, M. A. 2016a, A&A, 585, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Salmon, S. J. A. J., Reese, D. R., & Dupret, M. A. 2016b, A&A, 596, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Reese, D., & Dupret, M.A. 2017a, Eur. Phys. J. Web Conf., 160, 03005 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Reese, D. R., & Dupret, M. A. 2017b, A&A, 598, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Salmon, S. J. A. J., Noels, A., et al. 2017c, A&A, 607, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Reese, D. R., & Dupret, M. A. 2018, A&A, 609, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buldgen, G., Salmon, S., & Noels, A. 2019a, Front. Astron. Space Sci., 6, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Buldgen, G., Rendle, B., Sonoi, T., et al. 2019b, MNRAS, 482, 2305 [Google Scholar]
 Buldgen, G., Farnir, M., Pezzotti, C., et al. 2019c, A&A, 630, A126 [EDP Sciences] [Google Scholar]
 Buldgen, G., Farnir, M., Eggenberger, P., et al. 2022, A&A, 661, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1964, ApJ, 139, 664 [Google Scholar]
 Chandrasekhar, S., & Lebovitz, N. R. 1964, ApJ, 140, 1517 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2021, Liv. Rev. Sol. Phys., 18, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Clement, M. J. 1964, ApJ, 140, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Colgan, J., Kilcrease, D. P., Magee, N. H., et al. 2016, ApJ, 817, 116 [Google Scholar]
 Cox, J. P., & Giuli, R. T. 1968, Principles of Stellar Structure (New York: Gordon and Breach) [Google Scholar]
 Deal, M., Alecian, G., Lebreton, Y., et al. 2018, A&A, 618, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Brandão, I., Silva Aguirre, V., et al. 2016, A&A, 589, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dziembowski, W. A., Pamyatnykh, A. A., & Sienkiewicz, R. 1990, MNRAS, 244, 542 [NASA ADS] [Google Scholar]
 Eddington, A. S. 1959, The Internal Constitution of the Stars (New York: Dover Publication) [Google Scholar]
 Elliott, J. R. 1996, MNRAS, 280, 1244 [NASA ADS] [CrossRef] [Google Scholar]
 Farnir, M., Dupret, M. A., Salmon, S. J. A. J., Noels, A., & Buldgen, G. 2019, A&A, 622, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse, 15 [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Irwin, A. W. 2012, FreeEOS: Equation of State for Stellar Interiors Calculations [Google Scholar]
 Kosovichev, A. G. 2011, in Lecture Notes in Physics, eds. J. P. Rozelot, & C. Neiner (Berlin: Springer Verlag), 832, 3 [NASA ADS] [CrossRef] [Google Scholar]
 LyndenBell, D., & Ostriker, J. P. 1967, MNRAS, 136, 293 [Google Scholar]
 Otí Floranes, H., ChristensenDalsgaard, J., & Thompson, M. J. 2005, MNRAS, 356, 671 [CrossRef] [Google Scholar]
 Paquette, C., Pelletier, C., Fontaine, G., & Michaud, G. 1986, ApJS, 61, 177 [Google Scholar]
 Pijpers, F. P., & Thompson, M. J. 1994, A&A, 281, 231 [NASA ADS] [Google Scholar]
 Potekhin, A. Y., Baiko, D. A., Haensel, P., & Yakovlev, D. G. 1999, A&A, 346, 345 [NASA ADS] [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [Google Scholar]
 Reese, D. R., Marques, J. P., Goupil, M. J., Thompson, M. J., & Deheuvels, S. 2012, A&A, 539, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rendle, B. M., Buldgen, G., Miglio, A., et al. 2019, MNRAS, 484, 771 [Google Scholar]
 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1 [Google Scholar]
 Roweis, S. 1996, Levenbergmarquardt optimization, Notes, University Of Toronto [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schlattl, H. 2002, A&A, 395, 85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Scuflaire, R., Théado, S., Montalbán, J., et al. 2008a, AP&SS, 316, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Scuflaire, R., Montalbán, J., Théado, S., et al. 2008b, AP&SS, 316, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Sonoi, T., Samadi, R., Belkacem, K., et al. 2015, A&A, 583, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828 [Google Scholar]
 Turcotte, S., Richer, J., Michaud, G., Iglesias, C. A., & Rogers, F. J. 1998a, ApJ, 504, 539 [Google Scholar]
 Turcotte, S., Richer, J., & Michaud, G. 1998b, ApJ, 504, 559 [Google Scholar]
 Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, APJS, 45, 635 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: The r_{10} and r_{02} kernels
For the r_{02} ratios, we have:
where:
The kernel is found by switching a and b in the previous expression.
For the r_{10} ratios, we have:
where:
Again, the kernel is found by switching a and b in the previous expression.
Appendix B: Correlation matrix for the ratios
Let X_{i} be N random variables, and be their respective mean. We consider the random variable Z which is the general form of a frequency ratio:
In order to derive the error propagation, we use the first order Taylor expansion around the mean,
In the following, we introduce the short notation . At the first order, the correlation matrix is equal to
where we used the linearity properties of the covariance matrix, and that the covariance of a random variable with a constant is null.
For a set of N_{f} frequencies ν_{i}, whose uncertainties are described by the covariance matrix E_{ij}, and i being a short notation for the identification pair (n, l), a frequency ratio can be expressed as:
For example, for the r_{02}(n) ratio, the coefficients are:
where the index indicates the identification pair of the corresponding frequency.
The covariance between two ratios is given by
If the frequencies are uncorrelated, which is the typical assumption, one has , where σ_{i} is the uncertainty of the frequency ν_{i}. In this case, the covariance simplifies to
and the 1σ error bar on r is
In the SOLA method, we use relative differences of ratios,
where obs stands for observed and mod for modelled. Eq. (B.11) therefore becomes
The error term in the SOLA cost function is then
where , , is the square of the Frobenius norm of the covariance matrix of the ratios, and N_{r} is the number of observed ratios.
All Tables
All Figures
Fig. 1. Frequency kernels versus ratio kernels of the model 1.3MoAGSS from Sect. 3.1. For both panels, we show the kernel amplitude as a function of the position, and for the structural variable pair (c^{2}, ρ). Upper panel: loworder (l = 0, n = 11) and highorder (l = 0, n = 26) kernels for the individual frequencies, respectively in orange and blue. Lower panel: loworder (n = 11) and highorder (n = 26) kernels for the frequency separation ratios, respectively in orange and blue. 

In the text 
Fig. 2. Target function and entropy profiles for the calibration model. The observed and modelled entropy profiles are in blue and orange, respectively. The target function is in green. 

In the text 
Fig. 3. Impact of the surface effects on the frequencies. The frequencies before (orange) and after (purple) the addition of the surface effects are shown for different harmonic degrees l. 

In the text 
Fig. 4. Averaging kernel considering several sets of modes. Upper panel: averaging kernels over the full radius domain. Lower panel: xoom on the region that contributes to the inversion. 

In the text 
Fig. 5. Behaviour of the indicator in a controlled environment. Left panel: results for the lower masses (1.1 M_{⊙} to 1.2 M_{⊙}). Right panel: results for the higher masses (1.3 M_{⊙} to 1.5 M_{⊙}). 

In the text 
Fig. 6. Entropy profiles and hare and hounds results of targets 2 and 5. The entropy proxy S_{5/3} is in cgs units . (a) Entropy profiles of target 2. (b) Hares and hounds results of target 2. (c) Entropy profiles of target 5. (d) Hares and hounds results of target 5. 

In the text 
Fig. 7. Comparison of the relative r_{01} differences for models that can (right) and cannot (left) be used for a meaningful inversion. The relative differences above 0.5 are discarded because they are not appropriate for the linear formalism of the inversion. (a) Model LM1 of target 2, well fitted r_{01} ratios that are not suited for an inversion. (b) Model LM of target 2, relative r_{01} differences suited for an inversion. 

In the text 
Fig. 8. Entropy profiles and hare and hounds results of targets 1, 3, 4, and 6. The entropy proxy S_{5/3} is in cgs units . (a) Entropy profiles of target 1. (b) Hare and hounds results for target 1. (c) Entropy profiles of target 3. (d) Hare and hounds results for target 3. (e) Entropy profiles of target 4. (f) Hare and hounds results for target 4. (g) Entropy profiles of target 6. (h) Hare and hounds results for target 6. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.