© ESO, 2018

1. Introduction

Inversion techniques have been used for several decades in helioseismology to analyze the structure of the Sun. Amongst the greatest successes of this field, one finds the determination of the base of the solar convective envelope (Kosovichev & Fedorova 1991) and the helium abundance in this region (Antia & Basu 1994a; Kosovichev 1993), as well as the inversion of the sound speed profile (Antia & Basu 1994b). While for the Sun it is possible to determine a whole structural profile with linear techniques, the case of asteroseismic inversions is far more difficult. This is due to the small number of modes and the absence of oscillations of high harmonic degree, which can help to scan through the whole structure of the star. Initial studies have been performed to carry out inversions based on the variational principle for stars other than the Sun (See Gough & Kosovichev 1993a,b; Takata & Montgomery 2002; Basu 2003, for a few examples of such studies.).

Since the advent of space photometry missions such as CoRoT (Baglin et al. 2009) and Kepler (Borucki et al. 2010), we now have seismic data with sufficiently small uncertainties to attempt structural inversions to constrain models of solar-like stars. In the future, the TESS and Plato 2.0 missions (Rauer et al. 2014) will bring further data, continuing what is now named the space photometry revolution.

While for some of the best Kepler targets, full profile inversions could be attempted, they would require very specific conditions. For example, a convective parameter inversion as shown in Kosovichev (1993) would need a very accurate knowledge of the mass and radius of the target as shown in Buldgen et al. (2017b) from the analysis of the linear relations for various structural kernels. Secondly, an efficient way to deal with surface corrections without degrading the capabilities of the inversion technique would need to be found. As a workaround to these difficulties, and to allow more versatile applications of these techniques, Reese et al. (2012) adapted the classical inversion techniques so that they would focus on extracting global information such as the mean density from the oscillation spectra. This global information is defined by integrated quantities which are chosen for their particular ability to probe certain aspects of stellar structure (Buldgen et al. 2015b,a, see for example).

The strength of this approach is that it focuses all of the information provided by the seismic spectra on the determination of one piece of information at a time. Moreover, this information, which is a linear combination of frequency differences is related through the variational integral relations (Dziembowski et al. 1990) to structural characteristics. Our goal is to further correct seismic models which have been built using the classical forward modeling method used in asteroseismology. Once the integrated quantity is defined, one can use the SOLA inversion technique (Pijpers & Thompson 1994) to check whether it is possible to obtain corrections of the chosen structural indicator. The success of the operation depends on whether the integrated quantity behaves linearly and whether the target function can be easily fitted with the amount of seismic information available. Ultimately, one still faces the classical trade-off problem of inversion techniques and the fine-tuning of the parameters has to be done carefully if one wishes to efficiently extract the structural information from the seismic observations.

In this paper, we present results from tests on artificial data for new indicators based on an entropy proxy. We start by presenting the kernels associated with the structural quantity we use in the inversions. Furthermore, we show how this variable naturally reproduces the entropy plateaus in adiabatic convective regions and how its behavior can be used to probe both convective cores and envelopes with custom-built indicators. Probing convective regions and their surrounding layers is crucial as they are likely subject to extra-mixing.

These questions are illustrated by the current uncertainties on the solar tachocline (see Zahn 2007, and references therein) and emphasize the physical complexity linked to convective envelopes. In the Sun, indications of mixing can be seen in the relative sound speed differences or from rotation inversions, where the change from differential to solid body rotation is associated with this particular region of solar structure. In asteroseismology, convective penetration has also been observed in the CoRoT target HD 52265 (Lebreton & Goupil 2012). From a seismic point of view, glitch fitting techniques may help to position the acoustic depth of the convective envelope (see Monteiro et al. 1994, 2000; Verma et al. 2017). However, these techniques also require very high data quality and currently, the signal related to the base of the convective zone is considered by some authors to be too weak to be fully exploited (Verma et al. 2017). In that sense, our approach takes a different path, by focusing on global information rather than localized signatures to provide additional constraints.

Besides additional mixing processes, inaccuracies in the physical ingredients and numerical techniques used in stellar models can leave their mark on the boundaries of convective regions and therefore on the whole stratification of the model (Gabriel et al. 2014). This, in turn, will adversely affect age determinations, and to a lesser degree, the accuracy with which other stellar parameters are obtained. For example, the inaccuracy with which the extent of a convective core is determined induces uncertainties on the age of an observed target that far exceed 10%, emphasizing the importance of constraining convective regions in the current context of the Plato 2.0 mission.

In Sect. 4, we test the accuracy of both indicators for various targets and analyze their error contributions to see whether the inversion can be computed for the best solar-like targets at hand. Finally, we conclude by summarizing our results and comment on further analyses which have to be carried out to fully assess the potential of structural inversions in asteroseismology.

2. Kernels for the entropy proxy

Obtaining a new indicator with asteroseismic inversions means finding a new way to efficiently extract information from the frequencies. Due to the small number of observed modes, the target functions associated with the indicators must be constructed so that they can easily be fitted by structural kernels. Moreover, the choice of the structural variable must be physically motivated. For example, the use of the squared isothermal sound speed u in Buldgen et al. (2015a) to analyze deep regions is motivated by the approximate relation in the core. In Dziembowski et al. (1990), its use was motivated by the problem of measuring the helium abundance in the solar convective envelope.

