Aims. We use Hermite splines to interpolate pressure and its derivatives simultaneously, thereby preserving mathematical relations between the derivatives. The method therefore guarantees that thermodynamic identities are obeyed even between mesh points. In addition, our method enables an estimation of the precision of the interpolation by comparing the Hermite-spline results with those of frequent cubic (B-) spline interpolation. Methods. We have interpolated pressure as a function of temperature and density with quintic Hermite 2D-splines. The Hermite interpolation requires knowledge of pressure and its first and second derivatives at every mesh point. To obtain the partial derivatives at the mesh points, we used tabulated values if given or else thermodynamic equalities, or, if not available, values obtained by differentiating B-splines. Results. The results were obtained with the grid of the SAHA-S equation-of-state (EOS) tables. The maximum $lg P$ difference lies in the range from $10^{-9}$ to $10^{-4}$, and $\Gamma_1$ difference varies from $10^{-9}$ to $10^{-3}$. Specifically, for the points of a solar model, the maximum differences are one order of magnitude smaller than the aforementioned values. The poorest precision is found in the dissociation and ionization regions, occurring at $T \sim 1.5\cdot 10^3 - 10^5$ K. The best precision is achieved at higher temperatures, $T>10^5$ K. To discuss the significance of the interpolation errors we compare them with the corresponding difference between two different equation-of-state formalisms, SAHA-S and OPAL 2005. We find that the interpolation errors of the pressure are a few orders of magnitude less than the differences from between the physical formalisms, which is particularly true for the solar-model points.

The Sun is the most studied of stars and a laboratory of fundamental physics. However, the understanding of our star is stained by the solar modelling problem which can stem from various causes. We combine inversions of sound speed, an entropy proxy and the Ledoux discriminant with the position of the base of the convective zone and the photospheric helium abundance to test combinations of ingredients such as equation of state, abundance and opacity tables. We study the potential of the inversions to constrain ad-hoc opacity modifications and additional mixing in the Sun. We show that they provide constraints on these modifications to the ingredients and that the solar problem likely occurs from various sources and using phase shifts with our approach is the next step to take.

The Sun is the most studied of all stars. It is a reference for all other observed stars and a laboratory of fundamental physics helping us understand processes occuring in conditions irreproducible on Earth. However, our understanding of the Sun is currently stained by the solar modelling problem which can stem from various causes, such as the opacities, the equation of state and the mixing of chemical elements. We combine inversions of sound speed, an entropy proxy and the Ledoux discriminant with constraints such as the position of the base of the convective zone and the photospheric helium abundance. We test various combinations of standard ingredients for solar modelling such as abundance tables, equation of state, formalism for convection and diffusion and opacity tables and study the diagnostic potential of the inversions to constrain ad-hoc modifications of the opacity profile and additional mixing below the convective envelope. Combining inversions provides stringent constraints on the modifications on the models, far beyond what is achievable only from sound speed inversions. We constrain the form and amplitude of the opacity increase required and show that a 15% increase at log T=6.35 provides a significant improvement but that a more global increase of the opacity, within the uncertainties of current tables, coupled with an additional mixing at the bottom of the convective zone gives the best agreement for low metallicity models. We show that high metallicity models do not satisfy all the inversion results. We conclude that the solar problem likely occurs from various small sources, as ingredients such as the equation of state or the formalism of convection can induce small but significant changes and that using phase shift analyses combined with our approach is the next step for a better understanding of the inaccuracies of solar models just below the convective envelope.

