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.

The solar abundance controversy inspires renewed investigations of the basic physics used to develop solar models. Here we examine the correction to the proton-proton reaction rate due to dynamic screening effects. Starting with the dynamic screening energy from the molecular-dynamics simulations of Mao et al., we compute a reaction-rate correction for dynamic screening. We find that, contrary to static screening theory, this dynamic screening does not significantly change the reaction rate from that of the bare Coulomb potential.

In the hot, dense plasma of solar and stellar interiors, the Coulomb interaction is screened by the surrounding plasma. Although the standard Salpeter approximation for static screening is widely accepted and used in stellar modeling, the question of dynamic screening has been revisited. In particular, Shaviv and Shaviv apply the techniques of molecular dynamics to the conditions in the solar core in order to numerically determine the dynamic screening effect. By directly calculating the motion of ions and electrons due to Coulomb interactions, they compute the effect of screening without the mean-field assumption inherent in the Salpeter approximation. Here we reproduce their numerical analysis of the screening energy in the plasma of the solar core and conclude that the effects of dynamic screening are relevant and should be included in the treatment of the plasma, especially in the computation of stellar nuclear reaction rates.

This is the first paper in a series that deals with solar-physics applications of the equation-of-state formalism based on the formulation of the so-called "Feynman-Kac (FK) representation". Here, the FK equation of state is presented and adapted for solar applications. Its domain of validity is assessed. The practical application to the Sun will be dealt with in Paper II. Paper III will extend the current FK formalism to a higher order. Use of the FK equation of state is limited to physical conditions for which more than 90% of helium is ionized. This incudes the inner region of the Sun out to about .98 of the solar radius. Despite this limitation, in the parts of the Sun where it is applicable, the FK equation of state has the power to be more accurate than the equations of state currently used in solar modeling. The FK approach is especially suited to study physical effects such as Coulomb screening, bound states, the onset of recombination of fully ionized species, as well as diffraction and exchange effects. The localizing power of helioseismology allows a test of the FK equation of state. Such a test will be beneficial both for better solar models and for tighter solar constraints of the equation of state.

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.

Using a criterion proposed by Salpeter,and standard solid state physics,we have determined the Debye temperature of a solid in equilibrium with the electron gas surrounding it.The results obtained can have implications in studies of the interstellar medium.

Although 98% of the solar material consists of hydrogen and helium, the remaining chemical elements contribute in a discernible way to the thermodynamic quantities. The contribution of various heavy elements in a set of thermodynamic quantities has been examined. It has emerged that prominent signatures of individual elements exist, and that these effects are greatest in the ionization zones, typically located near the bottom of the convection zone. The main result is that part of the features found here depend strongly both on the given species (atom or ion) and its detailed internal partition function, whereas other features only depend on the presence of the species itself, not on details such as the internal partition function. The latter features are obviously well suited for a helioseismic abundance determination, while the former features present a unique opportunity to use the sun as a laboratory to test the validity of physical theories of partial ionization in a relatively dense and hot plasma. This domain of plasma physics has so far no competition from terrestrial laboratories. [shortened version of original abstract]

The generalized Fermi-Dirac functions and their derivatives are important in evaluating the thermodynamic quantities of partially degenerate electrons in hot dense stellar plasmas. New recursion relations of the generalized Fermi-Dirac functions have been found. An effective numerical method to evaluate the derivatives of the generalized Fermi-Dirac functions up to third order with respect to both degeneracy and temperature is then proposed, following Aparicio. A Fortran program based on this method, together with a sample test case, is provided. Accuracy and domain of reliability of some other, popularly used analytic approximations of the generalized Fermi-Dirac functions for extreme conditions are investigated and compared with our results.

Using a criterion proposed by Salpeter and standard solid-state physics,we have determined conditions for the occurence of the plasma-solid transition.Possible astrophysical applications are discussed.

Using a simple criterion proposed by Kirzhnitz,we have discussed the plasma-solid transition in two simple model cases: the Fermi-Dirac and Bose Einstein gas.Some possibilities of astrophysical applications of our results are also indicated.