In this study, we focus on convective regions. Specifically, we wish to be able to analyze the uncertainties linked to the detection of convective cores and extra-mixing at the boundaries of convective regions. However, we have to keep in mind the intrinsic limitations of inversions in asteroseismology. From our previous studies (Buldgen et al. 2017b), we know that kernels like the (A,Γ₁) or (A,Y) kernels, where , cannot be used without a very accurate and precise determination of the radius, which far exceeds observational uncertainties which are of the order of a percent in the case of interferometry. This means that we have to find a new structural pair which efficiently probes convective regions. In this section, we derive new kernels associated with the variable , with P the pressure and ρ the density, and justify their choice as an efficient probe of convective regions.

The choice of S_5/3 as a structural variable stems from its relation with the entropy of an ideal gas. Using thermodynamical relations (see e.g. Kippenhahn et al. 2012, for thermodynamical relations from which the following equation can be derive), one can show that (1)with k_B being the Boltzmann constant, μ the mean molecular weight and m_u the atomic mass unit. In this equation, f only depends on the mean molecular weight and various physical constants.

This quantity has the interesting property of forming a plateau in the adiabatic convection zones. The height of this plateau is related to the temperature and mean molecular weight gradients in the vicinity of the convective zone’s boundary and thus to the stratification of these poorly modeled regions. We illustrate this property in Sect. 3, when we derive the seismic indicators.

Due to the limited number of frequencies in asteroseismology, the second variable of the structural pair must be chosen so that the cross-term is naturally small. Two variables satisfy this condition: Γ₁, the adiabatic exponent defined as and Y, the helium abundance. This means that the structural pairs we are aiming for are the (S_5/3,Γ₁) and (S_5/3,Y) pairs, which can be derived from the (ρ,Γ₁) and (ρ,Y) pairs. One could also use the (u,Γ₁) and (u,Y) pairs to derive the differential equations without any further difficulties. However, since in practice the (u,Y) and (u,Γ₁) kernels are already obtained from the numerical resolution of a second order differential equation, it is wiser to use the (ρ,Γ₁) and (ρ,Y) kernels as a starting point to avoid multiplying the sources of numerical errors.

Using the direct method presented in Buldgen et al. (2017b), the equation obtained for the (S_5/3,Γ₁) kernels is a second order differential equation written as follows (2)with where r is the radial position, R the total radius of the star, and m the mass contained within the sphere of radius r, , . As was the case for the (u,Y) kernels, one can use exactly the same equation to obtain the (S_5/3,Y) kernels from the (ρ,Y) kernels by simply taking and .

The and are directly obtained from the following algebraic relations: Examples of such kernels are illustrated in Fig. 1.

Fig. 1 Structural kernels for the n = 15,ℓ = 0 mode associated with the (S5/3,Γ1) structural pair on the left-hand side and with the (S5/3,Y) pair on the right-hand side for Target 1 of Table 3. The kernels are presented in their non-dimensional form.

The integral equation used as a starting point for the problem is used to derive the additional boundary condition required to solve the second order equation. This equation is written (5)and the additional boundary condition imposed on the new structural kernels is that they satisfy this equation at one point. This can be done either iteratively or by decomposing the problem into a homogeneous component and a non-homogeneous component (see Buldgen et al. 2017b, for further details).

Once the kernels are derived, Eq. (5)also provides a first verification step to ensure that the kernels are in agreement with the initial steps of their derivation. The second verification is to ensure that the variational integral expressions are satisfied by the new kernels. In other words, ensure that we have: (6)Both verifications are illustrated in Fig. 2 for the (S_5/3,Γ₁) pair and Fig. 3 for the (S_5/3,Y) pair. The order of magnitude of the agreement is very similar to what is found for classical kernels such as the (ρ,c²) and the (ρ,Γ₁) structural pairs. However, small differences can always be seen when changing the structural pair, as presented in Buldgen et al. (2017b). We recall that to ensure the verification of the variational equations, the scaling method mentioned in Buldgen et al. (2017b) is of course necessary since the observed target and the reference model may not have the same radius.

Fig. 2 Left panel: verification of Eq. (5)for Target 1 of Table 3 for the n = 16,ℓ = 0 mode kernel , where is the right-hand side of this equation and ℛ is the residual. All quantities are presented in their non-dimensional form. Right panel: verification of Eq. (6)for modes of degree ℓ = 0,1,2,3 and various radial orders between Target 1 of Table 3 and a model of the same evolutionary sequence, 500My younger.

Fig. 3 Left panel: verification of Eq. (5)for the (S5/3,Y) structural pair for n = 16,ℓ = 0 mode kernel , where is the right-hand side of this Equation and ℛ is the residual. All quantities are presented in their non-dimensional form. Right panel: verification of Eq. (6)for the (S5/3,Y) structural pair for modes of degree ℓ = 0,1,2,3 and various radial orders. The models used are the same as in Fig. 2.

