Free Access
Issue
A&A
Volume 607, November 2017
Article Number A58
Number of page(s) 5
Section The Sun
DOI https://doi.org/10.1051/0004-6361/201731354
Published online 13 November 2017

© ESO, 2017

1. Introduction

In the last decades, helioseismology has led to striking successes. The precise location of the base of the convective envelope at a fractional radius of 0.713 ± 0.001 (Kosovichev & Fedorova 1991); the inversion of the solar sound speed, density, and rotation profiles (Antia & Basu 1994; Kosovichev et al. 1997); the determination of the helium mass fraction in the convective envelope at Y = 0.2485 ± 0.0035 (Basu & Antia 1995); and the outcome of the “solar neutrino problem” (Bahcall & Peña-Garay 2004) are amongst the greatest achievements in this field. In the 1990s, the internal structure of the Sun was extremely well reproduced by standard solar models (hereafter SSMs). Therefore, the physical ingredients of these numerical models, particularly the solar chemical element abundances (Grevesse & Noels 1993, hereafter GN93), were applied to stars other than the Sun and used to compute grids of stellar models. These grids are one of the basic components in various fields such as stellar population analysis, Galactic evolution, and exoplanetology.

Later on, the physical ingredients of the solar models, such as the equation of state (Rogers & Nayfonov 2002) or the heavy element abundances, were continuously refined (Grevesse & Sauval 1998, hereafter GS98), but the changes being quite small, the agreement of the models with helioseismology remained. However, two refinements with stronger impacts were more recently brought forward.

The first one addressed the solar chemical mixture with a revised set of heavy element abundances published by Asplund and collaborators (Asplund et al. 2004, 2005). The abundant heavy elements (C, N, O) saw a strong decrease in their abundances and the metallicity of the Sun was thus reduced by about 30%. Using these new results led to strong disagreements between SSMs and helioseismology (Serenelli et al. 2009). Further revision of the spectroscopic determinations (Asplund et al. 2009, hereafter AGSS09) led the metallicity to increase again slightly, but were insufficient to restore the agreement with helioseismology. These discrepancies were thought to originate in additional physical processes acting in the solar radiative zone (Kumar et al. 1999; Castro et al. 2007), but none of these attempts provided a clear and decisive answer to the issue. Simultaneously, other studies used seismology to estimate the solar metallicity. Some confirmed the GS98 values (Basu & Antia 2006), while others agreed with the AGSS09 values (Vorontsov et al. 2014), illustrating the current stalemate regarding this question.

thumbnail Fig. 1

Frequency ratios r02 and r13 for the Sun and two SSMs. The observations are given with their error bars in blue, while the dashed green line shows the results for an SSM built using the OPAL opacities, the Free equation of state and the AGSS09 abundance tables. The dashed orange line shows the values of these ratios for an SSM built with GN93 abundances, the OPAL opacities and the Free equation of state. The solid green and orange lines show the ratio values when using the OPLIB opacities instead of the OPAL opacities for the AGSS09 and GN93 abundances respectively.

The second important change was the revision of the stellar material opacity. The solar problem has been linked to the opacity at the base of the convective envelope and a process inducing a local increase of the opacity has recurrently been proposed as the solution to the controversy (Basu & Antia 2008). Until recently, the most commonly used opacities were the OPAL opacity tables (Iglesias & Rogers 1996) but an underestimate of the opacity in more massive stars was convincingly revealed by different studies (Salmon et al. 2012; Cugier 2012; Turck-Chièze et al. 2013). These findings initiated both innovative measurements with high-energy laser facilities and numerical computational efforts to improve theoretical calculations. The first experimental results for iron revealed an important discrepancy with theoretical expectations (Bailey et al. 2015). In parallel, two new sets of theoretical opacities were developed, one dedicated to the Sun from the OPAS consortium (Mondet et al. 2015) and the other covering the wide range of stellar conditions, by the Los Alamos National Laboratory (Colgan et al. 2016, hereafter OPLIB opacities), which could become commonly used in stellar models.

