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