3. Using the entropy proxy to obtain indicators of convective regions

In this section, we illustrate the sensitivity of the entropy proxy, S_5/3, to convective regions and we show how to build seismic indicators useful for asteroseismic targets using the newly derived kernels.

Inverting an integrated quantity using the SOLA inversion technique actually consists in finding the linear combination of frequency differences which best reproduces the correction to be applied to this integrated quantity. As noted before, the indicator is optimized to probe certain parts of the stellar structure. We first present the generic case of a quantity A, the linear perturbations of which are given by (7)with being the target function used for the inversion and s₁ a variable such as ρ, c² or S_5/3. We reiterate that a generic expression of the linear integral relations can be written. (8)with s₂ being the second variable of the integral relations and the being the kernel functions related to each structural variable.

The inversion process estimates the differences between the observed integrated quantity, A_Obs and that of the reference model, A_Ref, using the following relation (9)with δν = ν_Obs−ν_Ref and the c_i being the inversion coefficients determined from the minimization of the SOLA cost function (see Pijpers & Thompson 1994, and Eq. (10)). The inverted difference in Eq. (9)is not exactly the real difference in Eq. (7), due to the intrinsic limitations of the inversion.

The inversion process as a whole is a trade-off between various antagonistic terms. This is understood from the cost function used to carry out the SOLA inversion (10)which contains four terms. The first term is associated with fitting the target function: it ensures the accuracy of the inversion and includes the averaging kernel, defined as (11)with N the total number of individual observed frequencies. The second integral deals with the cross-term kernel, defined as (12)which is a pollution term by the variable s₂ whose contribution has to be damped. The third term is linked with the propagation of observational error bars of the individual frequencies, denoted here σ_i. The last term is linked to the additional normalization constraint on parameter k derived using homologous relations (see Sect. 3.3 and (Reese et al. 2012) for a description of this approach) and unlike β and θ, η is not a free parameter of the inversion but a Lagrange multiplier. The free parameters of the inversion, β and θ, are used to give more or less importance to each of the antagonistic terms and are thus called the trade-off parameters.

The optimal set of parameters is defined by analyzing the amplitude of the first three individual terms of Eq. (10)and visually verifying that the fit is reliable. In practice, this implies that the inversion will be a compromise between precision, accuracy, and cross-term. As stated in Sect. 2, the cross-term is damped by a suitable choice of variables which has been validated in previous studies (Basu 2003; Reese et al. 2012). Consequently, the main part of the trade-off problem is to find a suitable compromise between precision and accuracy. This is done by comparing the terms responsible for these characteristics in so-called trade-off curves (see Backus & Gilbert 1967; Pijpers & Thompson 1994; Rabello-Soares et al. 1999; Reese et al. 2012, for the full definition of this concept and various applications), where one plots the following quantities for various values of θ. In practice, a visual inspection of the agreement of the averaging kernels with their target is also informative on the quality and reliability of the inverted results. In practice, the set of inversion parameters will depend on the target function of the inversion and the observed seismic data.

In addition to these reliability assessments, further analyses can be performed, by separating the contribution to the linear correction derived from the frequencies using Eq. (9). Indeed, inversion techniques are subject to multiple error sources which can sometimes damp each other. To analyze this potential compensation, the linear correction of a given integrated quantity A can be decomposed into various contributions. To do this, we introduce the quantities A_Ref, A_Obs and A_Inv; the reference, observed and inverted values of the indicator, respectively. For these quantities, one has the following relations By computing the difference between the inverted and observed values of the indicator, one can isolate the errors in individual contributions, such that one as (17)where we define three main error contributions: ϵ_Avg, the error stemming from the mismatch between the averaging kernel and its target, which depends on the quality of the dataset and the value of the θ parameter; ϵ_Cross, the error stemming from the non-zero cross-term contribution which depends on the choice of the variable s₂ and the parameter β; and ϵ_Res, the residual error, which is defined as (18)This contribution is the most difficult to assess, since it can originate from surface effects, the linearization of the equation of state, or from non-linear effects. In inversions of observed targets, the residual errors would also be influenced by physical processes not included in the derivation of the variational relations and on systematics in the frequency determinations. However, since its calculation requires the knowledge of the structural differences between the target and the reference model, the value of this residual error is not accessible in practical cases. As we see in Sect. 4, these error contributions are used in combination with the classical trade-off analysis to determine the degree of reliability of the inversion and its accuracy.

3.1. Convective cores and deep regions

Convective cores represent one of the major difficulties when studying the evolution of stars with masses higher than approximately 1.2 M_⊙. Indeed, their presence can lead to large uncertainties in age determinations and can completely change the evolutionary track of a given model. Various studies focus on the uncertainties linked to convective cores; for example the recent studies by Deheuvels et al. (2016) and Claret & Torres (2016) to calibrate overshooting, or the derivation of dedicated seismic indices to the detection of a convective core in a given star (see Miglio & Montalbán 2005).