The solar issue impacts astrophysics as a whole since the “metallicity scale”, used to relate spectroscopic observations to the metallicity of stellar models, takes the Sun as its reference. To this day, the so-called “solar metallicity problem” remains a tedious issue which is not only linked to the metallicity, but to the whole micro- and macrophysical representation of the stellar structure. Indeed, asteroseismic results have already shown that our depiction of transport processes in stellar models is imperfect (e.g. Mosser et al. 2012; Deheuvels et al. 2014). Owing to the quality of the solar data, the Sun is still our best laboratory for testing the ingredients of stellar models. Consequently, providing new seismic diagnostics allowing a more in-depth probe of the solar structure is crucial. With this study, we provide a new diagnostic by performing structural inversions of an entropy proxy. The sensitivity of this inversion to the stratification just below the convective zone paves the way for a re-analysis of the importance of additional physical processes required in the description of the solar structure. In the following section, we show how our diagnostic sheds new light on the solar structure problem.

thumbnail Fig. 2

Comparison of the entropy profile between AGSS09 and GN93 SSMs. The red and blue crosses illustrate inversion results for the profile of the entropy proxy in the Sun for the GN93 SSM built with the OPAL and OPLIB opacities. The green and orange circles illustrate the effects of changing from the old OPAL opacities to the recent OPLIB opacities in AGSS09 SSMs.

2. Inversion of the solar entropy: a new seismic diagnostic

2.1. Inversion for standard solar models

The models considered in this study are SSMs, built with the Liège stellar evolution code (CLES; Scuflaire et al. 2008a). The frequencies were computed with the Liège oscillation code (LOSC; Scuflaire et al. 2008b). All models presented in this paper are computed using the Free equation of state (Irwin 2012), and either the OPAL or OPLIB opacity tables. In order to fully estimate the effects of a change in the heavy element abundances, we adopted two extreme mixtures, namely GN93 and AGSS09. The structural kernels and the inversions were computed with an adapted version of the InversionKit software (Reese et al. 2012) using the SOLA technique (Pijpers & Thompson 1994). We used the same solar seismic dataset as in Basu et al. (2009) and followed their definitions of the error bars for the inversion. We followed Rabello-Soares et al. (1999) to calibrate the free parameters of the SOLA technique and deal with the surface effects contributions.

Results of sound speed inversions for the new OPLIB opacities are shown in Guzik et al. (2015) and illustrate that SSMs built using the AGSS09 abundances display a slightly deeper convective envelope and slightly better agreement. However, these improvements are mitigated by a larger discrepancy with the helium abundance in the convective envelope found at 0.23. This reduction is also observed for GN93 models which now display a value of 0.24. This trend results from the lower values of the OPLIB opacities in most of the radiative region, which leads to calibrated SSMs with lower initial helium abundances. The improvements of the sound speed profile and the position of the base of the convective envelope do not result from an overall increase of the opacity, but from a steepening of its derivatives, which in turn leads to a steepening of the temperature gradient below the convective zone.

The issue becomes more intricate when the ratios of the small frequency separation to the large frequency separation are analysed. These ratios, denoted r02 and r13, are used to probe the solar core conditions (Roxburgh & Vorontsov 2003), and show a clear preference for the latest AGSS09 mixture when the OPLIB opacities are used, as is illustrated in Fig. 1. This constitutes a clear change from the previous SSMs with the OPAL opacities, which showed better agreement with higher metallicity abundances, such as the GN93 or GS98 tables (Chaplin et al. 2007). Hence, the situation is quite confusing since the sound speed inversion seems to favour the GN93 mixture when using the OPLIB opacities while the frequency ratios agree more closely with low metallicity models using the same opacity tables. It seems fair to admit that no clear solution emerges from classical helioseismic diagnostics. To shed new light on the solar problem, we propose a new seismic diagnostic consisting of inverting a solar entropy proxy, defined as S5/3=Pρ5/3\hbox{$S_{5/3}=\frac{P}{\rho^{5/3}}$}, where P is the pressure and ρ the density, which reproduces the behaviour of the entropy of the solar plasma. The kernels used are thus those of the (S5/31) pair, where Γ1=(lnPlnρ)S\hbox{$\Gamma_{1}=\left( \frac{\partial \ln P}{\partial \ln \rho}\right)_{S}$} is the adiabatic exponent. The constraining nature of this proxy originates in the plateau that it forms in convective regions. This plateau is due to the high efficiency of convection in the deep layers of the solar envelope where this phenomenon operates nearly adiabatically. In turn, the height of the plateau is a direct marker of the way we model the radiative zones of the Sun. In the layers below the convective envelope, the stratification is very sensitive to both opacity and chemical abundances. Consequently, a change in opacity, whatever its origin, or a variation of the abundances will impact both the temperature and the mean molecular weight gradients and thus the height of the plateau in a given solar model. Testing this height through seismic inversions offers a straightforward diagnostic, complementary to that of the sound speed inversions. Moreover, non-standard processes may also change the height of the plateau, making this diagnostic a very sensitive probe of the layers just below the convective envelope, which are precisely the ones where discrepancies are the largest and where the physical hypotheses of the SSMs are the most uncertain.

