Key words Neutrals in plasma, partially ionized plasma, pseudopotential model PACS 52.25.YaIn the framework of the BBGKY hierarchy the main features of interparticle interactions in a partially ionized hydrogen plasma are extensively studied. The theory developed is quite analogous to the Debye-Hückel approximation of a fully ionized plasma.Of interest in the theory of modern plasma physics is the study of thermodynamic properties of partially ionized plasmas. It is evident that the ionization degree is determined by the direct competition of the thermal ionization process and the opposite process of recombimation. In the thermodynamic equilibrium, though, the ionization degree does not depend on details of those processes and can principally be calculated from the pure thermodynamic point of view. The basic procedure here is to assume the free energy and, then, minimize it [1,2].We propose by this a somewhat different approach to the problem that is aimed at studying correlation phenomena in a partially ionized hydrogen plasma consisting of neutrals (atoms), protons and electrons that are referred to in the following by the subscripts n, p and e, respectively. We describe the system of interest by the Wigner-Seitz radiusand, then, introduce the number density parameterHere n = n n + n p is the total amount of protons in the system, a B =h 2 /m e e 2 signifies the first Bohr radius,h is the Planck constant, m e and −e refer to the electron mass and charge, respectively.To characterize the importance of interparticle interactions the coupling parameter is introduced aswhere T is the plasma temperature and k B designates the Boltzmann constant. Fixing two parameters r s and Γ does not determine completely the composition of the plasma medium and, thus, the following Saha equation is added n e n p n n = 2πm e k B T h 2 3/2with the ionization energy I = m e e 4 /2h 2 of the ground state of the hydrogen atom.