Looking at structural profiles of a model, the presence of a convective core can easily be seen in derivatives, where it introduces a discontinuity. This means that indicators based on derivatives, such as the t_u indicator presented in Buldgen et al. (2015a), defined as (19)or t presented in Buldgen et al. (2015b), defined as (20)will be extremely sensitive to convective cores. However, the target function associated with these indicators is impossible to fit with the structural kernels if a convective core is present in the model. Indeed, both of the target functions for t and t_u are strongly discontinuous whereas the structural kernels are not. Using S_5/3 as the structural variable for the kernels, we recover the sensitivity to convective cores without the need for a derivative. Moreover, as can be seen in Fig. 1, the structural kernels associated with S_5/3 have an increased intensity in the deep regions, around 0.2r/R, due to the mass dependency of the entropy proxy.

We define a new indicator for convective cores as follows: (21)with f(r) being the weight function associated with the indicator: (22)The target function is built to fit the lobe around 0.2r/R present in every kernel of either the (S,Γ₁) or the (S,Y) structural pair. It has seven non-dimensional free parameters that can be varied if required. We give an example of values for these parameters in Table 1. Their values will depend on the stellar type, because variations of the reference model of a given star will imply a slightly different behavior of the structural kernels. For example, more massive or evolved stars will tend to present a steeper decrease of the lobe around 0.2r/R, inducing a slight increase of the parameters a₂ and a₅. The asymmetry of the lobe in the kernels can also be fitted by varying the parameters a₁ and a₄, which is also a function of the stellar type. Moreover, the parameters a₃ and a₆ might be reduced because the maximum of the lobe will be located towards deeper regions.

The linear perturbation of this indicator leads to the following target function for the inversion: (23)This means that the target function of this inversion is defined as (24)This function is illustrated in blue in Fig. 4 along with the profile of 1 /S_5/3 in red. The product of both curves gives the argument of the integral defining the indicator. From visual inspection, it is clear that this argument probes the inner layers of the acoustic structure of the star. One should also note that this indicator is not restricted to stars with convective cores but can also be a complement to or a replacement for the t_u indicator.

Fig. 4 Target function of the core indicator SCore plotted alongside the structural function . The model used is Target 3 of Table 3. The target function is plotted in its non-dimensional form.

If the convective core is well established and has quite a high plateau with a steep entropy variation, a peak tends to appear in the target function, as illustrated in Fig. 5. This peak is in fact due to the boundary of the convective core. Due to the behavior of the structural kernels in the center, the entropy plateau of the convective core is erased, but if its height and radial extent are sufficient, some traces remain in the kernels. Consequently the target function defined in Eq. (21)can still be fitted and is well adapted to extract information about convective cores. From numerical tests, the linear relations between frequencies and structure (Eq. (2)) are still satisfied even in these particular cases. We see in Sect. 4 how efficient the inversion actually is when confronted with these effects. The weight functions can also be adapted to be more easily fitted if required. Ultimately, the diagnostic power is limited by the detected modes and their error bars.

Fig. 5 Left panel: structural profile of in red, showing the plateau in the convective core. In blue, an example of the adaptation of a target function to include the trace of the convective core and its border. The target function is plotted in its non-dimensional form. Right panel: structural kernels associated with S5/3 in the (S5/3,Γ1) structural pair, showing the trace of the entropy plateau in the convective core in the central regions. The model used is Target 6 of Table 3. The kernels are plotted in their non-dimensional form.

3.2. Convective envelopes

In this section, we show that it is possible to probe the upper radiative zone and the convective envelope of certain Kepler targets. However, fully isolating the convective envelope is done thanks to high ℓ modes, as in helioseismology. Since these modes are not accessible to asteroseismology a complete scan of the structure is not achievable. However, a few ℓ = 3 modes can still help with extracting information from the observations. The efficiency of this technique thus depends on the ability to detect octupole modes. However, such modes are difficult to observe, even for Kepler targets, and their frequencies are often determined with a rather low precision.

Our objective here is to probe the plateau in the convective envelope of S_5/3, shown in Fig. 6, and layers immediately below it. However, a compromise has to be found to exclude the uppermost region of the star, subject to surface uncertainties and a breakdown of the assumptions behind inversion techniques. Provided that the data quality is sufficient, the following indicator can be fitted (25)which means that the target function for this inversion will be (26)stemming from the relative linear perturbation of Eq. (25) and not forgetting the constant factor S_env in the denominator. The weight function g(r) is somewhat complicated and has been built to match the behavior of kernels while trying to extract information in the upper regions. To explain the choice of g(r), we separate it into three components: (27)

Fig. 6 Target function of the envelope indicator SEnv plotted alongside the structural function S5/3 for Target 1 of Table 3. The target function is plotted in its non-dimensional form.