Inversion results of the entropy proxy profile are given in Fig. 2 for solar models using either the former OPAL or the new OPLIB opacities. The orange and green circles illustrate the results for the AGSS09 abundances, while the blue and red crosses illustrate the results for the GN93 abundances. We notice that the plateau of the entropy proxy is shifted by about 2% due to the opacity changes between the OPAL and OPLIB opacity tables. While the agreement between the Sun and the GN93 SSMs is still of the order of 0.7%, which is quite good, the sign of the differences in the plateau has critical implications. A positive difference between the Sun and the GN93 model built using OPLIB tables means that the entropy plateau in the model is too low. If the GN93 abundances were to be reconciled with the entropy profile of the Sun, it would necessarily require some change inducing a less steep temperature gradient in order to raise the entropy plateau to the solar value. This appears to be in contradiction with the experimental results of Bailey et al. (2015) for iron in the physical conditions present at the base of the envelope, which would lead to a strong steepening of the temperature gradient in this region. Furthermore, theoretical calculations of iron spectral opacity in these conditions are still a matter of debate and could also change in the future (Iglesias & Hansen 2017; Nahar & Pradhan 2016; Blancard et al. 2016).

In contrast, the models built using the AGSS09 abundances and the OPLIB opacities are in better agreement than the GN93 OPLIB models in most of the radiative regions of the solar structure and still show negative differences in the convective envelope. These negative differences mean that a further steepening of the temperature gradient below the convective zone could improve the agreement with the Sun. The entropy inversion being very sensitive to the layers right below the convective zone, it could efficiently constrain non-standard processes. Indeed, additional mechanisms would alter both temperature and mean molecular weight gradients and the changes would be clearly seen in the variations in the height of the plateau.

In addition to the OPLIB opacities, we have also tested the OPAS opacity tables, which have been optimized for the base of the convective zone (Mondet et al. 2015). However, these tables do not cover the full solar conditions and are only available for the AGSS09 abundances, restricting their potential for comparisons using various physical ingredients. Standard solar models built using these tables lead to slightly larger discrepancies with the Sun.

We also tested the dependency of our proxy on the equation of state by using the OPAL 2005 equation of state instead of the Free equation of state. They induce differences that are around ten times smaller than those due to the change in opacities or abundances in the models. Therefore, most of the changes in the entropy plateau are to be expected from non-standard processes or updates in the opacity tables.

2.2. Analysis of the seismic diagnostic of the entropy proxy

Our entropy proxy, denoted S5/3=Pρ5/3\hbox{$S_{5/3}=\frac{P}{\rho^{5/3}}$}, comes from the Sackur-Tetrode equation for the entropy of a mono-atomic non-degenerate ideal gas, which reads S=3kB2(μmuln(Pρ5/3)+f(μ)),\begin{eqnarray} S=\frac{3k_{\rm B}}{2}\left(\mu m_{\rm u} \ln \left( \frac{P}{\rho^{5/3}} \right) + f(\mu)\right), \end{eqnarray}(1)with kB the Boltzmann constant, μ the mean molecular weight, mu the atomic mass unit, P the local pressure, ρ the density, and f(μ) a function that only depends on the mean molecular weight and physical constants. The most striking advantage of this proxy is its unambiguous behaviour towards opacity changes just below the convection zone. Indeed, if the derivative of the natural logarithm of S5/3 is taken with respect to the natural logarithm of P for an ideal gas, we obtain dlnS5/3dlnP=-23+53(dlnTdlnPdlnμdlnP)·\begin{eqnarray} \frac{{\rm d} \ln S_{5/3}}{{\rm d} \ln P}=\frac{-2}{3}+\frac{5}{3}\left(\frac{{\rm d} \ln T}{{\rm d} \ln P} - \frac{{\rm d} \ln \mu}{{\rm d} \ln P} \right) \cdot \end{eqnarray}(2)Now, for a given energy flux, an increased opacity below the convection zone induces a steeper temperature gradient against pressure. This in turn will increase the logarithmic derivative of S5/3 and will bring it closer to 0 since dlnS5/3dlnP\hbox{$\frac{{\rm d} \ln S_{5/3}}{{\rm d} \ln P}$} is negative. Therefore, the increase in entropy versus the radius is thus smaller as the pressure decreases and the height of the plateau is accordingly reduced with the steepening of the temperature gradient just below the convective envelope. This effect is illustrated in Fig. 2, where the near 2% shift results mainly from changing the opacity tables used in the SSMs. The steeper temperature gradient is a consequence of steeper dependence of the OPLIB opacities with temperature (see Colgan et al. 2016). The effect of a localized ad hoc opacity increase on this indicator has been observed in all test cases involving the past OPAL and the latest OPLIB and OPAS opacities. These tests on the structural models confirm the trends we have discussed here.

