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 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.

A reconnecting current layer is a `heart' of a solar flare, because it is a place of magnetic-field energy release. However there are no direct observations of these layers. The aim of our work is to understand why we actually do not directly observe current layers and what we need to do it in the future. The method is based on a simple mathematical model of a super-hot (T ~ 1E8 K) turbulent-current layer (SHTCL) and a model of plasma heating by the layer. The models allow us to study a correspondence between the main characteristics of the layer, such as temperature and dimensions, and the observational features, such as differential and integral emission measure of heated plasma, intensity of spectral lines Fe XXVI (1.78 and 1.51A) and Ni XXVII (1.59 A). This method provides a theoretical basis for determining parameters of the current layer from observations. Observations of SHTCLs are difficult, because the spectral line intensities are faint, but it is theoretically possible in the future. Observations in X-ray range 1.5--1.8 A with high spectral resolution (better than 0.01 A) and high temporal resolution (seconds) are needed. It is also very important to interpret the observations using a multi-temperature approach instead of the usual single or double temperature method.