with the following definitions: The function g₁ is used to probe deeper regions where the entropy is influenced by the way it is reconnected to the plateau of the convective envelope. Looking at Fig. 6, we see that this corresponds to the regions where a slope starts in the entropy profile. The g₂ component is a Fermi-Dirac distribution that reproduces the entropy plateau in the convective envelope. The g₃ function eliminates the surface regions as efficiently as possible through the hyperbolic tangent. This component is steep in order to avoid affecting the lower part of the profile. The target function combining all three components is represented in Fig. 6. These components can be further adjusted depending on the reference model and the observed modes. The total target function of the S_Env indicator thus counts 12 non-dimensional free parameters. We illustrate some examples of values for these parameters in Table 2. All of these parameters can vary depending on the stellar type. For instance, an increase in stellar mass will reduce the size of the convective envelope and thus lead to a more difficult fit of the target function of the S_Env indicator. This implies that the parameter b⁸ to b₁₁ will be influenced since the plateau of S_5/3 will be located at a different place. This also implies that the g₁ function has to be adapted. The peak has to be made larger and more asymmetrical with a less steep transition towards the convective envelope position, where the g₂ sees its amplitude increase. In practice, the parameters that vary the most are the parameters b₂ to b₇. It should be noted that the peak cannot be simply shifted, because of the low degree of the modes, the maximum of the target function will always be located around 0.4r/R which is the highest point in normalized radius where a form of localization of the kernels can be achieved.

For example, we see in Sect. 4, that octupole modes are required to carry out inversions of the S_Env indicator. The number of these modes and their error bars also influence the paramaters of the g₁ function. A larger number of precisely determined octupole modes allows for a higher position of the peak in the g₁ function and a more slowly decreasing slope towards the convective envelope, implying a better analysis of these regions. Besides the observed modes, the reference model also affects the building of the target function of the S_Env indicator. For example, a more massive star, having a shallower convective envelope, will require an adaptation of the g₂ and g₁ functions to avoid building a target that cannot be easily fitted with only low ℓ modes.

3.3. Non-linear generalization

As in Buldgen et al. (2015a) and Reese et al. (2012), one can improve the inversion by using a simple non-linear generalization of the inversion procedure. To do this, we follow the procedure in Sect. 3.2 of Buldgen et al. (2015a). From homology scaling, we can show that (31)This homology scaling is also crucial to know how the indicator value should be rescaled to be comparable for various reference models. We have seen in Buldgen et al. (2015a) that the inversion implicitly rescales the target model to the radius of the reference model without changing its mean density. This means that the quantity S_5/3 is then also rescaled. Consequently, when comparing inversions from various reference models, one should always be aware of this scaling and present them in a form proportional to the mean density. For example, using Eqs. (21)and (22), one can carry out a homology scaling of the core conditions indicator, leading to (32)meaning that the quantity GS_Core scales as ρ^{− 1/3}. For S_Env, homology scaling from Eqs. (25)and (30)leads to (33)Meaning that this quantity needs to be rescaled if it is to be compared for various reference models. In other words, one has S_Env/GR^4.5 ∝ ρ^1/3. To provide the coefficients for the non-linear generalization, we focus on the mass dependency of the quantities, since the radius is implicitly kept constant by the inversion. Looking at Eq. (31), this means that S_5/3 ∝ ν^2/3 in terms of the M dependency, thus leading to (34)for a homologous transformation which keeps the radius between models constant but changes their mean density through their mass. As stated in Sect. 3.2 of Buldgen et al. (2015a), the coefficient in this linear relation between perturbations of the frequencies and that of a structural quantity, denoted k in our previous paper, is crucial to derive the non-linear generalization. In this study, one has k = 2/3 in the relation between S_5/3 and the frequencies. Looking at Eq. (25), which defines S_env, we can see that since this indicator is proportional to S_5/3. Hence, one will have (35)and the value of 2/3 can be applied to Eqs. (31) and (34) of Buldgen et al. (2015a) defining the optimal value of S_Env and its associated errors bars in the framework of this non-linear generalization. Similarly, the value k will also be fixed to 2/3 in the additional condition on the inversion coefficients (fourth term of Eq. (10)) used to improve the regularization of the inversion process.

Using the same analysis, it is easy to show that the S_Core indicator will satisfy a relation of the opposite sign, due to its opposite mass dependency, (36)since it is proportional to . Consequently, the value −2/3 can also be used in Eqs. (31), (34) and (28) of Buldgen et al. (2015a) to derive the non-linear generalization of this indicator.

3.4. Relation between S_core, S_env, and stellar structure

In the preceding sections, we have shown how the entropy proxy could be used to obtain indicators of both convective cores and envelopes. In this section we will briefly show how they change with some specific aspects of stellar structure. We mention that these changes of course depend on the parameters used to build the target functions of the indicators. For example, placing the maximum of the Gaussian functions deeper in the S_Core target function will increase the changes due to convective cores.

We illustrate the relation between convective cores and the S_core indicator in Fig. 7 by plotting the rescaled profiles of two 1.3 M_⊙, 1.6 Gy models with a solar chemical composition. The model plotted in blue in Fig. 7, included a 0.15 pressure scale-height adiabatic overshoot, causing its convective core to be larger at a given age and inducing changes in its entropy plateau. Similarly, the peak in the S_core target function is much more pronounced in the blue profile than in the red, indicating a more extended convective core. The changes seen in the indicator thanks to this small variation are of around one percent of the total value of the indicator. This is of course quite small but results from the fact that both models have the same age and chemical composition, and only differ in one physical ingredient. In a more realistic case, where the models are selected using seismic constraints, the differences between reference model and target can be much larger.