2.3. Additional tests of the inversion techniques

We also performed further checks of the quality of the averaging kernels for the SOLA method. We illustrate in Fig. 3 the averaging kernels (Pijpers & Thompson 1994) of the SOLA method for various positions inside the Sun. The figure clearly shows that the target function in green is nicely reproduced at every depth, although some inaccuracies are present below 0.1 solar radii. This is expected since we lack very low degree and radial order modes able to probe the deepest layer of the solar structure efficiently.

thumbnail Fig. 3

SOLA averaging kernels and comparison with their target functions at various depths. The green curves show the Gaussian target functions of the SOLA inversion for various depths inside the Sun, while the red dashed curves show the averaging kernels.

3. Conclusion

The change of scenery caused by the use of the OPLIB tables in SSMs points out weaknesses for both high- and low-metallicity abundances tables. An intermediate metallicity value or an increased opacity at the base of the convective envelope could marginally restore the agreement for SSMs, but the discrepancies in helium seem to point out additional mechanisms, some physical ingredients that have to be included in the solar models whatever abundance tables are used. The uncertainties illustrated in this study and the sensitivity of the seismic diagnostic we developed lead us to advocate for a re-opening of the case of potential additional ingredients in helioseismic analyses using constraints such as the lithium abundance and the solar rotational profile in combined studies simultaneously using all the seismic information available. Changes in the physical ingredients of solar and stellar models will affect our determinations of stellar fundamental parameters. It is a necessary step if we want to bring these models to a new level of physical accuracy. For that purpose, seismic inversions of the entropy profile offer unprecedented opportunities to further test the structure of the Sun.