We present an example of an interpolation code of the SAHA-S equation of state that has been adapted for use in the stellar evolution code CESAM2k. The aim is to provide the necessary data and numerical procedures for its implementation in a stellar code. A technical problem is the discrepancy between the sets of thermodynamic quantities provided by the SAHA-S equation of state and those necessary in the CESAM2k computations. Moreover, the independent variables in a practical equation of state (like SAHA-S) are temperature and density, whereas for modelling calculations the variables temperature and pressure are preferable. Specifically for the CESAM2k code, some additional quantities and their derivatives must be provided. To provide the bridge between the equation of state and stellar modelling, we prepare auxiliary tables of the quantities that are demanded in CESAM2k. Then we use cubic spline interpolation to provide both smoothness and a good approximation of the necessary derivatives. Using the B-form of spline representation provides us with an efficient algorithm for three-dimensional interpolation. The table of B-spline coefficients provided can be directly used during stellar model calculations together with the module of cubic spline interpolation. This implementation of the SAHA-S equation of state in the CESAM2k stellar structure and evolution code has been tested on a solar model evolved to the present. A comparison with other equations of state is briefly discussed. The choice of a regular net of mesh points for specific primary quantities in the SAHA-S equation of state, together with accurate and consistently smooth tabulated values, provides an effective algorithm of interpolation in modelling calculations. The proposed module of interpolation procedures can be easily adopted in other evolution codes.

Evolution of the 7Li abundance in the convection zone of the Sun during different stages of its life time is considered to explain its low photospheric value in comparison with that of the solar system meteorites. Lithium is intensively and transiently burned in the early stages of evolution (pre-main sequence, pMS) when the radiative core arises, and then the Li abundance only slowly decreases during the main sequence (MS). We study the rates of lithium burning during these two stages. In a model of the Sun, computed ignoring pMS and without extra-convective mixing (overshooting) at the base of the convection zone, the lithium abundance does not decrease significantly during the MS life time of 4.6 Gyr. Analysis of helioseismic inversions together with post-model computations of chemical composition indicates the presence of the overshooting region and restricts its thickness. It is estimated to be approximately half of the local pressure scale height 0.5Hp) which corresponds to 3.8 % of the solar radius. Introducing this extra region does not noticeably deplete lithium during the MS stage. In contrast, at the pMS stage, an overshooting region with a value of approximately 0.18Hp is enough to produce the observed lithium depletion. If we conclude that the dominant lithium burning takes place during the pMS stage, the dispersion of the lithium abundance in solar twins is explained by different physical conditions, primarily during the early stage of evolution before the MS.

At intermediate and high degree $l$, solar p- and f modes can be considered as surface waves. Using variational principle, we derive an integral expression for the group velocities of the surface waves in terms of adiabatic eigenfunctions of normal modes, and address the benefits of using group-velocity measurements as a supplementary diagnostic tool in solar seismology. The principal advantage of using group velocities, when compared with direct analysis of the oscillation frequencies, comes from their smaller sensitivity to the uncertainties in the near-photospheric layers. We address some numerical examples where group velocities are used to reveal inconsistencies between the solar models and the seismic data. Further, we implement the group-velocity measurements to the calibration of the specific entropy, helium abundance Y and heavy-element abundance Z in the adiabatically-stratified part of the solar convective envelope, using different recent versions of the equation of state. The results are in close agreement with our earlier measurements based on more sophisticated analysis of the solar oscillation frequencies (Vorontsov et al. 2013, MNRAS 430, 1636). These results bring further support to the downward revision of the solar heavy-element abundances in recent spectroscopic measurements.

A detailed comparison is carried out between two popular equations of state (EOS), the Mihalas-Hummer-Dappen (MHD) and the OPAL equations of state, which have found widespread use in solar and stellar modeling during the past two decades. They are parts of two independent efforts to recalculate stellar opacities; the international Opacity Project (OP) and the Livermore-based OPAL project. We examine the difference between the two equations of state in a broad sense, over the whole applicable rho-T range, and for three different chemical mixtures. Such a global comparison highlights both their differences and their similarities. We find that omitting a questionable hard-sphere correction, tau, to the Coulomb interaction in the MHD formulation, greatly improves the agreement between the MHD and OPAL EOS. We also find signs of differences that could stem from quantum effects not yet included in the MHD EOS, and differences in the ionization zones that are probably caused by differences in the mechanisms for pressure ionization. Our analysis do not only give a clearer perception of the limitations of each equation of state for astrophysical applications, but also serve as guidance for future work on the physical issues behind the differences. The outcome should be an improvement of both equations of state.