Fig. 7 Left panel: scaled profile used in inversions for two 1.3 M⊙ models. The blue curve is associated with a model built using a 0.15 pressure scale height adiabatic overshoot. Right panel: target function of the Score indicator for both models.

Similarly, we illustrate in Fig. 8 the effects of opacity changes on solar models. We compare solar models using the AGSS09 and the GN93 abundances and different opacity tables. To further increase the differences, we do not include microscopic diffusion in the AGSS09 solar model. This leads to variations in the position of the entropy plateau of the convective envelope of the models, as can be seen in the left panel of Fig. 8. This emphasizes the direct link between the entropy plateau and the temperature gradient in stellar models. The right panel of this Figure illustrates the small changes in the target function of S_Env. The changes in the target function induce a change in the value of the indicator of less than one percent. Again this is quite small but is mainly due to the fact that the models have the same mass, the same age and are not selected on the basis of their seismic information. Section 4 shows a very different picture for its first target, a solar model including various physical processes fitted using individual small frequency separations without including these additional processes in the reference model.

Fig. 8 Left panel: scaled S5/3 profile used in inversions for two standard solar models. The blue curve is associated with a model built using the AGSS09 abundances while the red one uses the GN93 abundances. Right panel: target function of the Senv indicator for both models.

4. Tests on artificial data

In this section, we present the results of exercises carried out to test the accuracy and the capabilities of the indicators based on the entropy proxy. We test these capabilities in similar exercises as in Buldgen et al. (2015a), using targets built with certain specificities that the reference models of the inversion did not include. For example, some targets included strong overshooting or turbulent diffusion and were fitted with reference models without these processes. This approach attempts to simulate as best as possible the biases due to physical simplifications or approximations in our representation of stellar structure and evolution when carrying out seismic modeling of observed targets. It is also a good test of robustness and accuracy to see whether the linear approximation and the SOLA method are indeed capable of probing such aspects of stellar modeling.

In total, we built ten targets with various masses and ages, but only present six to avoid redundancy. The reference models were built by fitting the targets’ individual small frequency separations and their effective temperatures. The fit was carried out using a Levenberg-Marquardt algorithm for the minimization. The observational spectrum of the targets is given in Table 5. An error bar of 80 K was used on the effective temperature and the observational error bars on the frequencies were similar to those found for the best Kepler targets. The physical characteristics of the target models are given in Table 3. Some parameters have rather extreme values, in order to test the inversion and its limitations. Targets 1 to 5 were built using the OPAS opacities (Mondet et al. 2015; Le Pennec et al. 2015) and the Ceff equation of state (Christensen-Dalsgaard & Daeppen 1992) with the AGSS09 abundances (Asplund et al. 2009) while Target 6 was built using the OPAL opacities (Iglesias & Rogers 1996) and the OPAL equation of state (Rogers & Nayfonov 2002a) with its corrections (Rogers & Nayfonov 2002b). The reference models all used the Ceff equation of state and the OPAL opacities along with the GN93 abundances (Grevesse & Noels 1993). None of the reference models included turbulent diffusion. The formalism used to implement turbulent diffusion in the target models is presented in Miglio et al. (2007). All models were computed using the Liège Stellar Evolution Code (Cles, Scuflaire et al. 2008b) and the oscillations were computed using the Liège Oscillation Code (Losc, Scuflaire et al. 2008a). The reference models obtained from this first step of forward seismic modeling are presented in Table 4. As can be seen, most of the reference models, denoted “Ref”, still present significant structural differences from their respective target from Table 3. The goal of the exercises is now to determine whether the inversions of the structural indicators can detect these mismatches, and thus provide additional information which could lead to refinements of the seismic modeling in the study of an observed target¹.

We also tested the impact of supplementary observed modes on the accuracy of the inversions; for example ℓ = 4 and 5 modes but also modes with higher and lower n. Our goal was to assess which oscillations modes could help with extracting more information on stellar structure using inversions of integrated quantities. We discuss this at the end of Sect. 4.2. The error bars were usually taken to be similar to those of the best Kepler targets, that is around 3 × 10^-1μHz with slightly larger error bars for the lowest and highest frequencies, as expected in observed cases.

4.1. Results for S_Core inversions

In Table 6, we summarize the results of our test cases. Some kernels are illustrated in Fig. 9. In these tests, we consider observed quantities, denoted with the subscript “Tar”, to be those of the target models.

From the inversions for Targets 2 and 5, we can see that the S_Core indicator efficiently detects an inaccuracy in the stratification of the reference models, due to the absence of convective cores in these models. For Target 3, we can see that the inversion provides an accurate result, but that the error bars on the inverted result are too large to conclude that the model has to be rejected. From a modeling point of view, we can see from Tables 4 and 3 that Ref₃ is already a very good fit to Target 3. Adding more frequencies to this test case allowed to eliminate this problem but in practice, this would mean that for this particular case, the S_Core indicator would only be an additional check for the modeling. However, in most test cases, the observed differences between reference and inverted result are smaller than the error bars of the inverted result, meaning that we are still safe in terms of trade-off parameters.