The Mihalas-Hummer-Dappen (MHD) equation of state is a part of the Opacity Project (OP), where it mainly provides ionization equilibria and level populations of a large number of astrophysically relevant species. Its basic concept is the idea of perturbed atomic and ionic states. The MHD equation of state was originally restricted to the plasma of stellar envelopes, that is, to relatively moderate densities, which should not exceed rho < 10**-2 g/cm**3. However, helioseismological analysis has demonstrated that this restriction is much too conservative. The principal feature of the original Hummer & Mihalas (1988) paper was a simplified form of the Holtsmark microfield for randomly distributed ions. An improved MHD equation of state (Q-MHD) is introduced. It is based on a more realistic microfield distribution (Hooper 1966, 1968) that includes plasma correlations. Comparison with an alternative post-Holtsmark formalism (APEX) is made, and good agreement is shown. There is a clear signature of the choice of the microfield distribution in the adiabatic index gamma_1, which makes it accessible to present-day helioseismological analysis. But since these thermodynamic effects of the microfield distribution are quite small, it also follows that the approximations chosen in the original MHD equation of state were reasonable. A particular feature of the original MHD papers was an explicit list of the adopted free energy and its first- and second-order analytical derivatives. The corresponding Q-MHD quantities are given in the Appendix.

Misunderstandings have occurred regarding the conclusions of the paper by S. Arndt, W. Dappen and A. Nayfonov 1998, ApJ 498, 349. At occasions, its results were interpreted as if it had shown basic flaws in the general theory of dynamical screening. The aim of this comment is to emphasize in which connection the conclusions of the paper have to be understood in order to avoid misinterpretations.

The difference in the adiabatic gradient gamma_1 between inverted solar data and solar models is analyzed. To obtain deeper insight into the issues of plasma physics, the so-called ``intrinsic'' difference in gamma_1 is extracted, that is, the difference due to the change in the equation of state alone. Our method uses reference models based on two equations of state currently used in solar modeling, the Mihalas-Hummer-Dappen (MHD) equation of state, and the OPAL equation of state (developed at Livermore). Solar oscillation frequencies from the SOI/MDI instrument on board the SOHO spacecraft during its first 144 days in operation are used. Our results confirm the existence of a subtle effect of the excited states in hydrogen that was previously studied only theoretically (Nayfonov & Dappen 1998). The effect stems from internal partition function of hydrogen, as used in the MHD equation of state. Although it is a pure-hydrogen effect, it takes place in somewhat deeper layers of the Sun, where more than 90% of hydrogen is ionized, and where the second ionization zone of helium is located. Therefore, the effect will have to be taken into account in reliable helioseismic determinations of the astrophysically relevant helium-abundance of the solar convection zone.

We present stellar evolutionary models covering the mass range from 0.4 to 1 Msun calculated for metallicities Z=0.020 and 0.001 with the MHD equation of state (Hummer & Mihalas, 1988; Mihalas et al. 1988; D\"appen et al. 1988). A parallel calculation using the OPAL (Rogers et al. 1996) equation of state has been made to demonstrate the adequacy of the MHD equation of state in the range of 1.0 to 0.8 Msun (the lower end of the OPAL tables). Below, down to 0.4 Msun, we have justified the use of the MHD equation of state by theoretical arguments and the findings of Chabrier & Baraffe (1997). We use the radiative opacities by Iglesias & Rogers (1996), completed with the atomic and molecular opacities by Alexander & Fergusson (1994). We follow the evolution from the Hayashi fully convective configuration up to the red giant tip for the most massive stars, and up to an age of 20 Gyr for the less massive ones. We compare our solar-metallicity models with recent models computed by other groups and with observations. The present stellar models complete the set of grids computed with the same up-to-date input physics by the Geneva group [Z=0.020 and 0.001, Schaller et al. (1992), Bernasconi (1996), and Charbonnel et al. (1996); Z=0.008, Schaerer et al. (1992); Z=0.004, Charbonnel et al. (1993); Z=0.040, Schaerer et al. (1993); Z=0.10, Mowlavi et al. (1998); enhanced mass loss rate evolutionary tracks, Meynet et al. (1994)].