Transparency and stability of low density stellar plasma related to Boltzmann statistics and to thermal bremsstrahlung.

The stability of low density stellar plasma is analyzed for a star with a spherical symmetry which is in equilibrium between the gravitational attractive forces and the repulsive pressure forces of an ideal electron gas where the analysis is developed by the use of Boltzmann statistics. Fundamental results are obtained for the radius and total mass of such star and its gravitational forces are large due to the extreme large volume. The absorption and emission of radiation for extremely low density star plasmas is very small over the entire electro-magnetic spectrum.

transparency and stability of low density stellar plasma related to Boltzmann statistics and to thermal bremsstrahlung Y. Ben-Aryeh the stability of low density stellar plasma is analyzed for a star with a spherical symmetry which is in equilibrium between the gravitational attractive forces and the repulsive pressure forces of an ideal electron gas where the analysis is developed by the use of Boltzmann statistics. fundamental results are obtained for the radius and total mass of such star and its gravitational forces are large due to the extreme large volume. the absorption and emission of radiation for extremely low density star plasmas is very small over the entire electro-magnetic spectrum.
In the present work we are interested in the analysis of low density ionized stellar plasma with relatively high temperature for which the Boltzmann statistics is valid. We assume that the plasma behaves as a perfect gas characterized by the fact that the average interaction of the particles is smaller than the thermal energy 1 . For a perfect gas i.e., when the plasma is sufficiently rarefied, with electron density n e (cm −3 ) at temperature T, the electrons pressure P e is given by where k B is Boltzmann constant, T the absolute temperature and n e is the number of ionized-electrons per unit volume. Completely ionized plasma behaves as a perfect gas at low densities. In fact, plasma which is partly ionized behaves as a perfect gas at densities lower than a critical value. This critical value can be found by comparing the Coulomb interaction energy  E e r / In the present work we study the stability of a star under perfect gas conditions where the star stability is produced by the balance between the gravitational forces and the pressure produced by perfect gas conditions, and find the impact of such model on the properties of low density plasmas. While the star stability was studied under high electron densities using Fermi-Dirac statistics 2-6 we study here the problem for low densities with relatively high temperature using Boltzmann statistics. . We should take into account that the present use of Boltzmann statistics for the average velocity of the electrons and the radiation pressure of Eq. (1) are only approximate values which are used for estimating the order of magnitudes for the radius and mass of the low density star. In a more accurate description one needs to take into account the thermal escape of particles (which are in the high velocity region of the Boltzmann distribution) from the low density star. Such effect is similar to that obtained in evaporative cooling and has been treated by using a truncated Boltzmann distribution 12 . In the present work I do not consider this effect but I hope to study it in future work.  In the present work we treat thermal bremsstrahlung effects influencing the absorption and and emission of electromagnetic radiation by low density plasma. We find that such radiative processes are proportional to the products of the electrons and ions densities so that they become very small for very low density plasmas.
For simplicity of discussion we treat a star with spherical symmetry where its stability is related to gravitational forces without any other external perturbations e.g. magnetic fields etc. Magnetic fields might be important for high density stars 13 . But as shown quite long ago by Van Vleck 14 , following Bohr-Van-Leeuven theorem, when classical mechanics, including Boltzmann statistics and free electrons models are applied, the thermal averaged value of the magnetic fields vanishes. Such local classical kinetic equilibrium (LCKE) 15 is valid under the conditions: 1) Quantum effects are negligible in the plasma. 2) None of the distribution functions are time dependent. 3) The distribution function is Maxwellian. As such conditions are assumed in the present analysis we neglect magnetic fields effects. The stability conditions derived in the present work, for low density plasma, imply that as the plasma density is decreased the star radius is increased. As the gravitational forces are proportional to the star volume i.e. to R 3 , where R is the star radius, gravitational forces are significant also for low density plasmas.
Dark matter is thought in previous works to be non-baryonic, possibly being composed of some as-yet undiscovered particles [16][17][18] . There is not any evidence to the existence of such non-baryonic matter and its nature referred in the literature is obscure. In the preset work I discuss the possibility that the present analysis describes dark matter. Usually stars are observed either by radiation emitted by them or by reflected and transmittance of light incident on them from other sources. It is interesting to note by following the present analysis that stellar plasmas with very low densities do not emit any significant radiation and their transmittance is nearly equal 1 (reflectance nearly zero) so that they have the dark matter properties but with baryonic components. By following the conditions for stability of low density star plasma developed in the present work we find that for such star with spherical symmetry, strong gravitational effects occur due to the extremely large volume of such star. Therefore the existence of strong gravitational forces does not imply that such matter is non-baryonic especially as there is not yet any experimental evidence that such non-baryonic matter exists. Later in the discussion I describe certain properties of dark matter related to the present analysis and compare the present analysis with some other works describing dark matter.