Moreover, it should be noted that the values of the parameters in the definition of these indicators are not fixed. The number of frequencies available and the characteristics of the convective core (if present) imply that the parameters considered optimal in one case might be suboptimal in another case, depending on the mass, evolutionary stage or chemical composition of the model. For each of the targets in Table 6, we fine-tuned the parameters to analyze the diagnostic potential using the oscillation modes of Table 5. Further tests on targets of the Kepler Legacy sample (see Buldgen et al. 2017a) have also shown similar behaviors and proved that the method was indeed applicable to current seismic data.

From Table 6, we can conclude that the inversion is efficient and can provide a diagnostic of inaccuracies in our modeling of deep regions of solar-like stars with convective cores. Very slight compensation is seen for Targets 5 and 6, but further tests have shown these cases to be marginal. In addition, the averaging kernels illustrated in Fig. 9 fit their target function to an acceptable level of accuracy.

Fig. 9 Left panel: averaging kernel for the SCore indicator (blue) and the target function for the inversion (green). Right panel: cross-term kernel for the SCore inversion (blue), and target function which is 0 (green). These results are presented for Target 2 of Table 3.

4.2. Results for S_Env inversions

In Table 7, we present inversion results for the S_Env indicator. Some kernels are shown in Fig 10. We used the same targets, oscillation spectra and naming conventions as for the S_Core indicator.

From Table 7, we can see that reproducing the value of S_Env can be done efficiently for most targets. However, for Target 5, for example, the inversion is very difficult because of the proximity between the target and its reference model. The variation of the indicator is too small to be seen with the typical accuracy of asteroseismic data. In comparison, some standard solar models show larger differences from one to another than the difference seen between Target 5 and its reference model. In addition, Target 5 is more massive, meaning that its envelope does not go as deep as in a lower-mass star. In fact, we found the base of the envelope to be around 0.75 fractional radii in its reference model. Therefore, the sensitivity of the indicator has to be increased by pushing the weight function towards upper regions. However, the structural kernels, illustrated in the right panel of Fig. 5 for a massive star with a convective core (noticeable by the peak in the deeper layers) show the exact opposite trend, the more massive the model, the steeper the decrease towards upper regions. Consequently, the amount of seismic information required to probe the convective envelope of massive stars is higher than for low-mass stars, for which the entropy plateau goes deeper and the kernels are more suited to their purpose. We confirm this by carrying out supplementary test cases for additional low-mass targets and show that these can be probed with the S_Env indicator. We illustrate one of these tests in Fig. 11 and also test the effects of reducing the number of observed frequencies. For the S_Core indicator, reducing to 44 frequencies still works. However, by changing the parameters in the target function of the indicator, valuable information can still be gained with as low as 35 observed frequencies. For the S_Env indicator, using fewer than 45 frequencies may already lead to imprecise results and we noticed that ℓ = 3 modes were required to ensure an acceptable fit of the target function. Its range of application is therefore limited to the very best observed Kepler targets. In this particular case, we used a 0.9 M_⊙ model including microscopic diffusion and fitted the individual large and small frequency separations of this artificial target with models that did not include microscopic diffusion and had a different helium mass fraction than the target.

Including higher ℓ modes can help improve the fit of the target function, but not all radial orders are equivalent in this matter. As already observed for the t_u inversion, low n modes are better at fitting custom-made global quantities, while high n modes are barely used and associated with lower inversion coefficients. Indeed, at very high values of n, one reaches the asymptotic regime and the oscillation spectrum becomes very regular, thus the information to be extracted from the modes is degenerate. Moreover, the large error bars on these high frequencies might make them very difficult to use for precise determinations of seismic indicators and these modes are more affected by the surface effects, biasing their seismic diagnostic.

Fig. 10 Left panel: averaging kernel for the SEnv indicator (blue) and the target function for the inversion (green). Right panel: cross-term kernel for the SEnv inversion (blue), and target function which is 0 (green). These results are presented for Target 1 of Table 3.

Fig. 11 Left panel: illustration of the accuracy of the inversion of the SCore indicator illustrated in a plane for various numbers of observed frequencies. The reference values are given in red, the target values are plotted in green and the inverted values are plotted in various colors according to the number of observed frequencies. Right panel: same for the SEnv indicator.

5. Conclusion

From Sect. 4, we can see that both indicators are well suited to probing convective cores and envelopes. We can see from Table 6 that the residual error remains small compared to the variations of the S_Core indicator caused by the changes of the entropy plateau in the convective core. Consequently, we can conclude that the S_Core indicator is suitable for analyzing mismatches in the deep layers of solar-like stars, even if they present a convective core. The S_Env indicator, on the other hand, is quite efficient in its analysis of upper layers and can provide further constraints even if other indicators, such as the t_u indicator, are fitted within their error bars.