References

  1. Antia, H. M., & Basu, S. 1994, A&AS, 107, 421 [NASA ADS] [Google Scholar]
  2. Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C., & Kiselman, D. 2004, A&A, 417, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C., & Blomme, R. 2005, A&AS, 431, 693 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Bahcall, J. N., & Peña-Garay, C. 2004, New J. Phys., 6, 63 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bailey, J. E., Nagayama, T., Loisel, G. P., et al. 2015, Nature, 517, 3 [Google Scholar]
  6. Basu, S., & Antia, H. M. 1995, MNRAS, 276, 1402 [NASA ADS] [Google Scholar]
  7. Basu, S., & Antia, H. M. 2006, in SOHO-17, 10 Years of SOHO and Beyond, ESA Spec. Publ., 617, 36.1 [Google Scholar]
  8. Basu, S., & Antia, H. M. 2008, Phys. Rep., 457, 217 [NASA ADS] [CrossRef] [Google Scholar]
  9. Basu, S., Chaplin, W. J., Elsworth, Y., New, R., & Serenelli, A. M. 2009, ApJ, 699, 1403 [NASA ADS] [CrossRef] [Google Scholar]
  10. Blancard, C., Colgan, J., Cossé, P., et al. 2016, Phys. Rev. Lett., 117, 249501 [NASA ADS] [CrossRef] [Google Scholar]
  11. Castro, M., Vauclair, S., & Richard, O. 2007, A&AS, 463, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Chaplin, W. J., Serenelli, A. M., Basu, S., et al. 2007, ApJ, 670, 872 [NASA ADS] [CrossRef] [Google Scholar]
  13. Colgan, J., Kilcrease, D. P., Magee, N. H., et al. 2016, ApJ, 817, 116 [NASA ADS] [CrossRef] [Google Scholar]
  14. Cugier, H. 2012, A&A, 547, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. Vangioni-Flam, & M. Casse, 15 [Google Scholar]
  17. Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [NASA ADS] [CrossRef] [Google Scholar]
  18. Guzik, J. A., Fontes, C., Walczak, P., Wood, S. R., & Mussack, K. 2015, IAU General Assembly, 22, 2253707 [NASA ADS] [Google Scholar]
  19. Iglesias, C. A., & Hansen, S. B. 2017, ApJ, 835, 284 [NASA ADS] [CrossRef] [Google Scholar]
  20. Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
  21. Irwin, A. W. 2012, Astrophysics Source Code Library [record ascl: 1211.002] [Google Scholar]
  22. Kosovichev, A. G., & Fedorova, A. V. 1991, Sov. Astron., 35, 507 [NASA ADS] [Google Scholar]
  23. Kosovichev, A. G., Schou, J., Scherrer, P. H., et al. 1997, Sol. Phys., 170, 43 [NASA ADS] [CrossRef] [Google Scholar]
  24. Kumar, P., Talon, S., & Zahn, J.-P. 1999, ApJ, 520, 859 [NASA ADS] [CrossRef] [Google Scholar]
  25. Mondet, G., Blancard, C., Cossé, P., & Faussurier, G. 2015, ApJS, 220, 2 [NASA ADS] [CrossRef] [Google Scholar]
  26. Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Nahar, S. N., & Pradhan, A. K. 2016, Phys. Rev. Lett., 116, 235003 [Google Scholar]
  28. Pijpers, F. P., & Thompson, M. J. 1994, A&A, 281, 231 [NASA ADS] [Google Scholar]
  29. Rabello-Soares, M. C., Basu, S., & Christensen-Dalsgaard, J. 1999, MNRAS, 309, 35 [NASA ADS] [CrossRef] [Google Scholar]
  30. 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]
  31. Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [Google Scholar]
  32. Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Salmon, S., Montalbán, J., Morel, T., et al. 2012, MNRAS, 422, 3460 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  34. Scuflaire, R., Théado, S., Montalbán, J., et al. 2008a, Ap&SS, 316, 83 [NASA ADS] [CrossRef] [Google Scholar]
  35. Scuflaire, R., Montalbán, J., Théado, S., et al. 2008b, Ap&SS, 316, 149 [NASA ADS] [CrossRef] [Google Scholar]
  36. Serenelli, A. M., Basu, S., Ferguson, J. W., & Asplund, M. 2009, ApJ, 705, L123 [NASA ADS] [CrossRef] [Google Scholar]
  37. Turck-Chièze, S., Gilles, D., Le Pennec, M., et al. 2013, High Energy Density Physics, 9, 473 [NASA ADS] [CrossRef] [Google Scholar]
  38. Vorontsov, S. V., Baturin, V. A., Ayukov, S. V., & Gryaznov, V. K. 2014, MNRAS, 441, 3296 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Frequency ratios r02 and r13 for the Sun and two SSMs. The observations are given with their error bars in blue, while the dashed green line shows the results for an SSM built using the OPAL opacities, the Free equation of state and the AGSS09 abundance tables. The dashed orange line shows the values of these ratios for an SSM built with GN93 abundances, the OPAL opacities and the Free equation of state. The solid green and orange lines show the ratio values when using the OPLIB opacities instead of the OPAL opacities for the AGSS09 and GN93 abundances respectively.

In the text
thumbnail Fig. 2

Comparison of the entropy profile between AGSS09 and GN93 SSMs. The red and blue crosses illustrate inversion results for the profile of the entropy proxy in the Sun for the GN93 SSM built with the OPAL and OPLIB opacities. The green and orange circles illustrate the effects of changing from the old OPAL opacities to the recent OPLIB opacities in AGSS09 SSMs.

In the text
thumbnail Fig. 3

SOLA averaging kernels and comparison with their target functions at various depths. The green curves show the Gaussian target functions of the SOLA inversion for various depths inside the Sun, while the red dashed curves show the averaging kernels.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.