the Stability of Stellar plasmas With Low Densities Under Boltzmann Statistics
For simplicity of discussion we treat a star with spherical symmetry where its stability is related to gravitational forces without any other external perturbations. We simplify the treatment by assuming one ionic component with atomic number Z but the analysis can easily be generalized if we have more ionic components.
The gravitational force g at a distance r from the star center is due entirely to the mass M r interior to this distance: where G is the constant of gravitation. Assuming that φ is the gravitational potential then According to the hydrostatic equation ρ = dP g dr (4) where P is the star plasma pressure and ρ the star mass density, both are functions of distance r from the star center. Inserting Eq. (3) into Eq. (4), we get: This equation describes the decrease of the plasma pressure as we move to larger values of r balancing the attractive gravitational forces.
Assuming that the star plasma behaves as an ideal gas then the pressure P is given by Assuming also that the gradient of temperature is small relative to the gradient of the electron density n e (isothermal process) then we get: If there are κ nucleons for each ionized electron, then the mass density is given approximately by  (7) and assuming that at the star center the potential φ vanishes and its ionized electron density is n 0 we get: Here φ r ( ) is changing from zero to large values, as function of the distance r from the star center. Our aim in the following analysis is to find the change of n r ( ) e as function of the distance r from the star center. In these calculations we assume that n 0 is the density of ionized electrons in the star center which is taken as an experimental parameter. In order to find the dependence of the electron density n r ( ) e on the distance r from the star center we need to take into account the Poisson's equation for the potential ϕ which for a star with spherical symmetry has the form: here ρ r ( ) is proportional to the ionized electron density n r ( ) e , as given by Eq. (8), and G is the gravitational constant. Substituting the relation ρ κ = n m e N , from Eq. (8) into Eq. (10), and by using relation (9) for n r ( ) We note that on the right side of Eq. (11) appears an exponential function of the potential φ r ( ) with a very small coefficient given by . One should notice that the coefficient before the exponential is also a very small number for low electron densities. It is difficult to get analytical explicit solutions to this equation, as series expansion of this exponential function converges very slowly for large values of r. In the following analysis we transform Eq. (11) to differential equation of n r ( ) e . It will give after some calculations the change of n r ( ) e , from its initial value = = n r n ( 0) e 0 at the star center (taken as an experimental parameter) to smaller values as a function of the distance r from the star center.
According to Eq. (9) we get Multiplying Eq. (13) by n r ( ) e and rearranging its terms we get: We divide Eq. (15), for θ r ( ) by ξ. Then, this equation can be written as: Here m N is the mass of nucleon and κ is the ratio between the density of nucleons and that of electrons (e.g. for completely ionized Hydrogen plasma κ = 2). By neglecting the third term on the left side of Eq. (17) we obtain a differential equation which was solved numerically and was given as the Lane-Emden equation with = n 2 19 . It shows the change of the value of (where ξ in the present conditions is extremely small number) we find that the radius R of the star stretches over a very long distance which is of order Here we substituted the value of the star radius according to Eq. (18) and inserted an averaged value for the electron density n e which is proportional to its density n o in the star center. In Eq. (18) and (19) we can plug in CGS unit for T and n 0 , and the results for the radius and mass of the star are obtained, respectively. We find that the star mass is increasing with temperature as it is proportional to T Kelvin ( /(1 )) 3/2 . We find also the result that the star mass is inversely proportional to the square root of the density of the electrons in its center due to increase in its volume.
Alternative derivation for the star radius. We give here an alternative derivation for the star radius in order to support our results using a different approximation. We start from Eq. (11), where we substitute, in this equation the following definitions: and express the potential function as a series expansion. Then this equation can be written as: The solution for φ y ( ) can be estimated by using the following series expansion: Inserting Eq. (22) into Eq. (21), and performing the derivatives of φ y ( ) we get: We obtain the coefficients  a b c d , , , , by equating the sum of terms with the same exponent of y. We get: Since z is a small number we get for the solution of φ y ( ): the density of the star n r ( ) e is reduced to n e 1 0 so that the radius R of this star can be estimated as So this estimation for the star radius is larger by factor .  6 2 45 relative to the estimation made in Eq. (18) but both estimations lead to the same orders of magnitudes. optical properties of low density stellar plasma related to thermal bremsstrahlung. Radiative processes in low density plasma at high temperatures are obtained mainly by bremsstrahlung effects representing the acceleration of a charge in the Coulomb field of another charge (usually the acceleration of electrons by the nuclear charge). In the present work we are interested in thermal bremsstrahlung: that is we average the single speed expressions over a thermal distribution of speeds. A full analysis of this process requires a quantum treatment but in some regimes we can use a classical analysis and the quantum effect result as corrections (Gaunt factors 11 ) to the classical formulas. Following such approach the total power per unit volume emitted by spontaneous thermal emission integrated over frequencies is given in CGS units as 7 Thermal free-free bremsstrahlung power absorption per unit frequency f and unit volume, is given as 7 f f 3 where α f is the free-free absorption coefficient and B T ( ) f is the mean energy density at temperature T given by Planck law per unit frequency The absorption coefficient is evaluated as where g f is a small number under the condition  ≤ 1 ze v 2 We notice that both Eqs. (28) and (32) are proportional to the product n n e i so that for very low density plasmas the bremsstrahlung effects are very small. On the other hand by following the stability analysis made in the previous section we find that such low density stars have extremely large radius so that gravitational effects are significant. Although we made the stability analysis for a star with spherical symmetry we predict that low density plasmas with different structures will have also very large volumes.