Another aspect of inversions that has to be re-analyzed is linked to the surface effects. While inversions related to core-sensitive structural aspects are naturally less prone to showing problems, this might not be the case for all inversions. Actually, additional tests have shown that some inaccuracies, which appear as an increase of the residual error of the inversion, might be expected. Therefore, one needs to analyze how empirical corrections such as those of Sonoi et al. (2015) and Ball & Gizon (2014) might help solve the problem, since the classical surface-correcting method used in helioseismology cannot be used in asteroseismology due to the limited number of observed frequencies, this analysis will be carried out in future studies.

In conclusion, we have proved that the amount of seismic information found in the typical spectra of solar-like Kepler targets is sufficient to carry out more robust inversions of a core condition indicator, applicable to more massive stars. We have also shown that for the very best of these targets, for which octupole modes are observed, the inversion of a structural indicator probing the regions near the base of the convective envelope could be attempted, providing supplementary insight into the structural properties of the target. However, both S_Core and S_Env might require a slight adaptation of their target functions when applied to observations depending on the dataset and the reference model for the inversion. Nevertheless, they can be used alongside additional interferometric, spectroscopic or even seismic constraints such as, for example, glitch-fitting techniques. In turn, the use of these indicators will help improve the physical accuracy of stellar models and provide stringent constraints on stellar fundamental parameters required by other fields such as exoplanetology and Galactic archeology.

Acknowledgments

G.B. is supported by the FNRS (“Fonds National de la Recherche Scientifique”) through a FRIA (“Fonds pour la Formation à la Recherche dans l’Industrie et l’Agriculture”) doctoral fellowship. This article made use of an adapted version of InversionKit, software developed in the context of the HELAS and SPACEINN networks, funded by the European Commission’s Sixth and Seventh Framework Programmes. We thank the referee for his useful comments which have helped improve the manuscript.

References

All Tables

All Figures

	Fig. 1 Structural kernels for the n = 15,ℓ = 0 mode associated with the (S5/3,Γ1) structural pair on the left-hand side and with the (S5/3,Y) pair on the right-hand side for Target 1 of Table 3. The kernels are presented in their non-dimensional form.
In the text

Fig. 2 Left panel: verification of Eq. (5)for Target 1 of Table 3 for the n = 16,ℓ = 0 mode kernel , where is the right-hand side of this equation and ℛ is the residual. All quantities are presented in their non-dimensional form. Right panel: verification of Eq. (6)for modes of degree ℓ = 0,1,2,3 and various radial orders between Target 1 of Table 3 and a model of the same evolutionary sequence, 500My younger.

In the text

Fig. 3 Left panel: verification of Eq. (5)for the (S5/3,Y) structural pair for n = 16,ℓ = 0 mode kernel , where is the right-hand side of this Equation and ℛ is the residual. All quantities are presented in their non-dimensional form. Right panel: verification of Eq. (6)for the (S5/3,Y) structural pair for modes of degree ℓ = 0,1,2,3 and various radial orders. The models used are the same as in Fig. 2.

In the text

	Fig. 4 Target function of the core indicator SCore plotted alongside the structural function . The model used is Target 3 of Table 3. The target function is plotted in its non-dimensional form.
In the text

Fig. 5 Left panel: structural profile of in red, showing the plateau in the convective core. In blue, an example of the adaptation of a target function to include the trace of the convective core and its border. The target function is plotted in its non-dimensional form. Right panel: structural kernels associated with S5/3 in the (S5/3,Γ1) structural pair, showing the trace of the entropy plateau in the convective core in the central regions. The model used is Target 6 of Table 3. The kernels are plotted in their non-dimensional form.

In the text

	Fig. 6 Target function of the envelope indicator SEnv plotted alongside the structural function S5/3 for Target 1 of Table 3. The target function is plotted in its non-dimensional form.
In the text

	Fig. 7 Left panel: scaled profile used in inversions for two 1.3 M⊙ models. The blue curve is associated with a model built using a 0.15 pressure scale height adiabatic overshoot. Right panel: target function of the Score indicator for both models.
In the text

	Fig. 8 Left panel: scaled S5/3 profile used in inversions for two standard solar models. The blue curve is associated with a model built using the AGSS09 abundances while the red one uses the GN93 abundances. Right panel: target function of the Senv indicator for both models.
In the text

	Fig. 9 Left panel: averaging kernel for the SCore indicator (blue) and the target function for the inversion (green). Right panel: cross-term kernel for the SCore inversion (blue), and target function which is 0 (green). These results are presented for Target 2 of Table 3.
In the text

	Fig. 10 Left panel: averaging kernel for the SEnv indicator (blue) and the target function for the inversion (green). Right panel: cross-term kernel for the SEnv inversion (blue), and target function which is 0 (green). These results are presented for Target 1 of Table 3.
In the text

	Fig. 11 Left panel: illustration of the accuracy of the inversion of the SCore indicator illustrated in a plane for various numbers of observed frequencies. The reference values are given in red, the target values are plotted in green and the inverted values are plotted in various colors according to the number of observed frequencies. Right panel: same for the SEnv indicator.
In the text