Summary and Discussion
The present results can be summarized as follows: In low density plasmas the balance between the ideal gas radiation pressure and gravitational forces leads to a large star radius so that although the mass density is quite small the star volume is extremely large leading to a large star mass with strong gravitational forces. We showed also that the absorption and/or emission effects for such plasms can be treated by theories about thermal Bremsstrahlung. As these effects are proportional to the product of the density of electrons n e times the density of ionized ions n i the plasma becomes transparent for very low densities. It seems that the present analysis describes a new kind of dark matter. This idea raises fundamental questions and in the present discussion I give answers to them: (a) In astronomy, luminosity is defined as the total amount of electromagnetic energy emitted per unit time by a star, galaxy, or other astronomical object. One can use Hertzsprung-Russel diagram (HR diagram) for (2019) 9:20384 | https://doi.org/10.1038/s41598-019-56784-2 www.nature.com/scientificreports www.nature.com/scientificreports/ plotting the temperature of stars against their luminosity (see e.g. 20 and references included). The question arises what is the relation between such diagrams and the present analysis? Our answer is that the HR diagram is related to Stephan-Boltzmann law giving the total energy being emitted per unit time and per unit area at all wavelengths of black body radiation given as σT 4 where σ is known as the Stephan-Boltzmann constant and T is the temperature (in Kelvins). For getting the illumination one needs to multiply the electromagnetic energy emitted by Stephan-Boltzmann law by an effective radiation emitting area. This method is applied to an idealized object which is perfectly opaque and non-reflecting and for dense plasmas such assumption is approximately correct. For low density plasmas the emission of radiation depends on the emissivity which is the ratio between the radiation emitted by a certain body and that of black body radiation. For electromagnetic energy incident on plasma conservation of energy implies the relation: Incident energy = Absorption energy + Transmitted energy + Reflected energy. For any material which is in thermal equilibrium assuming that Boltzmann statistics \is not disturbed by the incident radiation then the emissivity of the plasma should be equal to its absorptivity. For low density plasma (or any other material) in which the absorption vanishes the emissivity tends also to zero so that Stephan-Boltzmann law and consequently HR diagram does not apply for such case. (b) Non-luminous matter which introduces strong gravitational waves (known as 'dark matter') was inferred already in 1937 by Swiss astronomer Fritz Zwicky 21 . The existence of this matter was confirmed later by Rubin, Thonnard & Ford 22 . Evidence for the existence of dark matter first came from the fact that luminous matter cannot account for velocities observed in rotation curves implying a significant invisible component. This fact can be related to the present low density plasma in the following way: Considering a small mass density ρ rotating in circular motion around a big spherical mass M with a radius R the centrifugal force per unit volume is given by ρv R / rot 2 where v rot is the rotational velocity. For extremely large distance R this force decays due to the 1/R dependence but the ideal low density plasma pressure given by Eq. (1) remains to be approximately constant over extremely large radial distances. For very large distance from the star center the ideal gas pressure becomes the dominant repulsive force. In conclusion we expect the luminous matter to be embedded in low density plasma which extends over extremely large volumes. There is a strong similarity between the present work and Dark Matter Halos 23 . A dark Halo is the inferred Halo of invisible material that permeates and surrounds galaxies. The difference between the present model and that of dark Halos includes a difference in the dependence on physical parameters (In addition to difference in the physical models). For dark Halos the dependence on physical parameters was obtained by empirical equations in which the density of dark Halos was found to vary as power law in radius 24 . For low density plasma the density remains approximately constant up to a large distance R given approximately by Eq. (18) and beyond such radius it has an exponential decay. The exact way in which it decays has only a secondary effect on the low density plasma physical parameters which depends mainly on the star radius.
In the preset work I discussed the properties of low density plasmas by using approximate analytical equations. I hope in a future work to study numerical solutions to the fundamental Eq. (11) and make comparisons with experimental observations (e.g., those of gravitational lensing 25 ) and with those obtained for dark Halos. Also I hope, in future work, to relate the present model to others cosmological models.